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1.0  nmtODDCTION 


This  report  describes  results  obtained  during  a  five-year  research 
program  on  the  Fundamental  Processes  in  Partially  Ionized  Plasmas  con¬ 
ducted  in  the  High  Temperature  Gasdynamics  Laboratory  at  Stanford  Uni¬ 
versity.  This  research  was  supported  by  a  grant  from  the  Air  Force 
Office  of  Scientific  Research  (AFOSR-83-0108)  and  was  conducted  under 
the  direction  of  Professors  Charles  H.  Kruger,  Morton  Mitchner,  and 
Sidney  Self.  Three  Ph.D.  candidates  are  currently  conducting  their 
doctoral  research  under  this  program,  while  one  student  completed  his 
Ph.D.  with  support  from  the  program. 

Several  space  power  and  propulsion  systems  of  potential  long-range 
interest  to  the  Air  Force  involve  partially  ionized  plasmas.  Such 
systems  include  MPD  thrusters,  both  open-  and  closed-cycle  MHD  power 
generation,  and  thermionic  energy  conversion.  Although  the  specific 
configurations,  the  exact  operating  conditions,  and  which  of  the  compet¬ 
ing  systems  will  prove  to  be  most  useful  in  the  long  term  remain  to  be 
established,  it  is  important  at  this  time  to  provide  a  broad  fundamental 
research  base  in  support  of  development  activity.  In  particular,  there 
are  a  number  of  key  issues  regarding  the  properties  and  discharge  behav¬ 
ior  of  partially  ionized  plasmas  and  the  interaction  of  discharges  with 
fluid  dynamics  that  need  to  be  understood  before  the  potential  and 
limitations  of  competing  systems  can  be  fully  evaluated.  In  addition, 
it  is  important  that  outstanding  young  applied  scientists  be  educated  in 
these  areas. 

The  present  research  on  partially  ionized  plasmas  is  discipline- 
rather  than  device-oriented  and  is  focused  on  three  major  areas; 

1 .  Recombination  in  molecular  plasmas 

2.  Discharge  effects:  plasma  electrode  interaction 

3.  Interaction  of  discharges  and  fluid  dynamics 

In  addition,  each  of  these  areas  involves  the  development  and  applica¬ 
tion  of  modern  plasma  diagnostic  techniques. 
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These  areas  are  overlapping  and  mutually  supportive.  Thus,  under¬ 
standing  of  plasm  properties  is  important  to  the  study  of  discharges 
and  their  interaction  with  fluid  dynamics.  In  the  same  spirit,  plasma 
diagnostic  methods  are  developed  so  that  they  can  be  applied  in  achiev¬ 
ing  our  research  objectives.  In  areas  1)  and  3),  the  research  is  pri- 
mrily  experimental  in  nature,  with  supporting  theoretical  studies  for 
the  planning  of  the  research  and  interpretation  of  the  data.  Area  2) 
has  so  far  been  restricted  to  theoretical  modeling  and  con5)utations . 

Progress  in  each  of  the  three  research  areas  is  described  in  Section 
2.0.  Publications  and  Presentations  resulting  from  this  work  are  cited 
in  Section  4.0,  and  Section  5.0  lists  the  personnel  who  have  contributed 
to  this  report. 
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2.0  PBOJBCT  DESCRIPTIONS 

Included  in  this  section  are  descriptions  of  progress  in  each  of 
three  project  areas.  Each  project  contains  the  following  subsections; 
(a)  Introduction;  (b)  Research  Objectives;  and  (c)  Status  of  Research. 
Additional  descriptions  may  be  found  in  the  publications  listed  in  Sec¬ 
tion  4.0. 

2.1  Reconibiiiation  in  Molecular  Plasaas 
2.1.1  Introduction 

The  equilibrium  thermodynamic  properties  and  the  quasi-equilibrium 
transport  properties  of  partially  ionized  plasmas  have  been  fairly  well 
worked  out,  at  least  with  regard  to  the  theoretical  aspects,  and  many 
in^ortant  results  are  given  in  the  book  Partially  Ionized  Gases  [1.1]. 
There  rerrain,  however,  a  nuiriber  of  areas  where  plasma  properties  are  not 
now  adequately  understood,  particularly  with  regard  to  nonequilibrium 
plasmas.  Since  a  knowledge  of  plasma  properties  is  fundamental  to  any 
system  description,  it  is  important  that  this  area  receive  appropriate 
attention.  The  electron  density  is  often  the  single  most  inportant 
thermodynamic  property;  in  many  nonequilibrium  situations  it  is  deter¬ 
mined  by  finite  ionization  and  recombination  rates.  These  rates  in  turn 
depend  on  the  constituents  of  the  gas  and  the  form  of  nonequilibrium. 

In  applications  of  interest  to  the  Air  Force  (e.g.,  MHD  generators 
and  space  propulsion)  where  a  high-pressure  (i.e.,  of  the  order  of  an 
atmosphere)  partially-ionized  plasma  flows  past  a  cooled  surface,  recom¬ 
bination  of  electrons  occurs  with  positive  alkali -metal  ions  under 
conditions  for  which  there  have  been  few,  if  any,  reported  measurements 
of  the  recombination  rate.  Typically,  these  conditions  are  for  rela¬ 
tively  low  electron  number  densities  n^  where  the  three-tody  recombina¬ 
tion  process 

e+M'^+e— >M+e,  (1.1) 

1 2 

usually  dominant  at  higher  electron  number  densities  (e.g.,  Ug  ^  10 

-3 

cm  ,  at  T  >  1000  K) ,  no  longer  plays  a  major  role.  (The  recombination 
coefficient  Og  for  this  process  is  defined  by  the  relation  dPg/dt  = 
~°*e”e”i '  denotes  the  ion  number  density  and  t  is  the  tine.) 
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The  goal  of  this  study  is  to  obtain  an  improved  understanding  of 
recombination  processes  under  the  aforementioned  conditions  of  low  n^ 
and  atmospheric  pressures,  by  developing  a  means  for  producing  a  highly 
nonequilibrium  plasma  in  a  mixture  containing  an  alkali  metal  vapor  and 
a  molecular  gas,  and  making  measurements  thereon.  From  an  engineering 
point-of-view,  it  is  the  actual  value  of  the  recombination  rate  which  is 
of  primary  importance,  independent  of  what  the  underlying  kinetic  mecha¬ 
nisms  may  be.  However,  a  knowledge  of  the  mechanisms  is  very’  useful  in 
that  it  enables  one  to  extrapolate  and  predict  values  of  the  recombina¬ 
tion  rate  for  conditions  that  go  beyond  those  for  which  measurements 
exist . 

In  addition  to  measurements  of  the  recombination  rates,  another 
objective  of  our  proposed  work  is  therefore  to  identify  the  operative 
recombination  mechanisms  and  to  make  comparisons  of  our  measurements 
with  existing  theories.  It  has  been  proposed  by  Bates  et  al.  [1*2]  that 


the  three-body  recombination  process 


e  +  +  X,, 


M  +  Xv 


(1.2) 


where  Xj^  denotes  a  molecular  species,  is  important  under  the  conditions 
of  interest.  The  recombination  coefficient  Oy.  for  this  reaction  (de¬ 
fined  by  the  relation  dn^/dt  =  “O^’^e'^i^  been  calculated  by  Bates  et 
al.  [1.2]  for  the  case  of  nitrogen  molecules,  but  no  experiments  have 
been  reported  to  date  that  would  enable  a  comparison  to  be  made  with 
theory . 

At  higher  pressures,  collisional  atomic  processes  dominate  radiative 

processes,  so  that  one  may  employ  the  important  simplifications  = 

n  cx'(T)  and  a-  =  n^uMT).  Here  a’(T)  and  oMT)  are  functions  of  temper- 
e  c  J  J  J  =  J 

ature,  and  nj  denotes  the  number  density  of  the  molecular  species  j. 
The  quantities  Og  and  Oj  are  also  referred  to  (loosely)  as  recombination 
coefficients. 

Experimental  studies  of  the  recombination  process  in  which  the  third 
body  is  an  atom  or  molecule  have  been  carried  out  using  (hydrocarbon) 
flames  [1.3],  where  the  products  of  combustion  contained  several  molecu¬ 
lar  species.  To  extract  information  about  the  individual  molecules  Xj, 


rm 
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the  assumption  was  made  that  the  contributions  from  the  different  mole¬ 
cules  in  the  mixture  were  additive .  Expressed  equivalently,  in  terms  of 
recombination  coefficients,  it  was  assumed  that  the  overall  rate  a 
(which  was  measured,  in  effect)  could  be  written  in  terms  of  the  indi¬ 
vidual  rates  Oj  by  the  formula 

a  =  F  X. a.  , 

•  11 
1 

where  Xj  =  n^/n  denotes  the  mole  fraction  of  species  X j .  Here  nj  is  the 
number  density  of  species  Xj  and  n  is  the  total  number  density. 

The  results  obtained  in  this  fashion  have  been  strongly  criticized 
by  Bates  [1.4],  who  claimed  that  the  recombination  coefficient  for  a 
molecular  mixture  is  not  additive ,  but  that  it  is  given  by  the  formula 


« =  - 11  — ^ — r’  • 

^  p  y  x.B.(p) 

j  ^  ^ 

The  quantities  Bj(p)  are  related  to  the  rates  of  de-excitation  of  M  by 
Xj  from  all  energy  levels  above  and  including  p,  to  all  levels  below  p. 
The  factor  n  is  given  in  terms  of  atomic  constants  and  the  temperature  T 
as  (2irmgkT/h^)^/^  =  KT^'^^.  Values  of  Bj(p)  have  been  calculated  by 
Bates  et  al.  [1.2]  for  the  molecular  species  H2  and  N2« 

For  a  single  species,  one  obtains  from  this  expression  for  a  the 
result 


a . 
1 


n  . 
3 


il 


1-1 


p 


It  is  therefore  apparent  that  the  value  of  a  for  a  mixture  of  molecules 
cannot  be  written,  in  general,  in  terms  of  the  individual  Oj  •  To  be 
able  to  calculate  a  for  a  mixture,  the  factors  Bj(p)  must  be  used. 
However,  before  being  able  to  rely  on  the  calculated  factors  Bj(p),  it 
is  necessary  to  test  the  theory  by  obtaining  data  in  which  only  a  single 
molecular  species  Xj  is  present  in  the  recombining  gas  together  with  the 
al)cali  metal  M.  It  is  clear  that  flame  experiments  are  unable  to  yield 
such  data. 
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2.1.2  Research  Objectives 

1.  To  design  and  construct  an  experimental  system  that  can  produce  a 

recombining  alkali  metal  plasma  under  conditions  where  the  molecular 

three-body  process  dominates. 

2.  To  examine,  in  particular,  the  recombination  of  cesium  ions  in  the 
presence  of  nitrogen. 

3.  To  measure  the  desired  recombination  rate  as  a  function  of  the 

background  gas  temperature. 

4.  To  measure  the  recombination  rate  as  a  function  of  the  partial 

pressures  of  Cs  and  N2. 

5.  To  examine  the  recombination  rate  of  interest  in  a  mixture  of  molec¬ 
ular  species. 

2.1.3  Status  of  Research 

The  remainder  of  this  section  describes  first  theoretical  calcula¬ 
tions  directed  at  analyzing  two  possible  experimental  systems  for  pro¬ 
ducing  the  nonequilibrium  conditions  necessary  for  studying  the  desired 
three-body  molecular  recombination  process.  One  concept  is  based  on  the 
use  of  a  supersonic  nozzle  to  rapidly  expand  a  flowing  high-temperature 
plasma  consisting  of  N2  and  thermally  ionized  Cs.  This  approach  would 
provide  steady-state  conditions.  Measurements  of  n^  at  several  spatial 
loca  'ons  (e.g.,  before  and  after  the  expansion)  would  provide  informa¬ 
tion  for  determining  the  recombination  coefficient.  The  second  concept 
would  employ  a  pulse  of  radiation  to  produce  ionizational  nonequilibrium 
by  photoionization  of  Cs  vapor  confined  with  N2  in  a  heated  test  volume. 
Measurements  of  n^  as  a  function  of  time  following  a  rapid  quenching  of 
the  photoionizing  radiation  would  enable  determination  of  the  recombina¬ 
tion  rate.  The  reasons  for  choosing  the  second  approach  are  presented. 

The  second  part  of  this  section  describes  the  experimental  work 
performed  in  developing  the  necessary  hardware  to  execute  the  pulsed 
photo ionization  experiments.  The  work  is  presented  in  four  main  parts; 
Reaction  Cell,  Gas  Handling,  Light  Sources,  and  Diagnostics. 


A.  Theoretical  Worfc 

A. 1  Supersonic  Expansion  Approach 

A. 1.1  Experimental  Configuration 

An  experimental  system  for  producing  a  nonequilibrium  recombining 
plasma  using  a  supersonic  expansion  approach  is  illustrated  schemati¬ 
cally  in  Fig.  1.2.  An  arc  heater  is  used  to  produce  a  flow  of  high- 
ten^erature  Ar  into  a  large  plenum  chamber.  An  alkali  metal  vapor  such 
as  cesium  or  potassium  is  added  to  the  gas,  together  with  various  molec¬ 
ular  species  introduced  either  singly,  or  in  controlled  combinations. 
The  purpose  of  the  plenum  is  to  permit  a  thorough  mixing  of  the  gases 
and  to  enable  the  mixture  to  come  to  thermodynamic  equilibrium. 

The  gas  mixture  then  flows  through  a  chamber  viewport  where  the 
stagnation  temperature  of  the  gas  can  be  measured  with  thermocouples  and 
by  spectroscopic  methods.  A  convergent-divergent  supersonic  nozzle  fol¬ 
lows,  in  which  the  gas  is  rapidly  cooled,  so  that  electron  recombination 


'  '  \  I 
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Fig.  1.2.  Experimental  system  based  on  supersonic  expansion  approach. 


is  induced.  The  electron  number  density  is  measured  spectroscopically 
or  by  an  electrical  conductivity  diacrnostic  at  the  exit  of  the  nozzle, 
thereby  providing  data  from  which  the  recombination  coefficient  can  be 
calculated.  The  gas  then  passes  through  a  supersonic  diffuser  and  is 
exhausted  to  the  atmosphere . 

A. 1.2  Recombination  Model 

To  design  the  apparatus,  it  is  necessary  to  determine  approximately 
the  optimal  experimental  conditions,  within  the  constraints  imposed  by 
the  operating  limits  of  the  arc  jet.  These  conditions  include  the 
plenum  pressure  and  temperature,  the  number  densities  of  the  potassium 
(which  we  shall  take  here  to  be  the  alkali  metal  vapor)  and  the  molecu¬ 
lar  gas  additives,  and  the  geometry  of  the  converging-diverging  nozzle. 
Since  three-body  recombination  will  occur  in  which  both  the  electrons 
and  the  molecules  serve  as  third  bodies,  it  is  desirable  that  experi¬ 
ments  be  conducted  for  conditions  such  that  the  latter  process  is  domi¬ 
nant.  Shown  in  Fig.  1.3,  are  the  values  of  for  N2/  as  calculated  by 
Bates  et  al.  [1.2],  and  the  values  of  Og  as  given  by  Hinnov  and 
Hirschberg  [1.5].  We  shall  take  ^2  as  the  molecular  third  body  in  the 
discussion  that  follows. 


The  quasi-one-dimensional  steady-state  electron  continuity  equation 
that  governs  the  behavior  of  n^,  is 

dn  ,  , 

e  du  d  .  • 

u  —  +  n  —  +un  —  £nA  =  n  , 

dx  e  dx  e  dx  e 


where  n^  denotes  the  net  electron  production  rate,  and  is  given  by  the 
expression 


n  =  (ns  -  n^)  (no*  +  n  oj  ) 
e  '  K  e-^  ^  e  e  ^2^2 


Here  s  =  s(T)  =  exp(-e^/kT),  where  is  the  ionization  energy  of 

potassium.  In  obtaining  n^,  the  ionization  rates  have  been  expressed  as 
Sg  =  s(T)ag  and  =  s(T)ojj  by  use  of  the  method  of  detailed  balanc¬ 

ing.  The  gas  flow  speed  u  (as  well  as  all  other  gas  properties)  is 
given  in  terms  of  the  cross-sectional  area  of  the  nozzle  A  =  A(x),  bv 
the  well-known  relations  of  gas  dynamics  [1.6]. 
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Fig.  1.3.  Temperature  dependence 
of  the  three-body  recombination 
coefficients  with  either  N2  or 
electrons  serving  as  the  third 
body. 
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Fig.  1.4.  Domain  of  experimental 
conditions  required  for  determining 
the  three-body  recombination  coef¬ 
ficient  with  N,. 


Before  discussing  some  of  the  results  of  integrating  the  continuity 
equation  for  n^Cx),  it  is  useful  to  introduce  two  concepts  that  serve  as 
a  guide  in  selecting  the  initial  condition  at  the  entrance  to  the 
nozzle.  We  denote  the  ratio  of  the  recombination  rate  by  N2  to  that  by 
electrons  as  r  =  n^  /(ngUg),  and  the  ratio  of  the  characteristic  time 
for  recombination  by  N2  to  the  residence  time  of  a  fluid  element  in  the 
nozzle  as  t  e  t^^  /t^^g.  Here  =  (o^  n^  ng)'"'  and  t^.^^  »  L/a,  where  L 

^  A  -  Z  £. 

is  the  nozzle  length  and  a  =  (yRT)  is  the  speed  of  sound. 

We  would  anticipate  that  desirable  upstream  gas  states  would  satisfy 
the  conditions  r  >  1  and  t  ~  1.  (For  t  «  1  the  flow  would  approach  the 
"equilibrium"  limit,  in  which  ionization  and  recombination  rates  are  ap¬ 
proximately  equal.  For  t  »  1  the  flow  would  be  approximately  "frozen," 
where  little  recombination  by  N2  occurs.)  The  condition  on  r  would 
ensure  that  recombination  rate  by  N2  would  exceed  that  by  electrons- 
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From  the  equation  defining  r  we  may  write 

2  2 


and  substituting  for  n  in  the  equation  defining  x,  we  obtain 


V  ■  [(f) 


(»;%. 


200  /t] 


(We  have  used  here  the  approximate  expression  t^^^g  =  0.005/»^,  corre¬ 
sponding  to  a  nozzle  length  of  0.1  m) .  Shown  in  Fig.  1.4  are  the  values 

of  plotted  as  a  function  of  temperature  for  values  of  the  parameter 

2 

(r/i)  equal  to  1  and  10,  and  the  corresponding  values  of  Ug  for  r  =  1 
and  10.  Also  shown  are  constant  partial  pressure  lines  for  N2  and  K 
with  pjj  =6  atm  and  Pj^  =  0.1  atm.  These  values  are  estimates  of  the 
approximate  maximum  values  of  the  nitrogen  and  potassium  partial  pres¬ 
sures  that  can  be  used  with  the  available  experimental  facilities. 

The  estimated  range  of  possible  experimental  values  for  n„  and  the 

^2 

temperature  fall  in  the  region  defined  by  T  ^  2000  K  and  by  the  curves 
labeled  (r/x)  =  1  and  p^^  =  6  atm.  It  appears  that  values  of  (r/x)  up 
to  about  10  can  be  achieved. 

To  obtain  a  more  accurate  picture  than  is  provided  by  the  preceding 
analysis,  a  program  was  written  to  integrate  the  electron  continuity 
equation  numerically.  The  purpose  of  this  program  was  to  estimate  the 
optimal  experimental  conditions  that  should  be  used  for  determining 
a^^(T).  The  variation  of  n^  with  streamwise  location  in  a  Mach  2 
converging-diverging  nozzle  is  shown  in  Figs.  1.5  and  1.6  for  various 
representative  values  of  the  upstream  stagnation  conditions  that  may  be 
achieved  with  the  steady-state  gas  flow  system  presently  available  in 
the  High  Temperature  Gasdynamics  Laboratory.  The  operating  characteris¬ 
tics  of  the  arc  jet  heater  in  that  system  limit  the  N2  partial  pressure 
to  about  four  atmospheres,  or  in  equivalent  terms  n^,  <  10  m  ^  for  a 
plenum  temperature  of  3000  K.  (Figures  1.5a  and  1.6b  are  for  conditions 
which  are  somewhat  more  extreme . ) 


I 


The  four  curves  in  each  figure  show  the  "frozen"  flow,  the  "equilib¬ 
rium"  flow,  the  recombining  flow,  and  the  contribution  to  the  recombina¬ 


tion  by  just  the  electron  three-body  process.  Figures  1.5  and  1.6  are 
for  stagnation  temperatures  of  2000  K  and  3000  K  respectively.  With 


reference  to  Pig.  1.5a,  Fig.  1.5b  shows  the  effect  of  increasing  n^j 


and  Fig.  1.5c  shows  the  effect  of  increasing  n  .  With  reference  to  Pig. 


1.6c,  Fig.  1.6b  shows  the  effect  of  increasing  n^  ,  and  Fig.  1.6c  shows 

"2 


the  effect  of  decreasing  n^.  The  term  "^2  effect"  for  each  figure  is 
the  fractional  decrease  in  n^  attributable  to  three-body  recombination 


with  N2,  that  is,  [ng(a^^  =  0)  -  n^l/ng.  These  results  show  that  in¬ 


creasing  the  values  of  n^  and  n^  are  generally  favorable  and  increasing 


T  is  unfavorable  with  respect  to  increasing  the  fractional  effect  of  N2» 


A. 1.3  Limitations  of  the  Steady-State  Experiment 


It  appears  from  Figs.  1.5  and  1.6  that  in  order  for  the  value  of  n 


at  the  exit  of  the  nozzle  to  depart  significantly  from  the  frozen  value. 


it  is  necessary  to  have  t  <  1,  but  this  requirement  then  precludes 


achieving  the  (desirable)  condition  r  ^  10.  If  decreased  so  as 
to  increase  r,  the  value  of  t  is  also  increased,  thereby  reducing  the 
effects  of  recombination  relative  to  frozen  flow.  This  behavior  may  be 


seen  more  clearly  from  the  relation 


obtained  by  combining  the  expressions  for  r  and  t  so  as  to  eliminate  n 


For  fixed  N2  partial  pressure  and  nozzle  dimensions,  the  right-hand  side 


is  a  function  only  of  the  stagnation  temperature.  Taking  Pj^  =4  atm 


and  tj.gg  =  0.005//T,  one  obtains  r/x  =  5.4  for  T  =  2000  K  and  r/x  =7.2 


for  T  =  3000  K. 


Another  limitation  of  the  rapid  expansion  experiment  stems  from  the 
practial  difficulties  of  providing  optical  access  for  diagnostic  mea¬ 


surements  of  the  expanding  gas  at  more  than  a  few  axial  locations  (e.g., 
the  plenum,  the  exit,  and  perhaps  one  intermediate  position) .  Since  the 
gas  temperature  decreases  by  approximately  1000  K  from  the  upstream  ple¬ 
num  to  the  downstream  exit,  the  accuracy  of  determining  the  temperature 
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dependence  of  ot|jjj^(T)  would  be  limited  by  the  relatively  coarse  resolu¬ 
tion  with  which  the  dependence  on  axial  location  of  n^  could  be  mea¬ 
sured.  The  use  of  the  available  arc  jet  introduces  a  further  possible 
complication  in  that  argon  atoms  will  be  present  in  the  recombining 
plasma  * 

A. 2  Pulsed  Photoionization  Approach 

A. 2.1  Experimental  Configuration 

An  experimental  system  for  producing  a  nonequilibrium  recombining 
plasma  using  a  pulsed  photoionization  approach  is  illustrated  schemati¬ 
cally  in  Fig.  1.7.  The  radiation  source  is  a  pulsed  short  arc  Xe  lamp 
directed  along  the  axis  of  a  tube  containing  the  test  gas  (e.g.,  Cs 
vapor  and  Nj).  The  Cs  vapor  and  N2  are  contained  in  a  one-inch-diameter, 
twelve-inch-long  Inconel  tube  fitted  with  sapphire  windows.  The  Cs 
pressure  is  one  to  ten  torr  and  the  nitrogen  pressure  is  about  one  atmo¬ 
sphere.  An  electric  tube  oven  heats  the  CS-N2  mixture  in  the  temperature 


Fig.  1.7.  Experimental  system  based  on  pulsed  photoionization  approach. 


range  400  to  1200  K.  A  low  flow  rate  of  the  gas  mixture  through  the 


tube  contributes  to  maintaining  high  gas  purity.  The  sapphire  window 
has  75%  transmission  of  radiation  down  to  200  nm  wavelengths  so  it  is 
transparent  to  photons  capable  of  singly  ionizing  ground  state  Cs  atoms, 
i.e.,  X  ^  318  nm,  and  opaque  to  photons  capable  of  ionizing  Nj,  i.e., 
A  ^  200  nm.  The  Xe  flashlanp  provides  a  2 -ms  pulse  of  ionizing  photons 
which  is  coupled  to  the  Cs  test  cell  through  an  ellipsoidal  reflector. 


nonequilibrium,  isothermal  Cs  plasm  is  produced  by  the  radiation 
pulse.  The  free  electrons  are  rapidly  thermlized  through  collisons 
with  the  N2.  As  the  electrons  begin  to  recombine  with  Cs'*’  in  the  pres¬ 
ence  of  the  nitrogen,  the  electron  number  density  would  be  determined  as 
a  function  of  time.  Possible  methods  for  obtaining  the  electron  density 
include  measurements  of  the  spectral  emission  from  free -bound  and  high- 
lying  bound-bound  transitions  in  Cs,  as  well  as  measurements  of  the 
electrical  conductivity  using  a  four-pin  conductivity  probe.  The  time 
history  of  the  electron  number  density  is  directly  related  to  the  rate 
coefficient  for  recombination  on  N2. 


Experimental  conditions  would  be  chosen  so  that  the  recombination  on 
N2  proceeds  at  a  much  greater  rate  than  electron  loss  by  recombination 
on  electrons,  or  by  diffusion  of  electrons  to  the  test  cell  walls. 
Under  such  conditions  dn^/dt  =  -o^  nj^  ng,  so  that 


eo 


Thus  the  slope  of  l/n^  plotted  vs  time  gives  (T)  multiplied  by  the 
known  concentration  of  N2.  The  temperature  dependence  of  Oj^  (T)  is 
en^irically  determined  by  performing  the  experiment  over  the  range  of 
temperatures  that  can  be  achieved  with  the  apparatus  previously  de¬ 
scribed. 


The  electron  number  densities  corresoonding  to  the  parameter  value 
r  =  10  are  given  in  Table  1.1.  Also  shown  in  the  table  are  character¬ 


istic  times  for  recombination  on  N_ 


~  <®N2"n  electron 

^diff  "■  ''  ^‘^a  “  ionizing 

light  pulse  to  decay  from  a  peak  to  a  minimxim  intensity,  tp]^g,  and 
for  thermalization  of  free  electrons  “  ^'’e  N  •10*(2mg/r[y,  Here 


diffusion  to  the  walls,  t^ 


r2/d. 


sec 


Table  1 . 1 

Characterization  of  Pulsed  Experiment 


Temp,  K 

ng,  cm"^ 

t^  ,  sec 

”2 

400 

axio"*® 

7x10“^ 

600 

7x10''° 

8x10"^ 

800 

2x10^^ 

7x10"^ 

^diff' 


R  =  1  cm  is  the  radius  of  the  test  cell,  D  is  the  ^unbipolar  diffusion 


coefficient,  K  is  a  constant  that  depends  on  atomic  properties  of  Cs 
and  N2  [1*7],  is  the  average  e,N2  collision  frequency  [1-8],  zmd 


10*(2m-/m„  )  is  the  average  fraction  of  electron  kinetic  energy  lost  in 


all  collisions  with  N2  [1.9].  A  typical  value  for  the  flashlamp  pulse 
forming  network  is  used  for  tpj^g. 

Table  1.1  shows  that  the  characteristic  times  are  ordered  as  << 

tpj^g  <  tjj^  «  t^^££,  and  therefore  that  the  experimental  conditions 


result  in  an  isothermal  plasma  dominated  by  recombination  on  N2.  The 
degree  of  excitation  of  ^2  by  encounters  with  photons  or  energetic  elec¬ 
trons  is  negligible.  The  N2  molecules  do  not  significantly  absorb 
photons  with  wavelengths  greater  than  200  nm  [1.10].  The  energy  corre¬ 
sponding  to  an  electron  liberated  from  the  Cs  ground  state  by  a  200  nm 
photon  is  2.1  eV.  To  excite  a  ground  state  N2  to  its  first  vibration- 
ally  excited  state  requires  approximately  0.22  eV  [1.11],  or  each  elec¬ 
tron  can  excite  at  most  10  N2  molecules.  An  estimate  of  the  maximum 


fraction  of  N2  molecules  which  could  be  excited  in  this  manner  is  there¬ 


fore  10*n_/n„  <  10' 


A. 2. 2  Reaction  Kinetics  Model 


Zn  addition  to  processes  (1.1)  ^tnd  (1.2),  a  loss  of  electrons  under 
conditions  of  interest  could  also  occur  through  a  sequence  of  reactions 
involving  fast  formation  of  the  ion  cluster  Cs''’»N2,  fast  two-body 
switching  to  form  Cs^/  followed  by  dissociative  recombination  of  the 
cesium  dimer  ion.  These  processes  can  be  written  as  follows: 


Cs^  +  N2  +  ^2 


Cs'^»N2  +  ^2 


Cs^*N2  +  Cs 


Cs2  +  ^2 


Cs  +  Cs 


{1.3a) 

(1.4a) 

(1.5) 


Reactions  analogous  to  (1.3a)  and  (1.4a)  can  occur  when  water  vapor 


impurities  are  present  and  H2O  replaces  N2 

These  reactions  are  written: 

a  clustering  partner. 

Cs'*’  +  H2O  +  N2  — ♦  Cs''’*H20  +  N2 

(1.3b) 

Cs'*'.H20  +  Cs  — »  CS2  +  H2O 

(1.4b) 

Other  possible  mechanisms  of  electron  loss  are  electron  dissociative 
recombination  with  Cs'^.N2/  electron  attachme.it,  and  ion-ion  recombina¬ 
tion  as  represented  respectively  by  the  reactions 

e  +  Cs'^*N2  — ^  Cs  +  N2 

(1.6) 

e  +  Cs  +  Nj  — ♦  cs“  +  Nj 

(1.7) 

Cs“  +  Cs"*"  — ♦  2Cs 

(1.8) 

In  what  follows,  we  present  a  detailed  examination  of  the  overall 
recombination  process  that  includes  both  the  ion-molecule  reaction 
channel  (1.3),  (1.4),  and  (1.5)  as  well  as  the  three-body  processes  (1.1 
and  1.2).  Using  estimates  of  the  recombination  rate  coefficients  and 
equilibrium  constants  (to  calculate  the  baclcward  rates)  constructed  from 
available  information  in  the  literature  (see  Table  1.2  for  references), 
we  show  that  experimental  conditions  at  higher  temperatures  can  be 
chosen  to  maintain  process  (1.2)  as  the  dominant  electron  loss  mecha¬ 
nism,  despite  the  fast  reactions  (1.3)  and  (1.4).  Later  we  show  that 
inclusion  of  reactions  (1.6)-(1.8)  does  not  affect  the  dominance  of 
reaction  (1.2)  in  electron  loss  frequency.  In  addition  we  also  show 
that  by  using  a  liquid  nitrogen  cold  trap  and  because  of  the  gettering 
action  of  Cs,  as  well  as  from  considerations  of  the  reaction  kinetics, 
the  deleterious  effects  of  a  possible  water  impurity  are  not  signifi¬ 
cant. 

a.  Partial  Equilibrium 

Although  reactions  (1.3)  and  (1.4)  proceed  in  a  forward  direction 
rapidly  relative  to  the  three  other  reactions  considered,  the  backward 
rates  for  these  reactions  are  also  large.  Reactions  (1.3)  and  (1.4)  are 
able  to  quickly  reach  a  partial  equilibrium  condition,  where  the  forward 


and  backward  rate  are  equal  in  magnitude.  This  partial  equilibrium  is 
maintained  while  the  recombination  reactions  slowly  change  the  ion  and 
electron  concentrations. 

Table  1.2  lists  the  preceding  ten  chemical  reactions.  Also  listed 
are  the  values  of  the  corresponding  forward  and  backward  rate  constants 
and  the  ratios  of  these  two  constants  (i.e.,  the  equilibrium  constants), 
used  in  all  subsequent  calculations.  The  constants  were  obtained  using 
the  references  shown.  The  photoionization  of  Cs  vapor  by  ultraviolet 
radiation  from  a  flashlan^i  can  be  represented  by  the  reaction 

Cs  +  hv - ^  Cs"*"  +  e  . 

The  photon  flux  is  modeled  as  a  rectangular  function  in  time,  and  the  Cs 
vapor  cam  be  shown  to  be  optically  thin.  Under  these  conditions  the 
rate  of  increase  of  the  electron  number  density  can  be  represented  by 
d[e]/dt  =  [Cs]ky  where  k^  is  the  integral  over  wavelength  of  the  product 
of  the  photon  flux  and  the  total  photoionization  cross  section.  The 
value  for  k^  used  in  these  calculations  is  1.0  sec”^.  This  value  corre¬ 
sponds  to  a  photon  flux  (wavelengths  <  318  nm)  of  about  10^®  photons/ 
cm  •  sec,  which  is  an  estimate  of  the  flux  attainable  with  our  pulsed 
short  arc  flashlamp  system. 

Shown  in  Figs.  1.8a,  1.9a,  and  l.lOa  are  the  time  histories  of  the 
number  densities  of  the  species  of  interest  for  three  different  test 
cell  temperatures.  The  Cs'*^*H20  species  occurs  in  concentrations  which 
are  off  the  lower  end  of  the  plotted  scale.  This  species  is  treated 
explicitly  in  the  section  on  impurities.  (Note  that  reactions  (1.6)- 
(1.8)  are  not  modeled  in  the  results  presented  in  these  figures.)  It  is 
assumed  that  the  flashlamp  is  energized  at  t  =  0  and  it  is  extinguished 
at  t  =  1  ms;  the  dynamic  relaxation  of  the  nonequilibrium  plasma  occurs 
for  times  t  >  1  ms.  These  plots  were  made  by  numerically  solving  the 
set  of  differential  equations  which  describe  the  kinetics  of  the  species 
involved  in  reactions  (1.1)  through  (1.5).  For  t  <  0,  the  plasma  is 
in  a  state  of  local  thermodynamic  equilibrium.  The  equilibrium  spe¬ 
cies  concentrations  shown  in  the  figures  are  determined  by  the  Cs 
vapor  pressure  (controlled  by  the  reservoir  temperature),  by  the  initial 


I  108 


?  107 


•0.02  0  0.02  0.04  0.06  0.08 

time,  sec 

Pig.  1.8a.  Species  concentrations  from  numerical  solution  of  kinetics 
equations  for  reactions  (1.1)-(1.5).  T  =  1000  K;  [Cs]  = 

4.5x10^4  ci^-3.  _  .jqI^  cm“^;  [H2O]  =  10^^  cm~^.  Two 

distinct  regimes  are  shown;  t  <  0,  complete  equilibrium; 
t  >  1  ms,  plasma  relaxation.  Photoionization  occurs  between 
these  regimes. 
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Fig.  1.8b.  Electron  loss  frequencies  from  numerical  solution  of  kinet¬ 
ics  equations  for  reactions  (1.1)-(1.5).  T  =  1000  K;  [Cs]  = 
4.5x10^^  cm"8;  [Nj]  =  10^^;  tH20]  =  10^^  cm“^.  1/T^(e)  = 

(net  electron  loss  rate)/[e] .  For  ambipolar  diffusion, 

1/T^^^£(e)  ^  10  sec"^. 
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N2  concentration,  and  by  the  temperature  of  the  reaction  cheunber.  Typi¬ 
cal  values  for  these  quantities  were  taken  to  be  [Cs]  =  4.5x10^^  cm“^, 
[N2]  =  10''^  cm"^,  [H2O]  =  10"'^  cm“^,  and  600  K  <  T  <  1500  K. 

For  a  general  reaction  i  where 

A  +  B  - C  +  D 

we  can  obtain  or der-of -magnitude  estimates  of  the  characteristic  times 
for  significant  decrease  or  increase  of  concentration  of  species  respec¬ 
tively  by,  for  example,  T+(A)  ~  (k^LB])"'*,  or  t^CC)  ~  [C]  (k^ [A]  [B]  ) . 
In  what  follows,  it  is  shown  how  relations  such  as  these  may  be  employed 
to  illustrate  the  essential  physics  underlying  the  numerical  solutions 
presented  in  Figs.  1.8a,  1.9a,  and  1.10a.  Impurity  reactions  (1.3b)  and 
(1.4b)  have  been  included  in  the  numerical  solutions  shown  in  Figs. 
1.8-1.10,  but  an  explanation  of  the  effect  of  including  these  reactions 
is  deferred  until  the  section  on  impurities. 

Using  the  equilibrium  concentrations  for  t  <  0  shown  in  Fig.  1.9a, 
characteristic  times  for  significant  changes  in  species  concentrations 
are  calculated  in  Table  1.3,  for  a  temperature  of  1200  K.  Considering  the 

Table  1.3 

Characteristic  Times 

t  =  0,  T  =  1200  K,  Species  Concentrations  from  Fig.  1.9a 


duration  =  10“^  sec 


= 

(  [Cs^]  •  [N2]  •k2)“^  =  ( lo^oo’^oo"^®)"''  =  1  sec 

T5(e) 

= 

([CS2]*k5)“’  =  (3x10®. SxlO"^)"''  =  10^  sec 

Ti(e) 

= 

(  [Cs'*’]  .  [e]  •k^)"''  =  dO^.IO^.IO"^^)”''  =  10^  sec 

T^a(Cs+) 

= 

([N2]*k3a)“^  =  ( 10^^.6x10"''0)"^  =  10“^  sec 

Trees'*’) 

= 

(  [e]  •  [N2]  •k2)"^  =  (10^. 10^^. 10"^®)“'’  =  1  sec 

tL<Cs2) 

= 

[CS2]/(  [Cs''‘.N2]  ‘[Cs]  .k4a) 

=  3x10®/(8x10®.5x10’'*.10"®)  =  10"^  sec 

T^(CS^) 

= 

([ej.kg)"^  =  (10^.5x10“®)"''  =  1  sec 

KT  *r.  *n  JWJtnJWKruniwi^. 


important  reactions  affecting  the  electron  concentration,  the  ordering 
Tpuise  "  10”^  sec  <<  t'|[  <  Tg  <  <iemonstrates  that  the  duration  of  the 
photoionization  reaction  (one  millisecond)  is  sufficiently  short  to 
throw  the  plasma  out  of  equilibrium-  The  recombination  reactions  (l-l), 
(1.2),  and  (1-5)  will  not  affect  the  change  in  ion  and  electron  concen¬ 
trations  during  and  immediately  after  the  flash  because  these  reactions 
deplete  charged  particle  concentrations  in  a  characteristic  time  which 
is  much  longer  than  the  pulse  duration. 

The  orderings  T^g^(Cs‘*’)  «  T2(Cs'^)  and  «  15(032)  [note  that 
these  x's  are  calculated  at  equilibrium  where,  for  any  species  s,  T^{s) 
==  1^(3)]  indicate  that  reactions  (1.3a)  and  (1.4a)  will  proceed  to 
establish  partial  equilibrium  concentrations  of  Cs''’»N2  and  Cs^  while  the 
recombination  reactions  change  the  ion  and  electron  concentrations  rela¬ 
tively  slowly. 

This  partial  equilibrium  is  embodied  in  the  relations 


[Cs 

[Cs"^]  [N^l 


(1.9) 


[Cs]  [Cs"^] 


(1.10) 


Equilibrium  relation  (1.10)  corresponds  to  the  reaction  Cs^  — ^  Cs  +  Cs^ 
obtained  by  adding  reactions  (1.3a)  and  (1.4a).  Once  reactions  (1.3a) 
and  (1.4a)  are  shown  to  be  fast  enough  to  maintain  partial  equilibrium, 
the  equilibrium  relations  need  not  correspond  directly  to  the  reaction 
channel.  Relations  (1.9)  and  (1.10)  have  been  chosen  to  express  the 
equilibrium  condition  in  order  to  illustrate  that  the  equilibrium  con¬ 
centrations  of  CS2  and  Cs''’*N2  are  approximately  decoupled. 

Since  N2  and  Cs  are  present  in  much  larger  quantities  than  Cs^,  Cs^, 
or  Cs"'’»N2  (as  can  be  justified  by  the  results  of  the  calculation  being 
described),  the  following  major  species  relations  can  be  written; 


[N2]  =  constant  =  10'^  cm”^  , 


(1.11) 


I 

1^ 

1 


I 

I 

I 

1 
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[Cs]  =  constant  =  4.5x10^^  cm”^ 


(1.12) 


From  eqs.  (1.9)  and  (1.11)  and  Table  1.2,  the  ratio  [Cs‘''«N2]  /  [Cs'*']  = 
[N2]  •1'3a/^^-3a  shown  to  be  less  than  1%  at  all  temperatures 

greater  them  500  K.  Additionally,  it  has  been  shown  that  the  recombina¬ 
tion  reactions  are  too  slow  to  significantly  affect  charge  concentration 
during  the  flash.  Neglecting  the  Cs‘^»N2  concentration  and  the  effects 
of  recombination,  charge  conservation  during  the  flash  can  be  written 

[Cs"^]  +  [CS2I  =  t+1  »  (1.13) 

where  [+]  denotes  the  concentration  of  positive  (or  negative)  charge 
produced  by  the  photoionization  process.  The  charge  prod'^-'ed  by  the 
flash  is 

[+]  =  [Cs]*kg.At  =  4.5x10''^»1«10"^  =  4.5x10’’^  cm"^  .  (1.14) 

Equations  (1.10),  (1.12),  (1.13),  and  (1.14)  are  sufficient  to  determine 
equilibrium  concentrations  of  Cs^  and  Cs'*’;  the  value  for  [Cs'*’]  can  then 
be  used  with  eqs.  (1.9)  and  (1.11)  to  determine  the  Cs''’*N2  concentra¬ 
tion.  The  species  concentrations  obtained  in  this  way  are  shown  by  the 
curves  plotted  in  Fig.  1.11,  and  compare  well  with  the  results  of  the 
numerical  solution  at  t  =  1  ms  which  have  been  plotted  as  individual 
points  over  these  curves.  The  results  shown  in  Fig.  1.11  can  be  used  to 
predict  the  relative  rates  of  the  recombination  reactions  which  follow 
the  flash. 

b.  Characteristic  Times  for  Recombination 

Shown  in  Fig.  1.12  are  the  characteristic  reciprocal  times  1/t  for 
significant  electron  recombination  by  each  of  the  processes  (1.1), 
(1.2),  and  (1.5),  at  t  =  1  ms,  immediately  following  termination  of  the 
flash.  The  quantity  1/t  is  a  measure  of  the  frequency  of  electron  loss 
via  a  given  recombination  reaction.  The  continuous  curves  were  calcu¬ 
lated  using  the  formulas  shown  in  Table  1.3  evaluated  for  the  concentra¬ 
tions  given  in  Fig.  1.11.  At  temperatures  above  1000  K,  reaction  (1.2) 
is  clearly  the  dominant  loss  mechanism  immediately  after  the  flash. 


#,v 


The  numerical  solutions  shown  in  Figs.  1.8b,  1.9b,  and  1.10b  give 

detailed  information  about  electron  loss  mechanisms  for  extended  times 
after  the  flash.  These  plots  correspond  in  temperature  and  time  to  Figs. 
1.8a,  1.9a,  and  1.10a.  The  frequency  of  electron  loss  by  reaction  i, 

1  _  ^forward  -  backward  reaction  rate) 

_  - . - 


+,  ,  [e] 

T^(e) 

is  shown  for  recombination  reactions  (1.1),  (1.2),  and  (1.5)  during 
plasma  relaxation.  The  net  reaction  rate  (forward  -  backward  rate)  must 
be  considered  for  t  >  1  ms  when  the  reactions  can  approach  equilibrium. 
At  t  =  1  ms  the  net  rate  and  forward  rate  are  approximately  equal. 
Generally  the  electron  concentration  and  the  net  reaction  rate  change 
with  time,  so  that  this  frequency  is  usually  time  dependent. 

A  recombination  reaction  is  considered  to  be  the  dominant  electron 
loss  mechanism  when  its  electron  loss  frequency  exceeds  the  frequencies 
of  all  other  competing  loss  mechanisms.  If  measurement  of  electron  con¬ 
centration  decay  is  to  be  useful  in  determining  the  rate  of  reaction 
(1.2),  experimental  conditions  must  be  chosen  so  that  reaction  (1.2)  is 
the  dominant  electron  loss  mechanism.  For  all  temperatures  greater  than 
1000  K,  and  all  times  shown,  reaction  (1.2)  dominates  reactions  (1.1) 
and  (1.5)  in  electron  loss.  The  electron  loss  frequency  due  to  ambi- 
polar  diffusion,  VT^j^^^(e),  (at  about  1200  K  and  one  atmosphere  of  N2) 
is  approximately  constant  at  10  sec”^.  As  the  electron  concentration 
decays  with  time,  1/T2(e)  decreases  because  the  reactions  move  toward  an 
equilibrium  balance  where  the  forward  and  backward  rates  are  equal. 
Electron  concentration  measurements  are  only  useful  while  the  resulting 
reduced  electron  loss  frequency  remains  larger  than  VT^iff^®^'  Hence 
for  the  conditions  of  Figs.  1.8-1.10  electron  concentration  measurements 
should  be  made  in  approximately  the  first  10  ms  of  recombination.  Elec¬ 
tron  loss  frequencies  at  1  ms,  which  are  taken  from  numerical  solutions 
similar  to  those  shown  in  Figs.  1.8b,  1.9b,  and  1.10b,  are  plotted  as 
discrete  points  with  the  analytical  solution  shown  in  Fig.  1.12.  Both 
solutions  agree  that  at  temperatures  above  1000  K,  during  the  first 
10  ms  after  photoionization,  measurement  of  electron  concentration  decay 
will  provide  direct  information  about  reaction  (1.2). 
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Reaction  1.6  competes  directly  with  reaction  1.2  for  electrons,  but 
the  electron  loss  frequency  can  be  expected  to  be  low  due  to  its  de¬ 
pendence  on  the  low  concentration  of  Cs^»N2.  Reactions  1.7  and  1.8 
together  constitute  a  channel  for  electron  loss  via  recombination  with 
Cs'*’.  These  three  reactions  have  been  included  in  an  extended  computer 
model  of  the  kinetics.  The  results  of  the  numerical  calculation  of  con¬ 
centrations  and  electron  loss  frequencies  based  on  reactions  (1.1)-(1.8) 
are  shown  in  Figs.  1.13a  and  1.13b  (and  can  be  compared  with  Fig.  1.9). 
Figure  1.13b  indicates  that  reactions  (1.6)  and  (1.7)  are  significant 
secondary  electron  loss  mechanisms,  but  the  primary  electron  loss  mecha¬ 
nism  remains  reaction  (1.2).  Measurement  of  the  rate  of  electron  con¬ 
centration  decay  will  be  a  valid  indication  of  the  rate  of  reaction 
(1.2). 

c .  Impurities 

The  species  Cs,  Cs^,  e,  and  N2  are  necessary  for  the  experiment; 
other  chemical  species  which  may  be  formed  from  these  cannot  be  avoided. 
In  addition,  the  presence  in  the  test  cell  of  impurities  due  to  outgas- 
sing  or  impure  reactant  gases  must  be  considered.  Impurity  concentra¬ 
tions  can  easily  be  kept  small  enough  to  ensure  that  N2  is  the  dominant 
neutral  partner  in  all  reactions,  but  small  amounts  of  radical  species 
could  substantially  affect  ion-molecule  chemistry. 

Water  vapor  is,  of  all  impurities,  most  likely  to  be  both  suffi¬ 
ciently  reactive  and  plentiful  to  have  a  significant  effect.  The  amount 
of  H2O  present  in  the  N2  charge  will  be  minimised  by  passing  the  gas 
through  a  liquid  nitrogen  trap.  The  partial  pressure  of  the  H2O  will 
roughly  correspond  to  the  vapor  pressure  at  100  K,  ~  10”^®  torr 

[1.30].  At  300  K  this  pressure  corresponds  to  [H2O]  =  5x10®  cm”^.  For 
a  baked  and  pumped  metal  system  a  typical  outgassing  rate  is  10~^® 
torr* liter/sec .cm  [1.31].  Assuming  H2O  is  a  major  constituent  of  the 
outgas  and  using  design  dimensions  of  the  test  cell  with  a  period  of 
about  30  hours  between  the  end  of  pumping  and  the  end  of  the  experiment, 

Pr  o  “  10”^®(torr*liter/sec*cm^) *60  cm^*10®  sec/1  liter  =  6x10“^  torr. 

2  1  "i  -3 

At  300  K  this  pressure  corresponds  to  [H2O]  =  2xlo  cm  .  For  one 

atmosphere  of  N2  this  H2O  concentration  corresponds  to  an  impurity  level 
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Fig.  1.13a.  Species  concentrations  from  numerical  solution  of  kinetics 
equations  for  reactions  (1.1)-(1.3).  T=  1200  K.  The  fol¬ 
lowing  initial  species  concentrations  remain  roughly  con- 
"  4e5xl0  ^  cm  (^23  =  10^^  cm"^;  [H2O]  ^  10^^ 
cm  .  distinct  regimes  are  shown;  t  <  0,  complete 

equilibrium;  t  >  1  ms,  plasma  relaxation.  Photoionization 
occurs  between  these  regimes. 
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Fig.  1.13b.  Electron  loss  frequencies  from  numerical  solution  of  kinet 
ics  equations  for  reactions  (1.1) -(1.8).  T  =  1200  K.  Th 
following  initial  species  concentrations  remain  roughly  con 
stant:  [Cs]  =  4.5xio  ^  cm"^;  [N,]  =  1o’®;  [H-O]  =  10^^  cm”^ 
1/t,  =  (d[e]/dt),/[e] . 
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of  one  part  per  million.  However,  the  water  present  will  be  very  effec¬ 


tively  scavenged  by  reducing  reactions  with  Cs  vapor.  A  reduction  in 


H2O  concentration  by  a  factor  of  300  has  been  achieved  by  introducing  Cs 
vapor  into  a  10”^  torr  vacuum  [1.32].  In  light  of  scavenging  [H2O]  = 
10^^  cm”^  appears  to  be  a  conservative  worst-case  estimate. 


Even  assuming  an  amount  of  water  vapor  impurity  equal  to  this  worst- 


case  estimate,  it  can  be  shown  that  ion-molecule  chemistry  of  water 


vapor  in  the  plasma  will  not  interfere  with  the  primary  experimental 
goal  of  measuring  the  rate  of  reaction  (1.2).  The  association  and 


switching  reactions  (1.3)  and  (1.4)  can  occur  with  N2  replaced  by  H2O; 

related  bond  energies  are  E^  4.  „  =  6.7  Iccal/mol  [1.12],  and  E  4.  = 

Cs^»N2  Cs^*H20 


13.7  kcal/mol  [1.13] 


In  relation  to  the  measurement  of  the  rate  of 


reaction  (1.2),  reactions  (1.3)  and  (1.4)  provide  a  channel  for  the  pro¬ 
duction  of  CS2  (which  subsequently  dissociatively  recombines)  and  a  loss 
mechanism  for  Cs'*’.  It  has  been  shown  for  the  relatively  fast  reactions 
involving  the  Cs''’*N2  cluster  that  a  partial  equilibrium  calculation  can 
be  used  to  determine  species  concentrations  before  recombination  begins. 
Partial  equilibrium  represents  a  worst  case  since  it  allows  the  greatest 
advancement  down  the  reaction  channel  which  produces  Cs^  prior  to  the 
onset  of  recombination.  Partial  equilibrium  may  therefore  be  conserva¬ 
tively  assumed  for  the  reactions  involving  the  Cs''’»H20  cluster  in  a 
calculation  which  is  analogous  to  the  Cs^»ll2  partial  equilibrium  compu¬ 
tation.  Agreement  between  analytical  and  numerical  solutions  for  Cs'*’*H20 
concentration  in  Fig.  1.11  supports  the  validity  of  the  partial  equilib¬ 
rium  assumption.  For  equilibrium  involving  Cs'*’»H20  the  governing 
equations  are  (1.10),  (1.12),  (1.13),  (1.14), 


[Cs  .H^O] 


(1.15) 


[Cs  ] [H^O] 


[H2O]  =  constant  =  10 


13  ,^-3 
crn 


(1.16) 


Equation  (1.13)  still  accurately  describes  charge  conservation  since 
from  eqs.  (1.15)  and  (1.16)  and  Table  1.2,  [Cs"*" •H2O] / [Cs^]  =  [H2O]  • 

k2]j/k_3jj  <  0.01%  at  all  temperatures  greater  than  500  K.  Both  the 
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Cs'*’*N2  and  Cs'''*H20  clusters  have  been  shown  to  occur  in  such  small  con¬ 
centrations  ( relative  to  Cs'*’)  that  the  Cs^  equilibrium  concentration  can 
be  determined  independently.  In  Pig.  1.11  the  values  for  the  Cs^  and 
Cs'*’  concentrations  have  been  calculated  first  and  then  used  to  determine 
the  Cs''’*N2  and  Cs^*H20  concentrations.  The  Cs^*H20  cluster  will  not 
provide  a  significant  loss  mechanism  for  Cs'*’  and  the  presence  of  water 
vapor  will  not  affect  the  equilibrium  concentration  of  Cs^  calculated 
from  eqs.  (1.10),  (1.12),  (1.13),  and  (1.14).  Since  the  Cs^  concentra¬ 
tion  is  unchanged,  previous  statements  about  the  negligible  effect  of 
dissociative  recombination  of  Cs^  still  hold  in  the  presence  of  water 
vapor  impurities . 

d.  Conclusions  from  Reaction  Kinetics  Model 

A  model  has  laeen  constructed  to  assess  the  role  of  several  ion- 
molecule  reactions  that  can,  in  principle,  be  present  under  conditions 
of  relatively  low  electron  number  density  and  high  ambient  pressure, 
where  experiments  will  be  performed  to  measure  the  electron  recombina¬ 
tion  rate.  The  results  of  this  study,  based  on  numerical  and  analytical 
solutions  of  a  full  set  of  )cinetic  equations  with  rate  coefficients 
estimated  from  information  in  the  literature,  indicate  that  for  tempera¬ 
tures  in  the  range  of  1000  K  to  1500  K,  electron  decay  rate  measurements 
during  the  first  10  ms  of  plasma  relaxation  can  be  identified  with  the 
three-body  process  e  +  Cs^  +  N2  — >  Cs  +  N2,  and  can  therefore  be  used  to 
test  theoretical  predictions  for  this  process.  0\ir  analysis  of  possible 
difficulties  associated  with  the  presence  of  a  water  vapor  impurity  indi¬ 
cates  that  in  this  temperature  range  the  impurity  effect  is  negligible. 

It  should  be  noted  that  considerable  uncertainty  exists  for  the 
values  of  many  of  the  rate  constants  used  in  this  study.  However,  the 
study's  conclusions  are  based  on  robust  argroments  about  partial  equilib¬ 
ria  and  the  ordering  of  time  constants  which  generally  demonstrate  that 
two  reaction  rates  or  two  characteristic  times  differ  by  many  orders  of 
magnitude.  The  strength  of  the  inequalities  involved  malte  these  conclu¬ 
sions  insensitive  to  reasonable  uncertainties  in  rate  constants. 


In  essence,  what  is  needed  to  determine  (T)  is  a  neasurement  of 
the  recombination  rate  in  a  plasma  which  is  out  of  ionizational  equilib¬ 
rium  and  which  contains  only  a  single  molecular  species  that  can  act  as 
a  third  body  for  recombination.  The  pulsed  photoionization  experiment 
satisfies  these  requirements  while  avoiding  many  of  the  technical  diffi¬ 
culties  inherent  in  the  steady -state,  rapid  expansion  experiment.  Since 
the  plasma  is  stationary  in  the  pulsed  experiment,  diffusion  (not  con¬ 
vection)  determines  the  residence  time  of  a  plasma  differential  volume 
element.  The  characteristic  time  for  the  ambipolar  diffusion  of  a  Cs 
plasma  in  a  one -inch-diameter  tube  containing  N2  at  a  pressure  of  one 
atmosphere  is  about  500  ms,  at  600  K,  or  4000  times  longer  than  the 
convection-limited  residence  time  in  a  10 -cm  long  sonic  nozzle  at 
2000  K.  If  the  previously  derived  formula  for  r/t  is  used  with  tj.gg 
replaced  by  one  obtains  r/t  =  700  at  T  =  400  K  and  r/x  =  250  at 

T  =  800  K.  For  the  photoionization  experiment,  it  is  necessary  to  have 
T  [=(01^  ”e^~^/^diff^  ^  contrast  to  the  condition  t  =  1  for  the 

rapid  expansion  experiment).  It  is  possible  to  have  this  condition  on  t 
satisfied  and  at  the  same  time  have  r  ^  10  in  the  pulsed  experiment. 

The  photoionization  experiment  creates  an  isothermal  plasma  because 
the  liinetic  energy  of  free  electrons  after  photoionization  is  quickly 
thermalized.  Recombination  occurs  at  a  constant  temperature,  and  the 
temperature  dependence  of  (T)  is  measured  directly  by  simply  repeat¬ 
ing  the  experiment  at  different  tenperatures .  The  large  r/x  ratio 
obtainable  in  the  pulsed  experiment  allows  both  r  and  x  to  be  chosen 
independently  to  eliminate  ambiguity  alxjut  the  role  of  the  conpeting 
recombination  on  electrons  and  about  the  role  of  diffusion.  In  these 
respects  the  pulsed  photoionization  experiment  appears  to  be  more  favor¬ 
able  than  the  steady -state  rapid  expansion  experiment. 


irimental  Work 


Experimental  work  on  the  recombination  project  has  been  oriented 
toward  solving  the  technical  problems  associated  with  a  reaction  cell, 
handling  of  reactants,  photoionization  light  sources,  and  diagnostics. 


Our  experimental  efforts  in  each  of  these  areas  are  summarized  below 
along  with  pertinent  information  from  the  literature.  Based  on  this 
work  an  experimental  system  is  described  in  which  rates  for  electron-ion 
recombination  on  molecular  neutrals  can  be  measured. 

B. 1  Reaction  Cell 

The  photoionization  experiment  imposes  demanding  specifications  on 
the  environment  in  which  the  recombination  reaction  should  be  studied. 
Nitrogen  number  density  on  the  order  of  10  cm  ,  with  cesium  vapor 
concentrations  of  about  one  part  per  million  and  inpurity  concentrations 
of  much  less  than  one  ppm  must  be  maintained.  Optical  access  for  the 
photoionization  of  the  cesium,  and  access  for  species  concentration 
measurements  must  be  provided.  Pressure  differences  between  the  reac¬ 
tion  and  ambient  environments  of  about  one  atmosphere  must  be  supported. 
An  isothermal  region  must  be  heated  up  to  temperatures  as  high  as 
1200*C.  This  section  summarizes  our  investigation  into  the  technologies 
of  materials,  thermal  design,  seals,  and  windows  that  can  provide  the 
specified  reaction  environment. 

B.1.1  Materials 

Table  1.4  presents  a  selection  of  materials  which  are  known  to  be 
compatible  with  cesium.  Next  to  the  listed  naterials  are  the  approxi¬ 
mate  maximum  working  temperatures  for  exposure  to  cesium,  the  informa¬ 
tion  reference,  and  pertinent  comments.  In  general  glasses  are 

available  for  work  with  cesium  at  low  temperatures,  but  metals  or 
ceramics  must  be  used  at  higher  temperatures.  As  temperature  is 
increased,  cesium  reduces  the  Si02  in  glass.  Even  at  low  temperatures 
the  oxides  in  glass -to -metal  seals  are  rapidly  attacked  by  cesium  vapor. 
Some  metals  (commonly  used  as  brazing  alloys',  which  either  react  rapidly 
or  have  high  solubility  in  alkali  metals,  arc  -.ntimony,  bismuth,  cadmium, 
calcium,  gold,  lead,  silver,  selenium,  sulfur,  platinum,  and  tin  [1.33]. 
Ceramics  can  be  useful  as  dielectric  materials  for  the  reaction  cell 
walls,  and,  in  crystalline  form,  as  optical  window  materials.  Single¬ 
crystal  sapphire  has  excellent  optical  transmission  qualities,  high 
hardness,  and  it  is  commercially  available.  We  have  used  it  for  cesium 
absorption  measurements  and  photo ionization  experiments.  Lucalox  [1.34] 


Table  1 . 4 

Cesium  Compatible  Materials 


Maximum 

Material 

Temp . ,  ®C 

Ref. 

Comment 

Glasses 

Aluminosilicates  500 

Corning  1720 
Corning  1723 

Borosilicates  400 

Corning  7740 
Schott  8329 
Other  glasses 


Gehlinite  600 

Calcium  Aluminate  700 

Schott  8436  540 

Refractory  Metals  >  900 

Molybdenum  550 

Moly-50  Rhenium  1870 

Tantalum  2000 

Niobium  1400 

Ferrous  Metals 

304  S.S.  870 

Iron,  Steel,  Kovar  >  700 

Nickel  and  Alloys  >  700 

Haynes-25,  TD  Nickel  1100 

Copper  and  Alloys  5-600 

Ceramics 

Lucalox  1200 

Polycrystalline  AI2O2 
Single-Crystal  AI2O3  700 


[1.35]  Lower  temp,  for  Na  and  K. 

[1.36]  Temp,  dependent  transmis¬ 
sion  below  400  nm. 

[1.35]  (Pyrex).  Slowly  discolor 
at  higher  temp. 

[1.35] 

Poor  availability. 

Limited  shapes. 

400  nm  transmission  cutoff. 

[1.35] 

[1.37] 

[1.37] 

[1.38] 

[1.37] 

[1.37]  Superior  to  316  S.S. 

[1.35] 

[1.39]  Inconel,  Nichrome,  Monel. 

[1.39] 

[1.35] 

[1.35] 

[1.40]  High-pressure  sodium  lamps. 

[1.41]  Temp,  could  be  higher. 

[1.35] 


is  a  translucent,  polycrystalline  alumina  product  which  is  far  less 
expensive  than  sapphire.  Imaging  through  Lucalox  is  severely  impaired 
by  scattering,  but  Lucalox  can  be  used  as  a  window  material  for  photo¬ 
ionization  where  diffuse  radiation  is  acceptable.  Over  90%  of  the  light 
from  a  source  inside  a  standard  Lucalox  tube  is  claimed  to  be  transmitted 
through  the  tube  wall.  Other  ceramic  crystals  which  exhibit  good  alkali 
metal  compatibility  include  CaF2,  BaF2,  BeO,  ^2^3'  [’.35]. 
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Cesium  compatibility  may  not  be  the  limiting  factor  in  materials 
selection.  The  reaction  rate  of  some  naterials  with  other  reactant  or 
ambient  gases  can  be  excessive  at  temperatures  below  their  cesium  com¬ 
patibility  limit.  For  exan^le,  niobium  is  oxidized  rapidly  above  400® C, 
reacts  with  water  at  500®  C,  and  is  embrittled  by  Nj  and  by  H2  above 
400®C  [1.42,  1.43].  Additional  factors  which  affect  material  selection 
are  mechanical  strength  bo  support  typical  loads  at  operating  condi¬ 
tions,  vapor  pressure  at  elevated  temperatures,  and  ease  of  fabrication. 

Our  first  reaction  cell  design,  as  shown  in  Fig.  1.7,  provided  an 
example  of  conplicating  factors  in  material  selection.  Inconel  601  was 
initially  chosen  for  the  cell  wall  because  of  its  cesium  conpatibility 
and  its  excellent  resistance  to  oxidation  in  air  at  tenperatures  up  to 
1200®C.  However,  Inconel  proved  to  be  very  difficult  to  weld  to  commer¬ 
cial  stainless  steel  vacuum  flanges.  This  problem  was  circumvented  by 
brazing  the  flanges  to  the  Inconel.  High-temperature  vacuum  bakeout  of 
the  Inconel  resulted  in  dark,  sooty  deposits  on  the  windows  and  cooler 
parts  of  the  reaction  cell.  Electron  dispersive  spectroscopic  analysis 
of  this  deposit  indicated  it  was  composed  of  chromium  and  manganese  or 
their  oxides.  These  elements  are  present  in  Inconel  in  respective 
proportions  of  20%  and  1%.  These  deposits  could  be  avoided  for  reaction 
cell  tenperatures  above  450®  C,  only  when  at  least  an  atmosphere  of 
nitrogen  is  present.  For  later  experiments  we  used  a  stainless  steel 
reaction  cell.  Possible  alternatives  include  nickel  plating  the  inside 
of  the  reaction  cell  or  changing  to  a  ceramic  material. 

B.1.2  Seals 

The  fabrication  of  our  alkali  metal  vapor  system  required  seals 
between  conponents  such  as  the  reaction  cell,  windows,  the  metal  reser¬ 
voir,  and  electrical  feedthroughs.  Often  the  conponents  were  of  differ¬ 
ent  materials,  and  demountable  seals  were  desired.  The  seals  had  to  be 
gas  tight  and  resist  cesium  attack  while  operating  at  elevated  tenpera¬ 
tures.  Sealing  has  been  a  challenging  technical  obstacle  in  all  our 
system  designs. 

Metal  to  metal  joining  can  usually  be  satisfactorily  accomplished  by 
tungsten  inert  gas  welding  and  by  copper  or  nickel  brazing.  These 
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joints  will  remain  leak  tight  until  the  materials  begin  to  react  with 
the  gases  present  or  until  they  melt.  A  convenient  demountable  joint 
can  be  made  with  commercially  available  stainless  steel  vacuum  flanges 
(e.g.,  from  Varian  Corp.).  These  flanges  use  opposing  knife-edges  which 
deform  a  copper  metal  gasket.  At  400“C  these  flanges  performed  well  in 
our  experiments  with  cesium.  Bolt  thread  galling  was  prevented  with 
nolybdenum  disulfide  lubricant  and  silver-plated  nuts.  Upon  disassem¬ 
bly,  it  was  sometimes  necessary  to  employ  a  wedge  to  separate  the 
flanges  and  then  the  copper  gasket  had  to  be  pryed  off  a  flange.  These 
vacuum  flanges  have  been  reported  to  seal  well  up  to  650*C  [1.35].  An 
excellent  demountable  seal  for  snail -diameter  metal  tubing  is  the  VCR 
coupling.  (VCR  is  a  fitting  manufactured  by  the  Cajon  Coii5>any  [1.44].) 
With  nickel,  stainless  steel,  or  copper  gaskets  VCR  couplings  are  rated 
to  5370C. 

Alumina  can  be  joined  to  metal  by  metalization  and  brazing,  by  an 
oxide  frit  seal,  or  by  demountable  o-rings  and  gaskets.  A  brazed  joint 
is  nade  by  applying  metalizing  paint  to  the  alumina  and  firing  it  at 
high  temperature.  Alkali-resistant  paints  should  not  contain  silica; 
typically  they  are  tungsten  based  with  2%  yittria  oxide,  and  they  fire 
at  1550*  to  1825*C.  The  paint  is  then  plated  with  nickel  and  joined  to 
the  metal  part  with  a  copper  or  nickel  braze  material.  These  joints  are 
reliable  to  about  700*C,  but  the  seal  must  be  designed  to  minimize 
stress  due  to  differences  in  thermal  expansion  [1.45].  Sapphire  windows 
brazed  to  kovar  are  commercially  available  in  vacuum  flanges.  We  have 
used  two  such  windows  [1.46]  at  400®C  for  over  100  hours  of  exposure  to 
cesium  liquid  or  vapor.  One  of  the  windows  developed  a  small  leak  at 
the  brazed  joint  after  about  50  hours;  the  other  window  is  still  vacuum 
tight,  yet  both  windows  are  brazed  with  silver  (which  is  not  coii?)atible 
with  cesium) .  Varian  Corporation  is  laelieved  to  nake  a  sapphire  vacuum 
window  which  is  alkali  resistant. 

An  alkali -resistant  frit  seal  between  niobium  and  alumina  has  been 
developed  for  sodium  lan5)s.  The  frit  material  is  coti^osed  of  oxides  of 
calcium,  aluminum,  and  magnesium.  Niobium  is  used  Isecause  it  is  a  good 
thermal  expansion  match  for  alumina  and  will  not  melt  at  the  1460»C 
temperature  required  to  nake  the  seal  [1.35].  Frit  seals  can  run  for 
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thousands  of  hours  at  700*C  [1.40].  lypically  niobium  is  operated  in  a 
controlled  atmosphere  because  of  its  reactions  with  oxygen,  water, 
nitrogen ,  and  hydrogen . 


Metal  o-rings  or  gaskets  can  be  used  to  make  a  demountable  seal 
between  metal  and  sapphire.  Figxire  1.14  illustrates  an  o-ring  design  we 
have  tried.  A  sapphire  window,  two  inches  in  diameter,  and  0.060" 


thick,  was  compressed  between  two  metal  o-rings  which  were  located  in 


standard  vacuum  flanges.  The  o-rings  were  fashioned  from  copper- 


jacketed  Inconel  springs  which  were  designed  to  operate  at  temperatures 
up  to  450*C.  They  required  2000  Ibs/linear  inch  loading  to  seal  [1.47]. 
The  sapphire  easily  withstood  the  compressive  stress,  but  great  care  was 


needed  to  avoid  any  asymmetric  loading  which  could  create  bending 


moments  and  tensile  stress.  Windows  were  broken  in  assembly  and  during 


operation.  When  the  window  survived,  the  seal  usually  failed  after  one 
temperature  cycle  and  required  o-ring  replacement.  We  cannot  recommend 
this  metal  o-ring  seal  as  a  reliable  technique.  Sapphire  windows  have 
been  sealed  between  copper  gaskets  in  a  VCR  fitting  [1.35].  This 
arrangement  had  a'  leak  rate  which  increased  with  the  operating  tempera¬ 
ture  and  with  the  number  of  temperature  cycles. 
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Fig.  1.14.  High-temperature,  cesium  compatible,  vacuum  window  design. 
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Vacuum-tight,  heated,  demountable  seals  between  dissimilar  materials 
are  easily  achieved  with  elastomeric  o-rings  if  the  joints  can  be 
designed  to  operate  at  less  than  200*c  (Viton)  or  SOO^C  (Kalrez)  and  to 
be  free  from  alkali  metals.  Kalrez  costs  about  40  times  as  much  as 
Viton.  For  low-temperature  operation  in  the  presence  of  alkali  metals, 
the  best  o-ring  materials  have  been  found  to  be  Kel-F  3700  and  Buna-N 
tested  in  120*C  sodixun  vapor  for  nine  days,  and  Neoprene  and  Buna-N 
tested  in  180®C  sodium  vapor  for  seven  days.  Kel-F  and  Neoprene  are 
halogenated  and  they  could  potentially  react  explosively  with  alkali 
metals  [1.48] . 

B .  1 . 3  Temperat\ire 

If  cesium  vapor  contacts  regions  of  the  cell  which  are  colder  than 
the  saturation  (dew  point)  temperature  of  the  vapor,  cesium  will  con¬ 
dense  there.  The  saturation  tempera t\ire  for  our  experiments  is  typi¬ 
cally  in  the  range  200®-300®C.  To  prevent  the  condensation  and  loss  of 
cesium,  every  part  of  the  cell  must  be  heated  to  at  least  the  dew  point 
temperature,  as  in  a  hot-cell  design,  or  the  vapor  must  be  prevented 
from  reaching  the  cold  regions,  as  in  a  heat  pipe  oven  design.  Avail¬ 
able  sealing  technology  limits  the  maximum  sapphire  window  or  ceramic 
feedthrough  temperature  to  less  than  700®C.  A  reaction  environment 
maximum  temperature  of  1200*C  implies  that  temperature  gradients  will 
necessarily  exist  between  the.  windows  or  feedthroughs  and  the  region 
where  recombination  is  measured. 


A  hot  cell  was  constructed  in  the  course  of  the  present  work.  The 
heating  system  was  composed  of  three  separate  ovens.  Referring  to  Fig. 
1.15,  these  were  F4,  the  reservoir  over  (25®-400®C)  ,  F8,  the  reaction 
cell  oven  (400* -1000®C) ,  and  F9,  the  box  oven  (400®C)  The  respective 
purposes  of  the  ovens  were  to  control  the  cesium  vapor  pressure,  to  heat 
reactants  to  elevated  temperatures,  and  to  maintain  vacuum  windows  and 
seals  at  a  temperature  above  that  of  the  reservoir  but  below  the  operat¬ 
ing  limit  of  450®C.  Figure  1.7  shows  a  typical  high-temperature  profile 
through  the  reaction  cell.  The  experiment  was  designed  to  measure  reac¬ 
tion  rates  in  the  central  isothermal  region  of  the  cell. 
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VALVES 

1.  regulator  shut  off,  i  loke 

2. needle  0-400  cc/min,  Porter 
O.bakeable  bellows,  Nupro 

4. bakeable  bellows,  Nupro 

5. needle  0-400  cc/min.  Porter 
6.on-off  inline  o-ring 

7. bellows  o-ring 

8.  bakeable  bellows  1  1/2',Varian 

9.  bellows  o-ring 

10. gate  valve  o-ring 


METERS 

1  .flowmeter  50-300  cc,min.  Porter 

2.  diff.  press,  gage,  800  torr,  W&T 

3.  ion  gage 

4.  thermocouple  gage 


FIXTURES 
1  .Cs  trap 

2.  Cs  ampuole  holder 

3.  Cs  trap 

4.  reservoir  oven 

5.  diffusion  pump 

6.  zeolite  trap 

7.  roughing  pump 

8.  hi  temp  oven 

9.  med.  temp  oven 


Fig.  1.15.  Heating,  vacuum,  and  gas  flow  systems  for  a  hot  cell  design. 

Heat  pipe  ovens  [1.49]  utilize  evaporation,  condensation,  and  a 
metal  wick  to  maintain  cold  cell  ends  in  which  a  nonreactive  buffer  gas 
is  confined,  and  a  heated  interior  in  which  a  homogeneous,  steady-state 
concentration  of  alkali  metal  vapor  is  confined.  In  typical  operation, 
the  metal  vapor  pressure  is  equal  to  and  controlled  by  the  pressure  of 
the  buffer  gas.  The  vapor  and  gas  do  not  mix  except  in  a  thin  inter¬ 
face.  Materials  such  as  o-rings  vriiich  only  operate  at  low  temperature 
without  alkali  contact  can  be  used  on  the  ends.  These  ovens  have  been 
operated  in  a  regime  where  the  metal  vapor  pressure  was  determined  by 
tes^rature  control  of  the  wick,  so  the  buffer  gas  pressure  could  be 
arbitrarily  selected  [1.50].  The  central  region  of  the  heat  pipe  has 
been  operated  without  a  wick  and  it  has  been  heated  to  create  a  region 
of  superheated  vapor  and  buffer  gas  of  arbitrary  number  density  [1.51, 
1.52].  Turbulence  and  alkali  mist  or  rain  has  been  observed  in  the 
condensation  region  of  heat  pipe  ovens  which  were  operated  at  high 
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buffer  gas  and  alkali  vapor  pressures.  The  turbulence  and  mist 
adversely  affect  optical  transmission  through  the  pipe,  but  conditions 


at  the  center  and  ends  are  not  affected  [1.35]. 

B.2  Gas  Handling 

The  design  of  a  system  to  supply  gases  to  the  reaction  environment 
is  Influenced  by  chemical  reactions  with  inpurlties  and  by  mass  trans¬ 
port.  These  reactions  can  have  two  effects:  First,  cesium  atom  reac¬ 
tions  with  in5>urities  can  reduce  the  concentration  of  cesium  vapor  which 
is  available  for  photoionization,  and  secondly,  reactions  between 
charged  species  and  inpurities  duuring  recombination  can  interfere  with 
the  intended  measurement  of  the  Cs'*'-e~-N2  recombination  reaction.  Mass 
transport  is  necessary  to  move  cesium  to  the  reaction  zone  as  it  vapor¬ 
izes  from  a  liquid  state.  Other  mass  transport  processes,  such  as 
diffusion  to  a  colder  region,  may  simultaneously  reduce  the  vapor  con¬ 
centration.  The  net  itmss  transport  must  pro(^ce  a  steady-state  vapor 
concentration  in  the  reaction  zone. 

The  cesium  vapor  design  concentration  is  nominally  1  ppm  in  10^^ 
cm“^  of  nitrogen,  but  the  design  figure  could  vary  between  0.1  ppm  and 
10  ppm.  Cesium  vapor  is  a  strong  reducing  agent,  and  it  readily  reacts 
with  oxygen  and  water  vapor.  The  concentrations  of  these  common  inpuri¬ 
ties  must  be  maintained  far  below  the  lowest  cesium  concentration  to 
prevent  substantial  loss  of  cesium  (in  this  case,  much  less  than 
0.1  ppm).  With  respect  to  the  second  effect  of  impurities,  kinetics 
modeling  of  the  role  of  water  vapor  as  a  participant  in  an  undesirable 
electron  loss  mechanism  indicated  that  water  concentrations  of  up  to 
1  ppm  would  be  acceptable.  In  the  case  of  water  vapor,  reactivity  with 
cesium  determined  the  upper  limit  of  acceptable  impurity  concentration. 
Potential  inpurities,  such  as  hydrocarbons  or  chromivim,  may  not  react 
directly  with  cesium  vapor,  but  they  miight  form  negative  ions,  clusters 
with  cesium  ions,  or  otherwise  react  with  charged  species  during  recom¬ 
bination.  We  have  not  modeled  the  effect  of  such  impurities  on  the 
recombination  kinetics,  but  have  assumed  that  a  reaction  environment 
which  is  pure  enough  to  preserve  low  cesium  concentrations  will  not 
introduce  undesirable  charged  species  reactions. 
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Impurities  can  enter  the  system  through  leaks,  from  outgassing  sur¬ 
faces,  and  with  entering  reactants*  Leak  rates  can  usually  be  kept 
acceptably  low  with  proper  seal  design.  If  leaks  are  unavoidable,  the 
reaction  cell  pressure  should  be  higher  than  ambient  pressure,  so  reac¬ 
tants  will  flow  out  and  impurities  will  not  flow  in.  Outgassing  can  be 
reduced  by  preparatory  cleaning.  Standard  vacuum-cleaning  procedure, 
including  degreasing,  acid  pickle,  deionized  water,  and  alcohol  wash 
should  be  followed  [1.53].  The  system  should  be  baked  out  at  operating 
temperature  or  at  least  200*C,  in  vacuum  better  than  10“^  torr  or  under 
flowing  gas  of  high  purity,  before  cesium  is  admitted.  A  flowing  system 
which  continually  flushes  the  reaction  cell  with  a  fresh  supply  of  reac¬ 
tants  can  reduce  the  concentration  of  impurities  due  to  leaks  and  out¬ 
gassing  by  increasing  the  flow  rate.  Because  cesium  vapor  is  constantly 
mixed  with  new  nitrogen  in  a  flowing  system,  the  nitrogen  purity  must  be 
very  high  to  avoid  depletion  of  the  cesium.  Nitrogen  with  less  than  5  ppm 
total  impurities  is  commercially  available,  and  purification  to  <  1  ppm 
impurities  is  commonly  done  in  industry,  but  it  can  be  costly.  A  con¬ 
venient  gettering  technique  for  impurities  which  react  with  alkali  metals 
is  to  bubble  the  nitrogen  through  hot  NaK  [1.39] .  Gettering  is  inherent 
in  a  static  system,  such  as  a  heat  pipe  oven,  which  is  charged  with  a 
fixed  amount  of  nitrogen  and  cesium.  The  nitrogen  is  purified  by  reac¬ 
tions  with  cesium  before  the  cesium  charge  is  depleted.  Lower  purity 
standards  for  the  nitrogen  charge  are  acceptable  in  a  static  system  than 
in  a  flowing  system.  Additionally,  a  heat  pipe  tends  to  move  impurities 
away  from  the  reaction  zone  with  constant  evaporation,  condensation,  emd 
wicking  action.  Impurities  in  the  cesium  are  generally  other  alkali 
metals,  and  the  most  common  purification  technique  is  vacuum  distillation 
[1.39].  To  prevent  contamination  during  installation  the  cesium  should 
be  handled  under  positive  pressure  of  cover  nitrogen  or  under  vacuum. 

Techniques  of  mass  transport  vary  with  system  design,  but  diffusion 
velocity  is  usually  an  important  parameter.  Molecular  diffusion  of 
relatively  low  concentrations  of  cesium  vapor  through  nitrogen  has  a 
binary  diffusion  coefficient  of  [1.53] 
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Here  denotes  the  cesium  mean  free  path,  is  the  cesium  average  j 

thermal  speed,  d  (=  4.5  A)  is  the  average  atomic  di^uneter  for  Cs  and  N2/  | 

p  is  the  pressure  of  Nj,  k  =  1.38x10“^^  J/K,  and  m^g  is  the  mass  of  a  | 

2  * 

cesium  atom.  For  our  conditions,  Dcs-N2  is  on  the  order  of  0.1  cmVsec,  i 

euid  taking  the  characteristic  dimension  for  density  gradients  to  be  one  j 

centimeter,  maximum  diffusion  velocities  are  on  the  order  of  0.1  cm/ sec. 

gn  i 


B.2.1  Flowing  System  Desi 


from  the  apparatus.  A  minimum  flow  rate  was  chosen  to  sweep  out  impuri¬ 
ties  which  might  have  entered  the  system  at  a  rate  corresponding  to  the 
measured  vacuum  leak  rate.  The  residence  time  in  the  region  of  the  cell 
where  recombination  was  to  be  observed  needed  to  be  much  greater  than 
the  predicted  two -millisecond  reaction  time,  so  the  maximum  flow  rate 
through  the  2-cm-diameter  reaction  cell  was  limited  to  10  liters/min. 

The  tightest  constraint  on  gas  flow  rate  was  calculated  to  be  the 
time  required  to  saturate  nitrogen  bubbles  with  cesium  in  the  reservoir. 
This  calculation  was  based  on  the  time,  necessary  for  saturated 
cesium  vapor  to  diffuse  from  the  outer  boundary  of  a  nitrogen  bubble  of 
radius  r  to  the  center.  Diffusion  was  considered  complete  when  the 
cesium  number  density  at  the  bubble  center  was  96%  of  the  value  at  the 
bubble  boundary.  The  diffusion  time  was  given  by  [1.54] 

^diff  =0*4  i^^/Dcs-N2  * 

The  residence  time,  of  the  bubble  in  the  cesium  was  found  by  solv¬ 
ing  the  bubble's  equation  of  motion  for  its  terminal  velocity  [1.55], 
and  by  estimating  a  1-cm  path  through  the  cesium.  The  bubble  diameters 
expected  from  the  0.020"  holes  in  the  nitrogen  nozzle  (Fig.  1.16)  were 
in  the  range  0.04"-0.1".  These  bubbles  were  found  to  have  Reynolds  num¬ 
bers  of  about  1000,  and  a  pressure  drag  coefficient  of  0.4  was  estimated 
[1.55].  The  0 . 1 "-diameter  bubbles  are  the  largest  acceptable  due  to  the 
requirement  _>  1.  This  ratio  is  proportional  to  r“^’^.  Gas 
flow  rates  were  limited  by  the  maximum  diameter  bubble  and  minimum  resi¬ 
dence  time  to  less  than  40  std  cc/min.  Convective  transport  of  vapor 
through  the  bubble  was  not  considered  here,  but  it  could  be  a  dominant 
transport  mechanism  [1.55].  Hence  the  flow-rate  limit  set  by  this  bub¬ 
ble  saturation  calculation  is  thought  to  be  conservative. 

Aerosols  of  liquid  cesium  were  found  to  be  generated  by  the  bubbling 
in  the  cesium  reservoir.  The  aerosols  were  detected  with  scattering 
from  a  laser  bean^.  They  would  control  the  cesium  vapor  number  density 
based  on  the  cell  temperature,  and  the  reservoir  ten^rature  could  not 
be  used  as  a  direct  indication  of  number  density.  In  order  to  prevent 
the  aerosols  from  leaving  the  reservoir,  the  reservoir  diameter  above 


the  cesium  was  Increased  three-fold  to  1.5",  and  packed  copper  wool  and 
stainless  steel  mesh  were  inserted  1.5"  above  the  cesium  level.  Trap¬ 
ping  in  this  manner  was  not  effective. 


To  avoid  the  aerosol  problem,  an  alternative  approach  would  be  to 
replace  bubbling  with  evaporation  from  a  wick.  The  wick  serves  to 
increase  the  ratio  of  surface  area  to  volume  for  the  liquid  cesium, 
thereby  increasing  the  evaporation  rate.  Capillary  action  limits  the 
wick  vertical  dimension.  Metal  wicks  have  been  used  successfully  for 
evaporation  in  alkali  metal  heat  pipes  [1.49].  In  a  heat  pipe  oven  the 
wick  is  adjacent  to  the  reaction  zone;  the  wick  may  also  be  placed  in  a 
separate  reservoir  for  a  flowing  syst«n.  Diffusion  and  natural  convec¬ 
tion  transport  vapor  in  a  heat  pipe,  so  bulk  velocities  are  low  and  wick 
size  can  be  smaller  than  in  a  flowing  system  which  uses  forced  convec¬ 
tion  with  higher  bulk  velocities. 

B.2.2  Cesium  Number  Density  Measurement 

Ideally  the  value  of  the  cesium  number  density  in  the  reaction  cell 
should  be  given  from  equilibrium  and  saturated  vapor  pressure  relations 
[1.39]  as 
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Here  and  Tj.gg  are  respectively  the  temperatures  (in  Kelvins)  of 

the  vapor  in  the  reaction  zone  and  the  hottest  portion  of  cesium  liquid 
in  the  system.  (The  units  of  number  density  are  cm"^.)  The  actual  num¬ 
ber  density  might  vary  from  this  equilibrium  prediction  because  of  cold 
spots,  hot  spots,  aerosols,  adsorption,  desorption,  or  convection.  To 
test  that  a  system  is  performing  as  expected,  a  direct  measurement  of 
cesium  concentration  is  required.  Possible  techniques  for  measuring  the 
cesium  number  density,  n^g,  can  be  divided  into  intrusive  and  non- 
intrusive  categories.  The  intrusive  techniques  require  collection  of 
the  Cs  vapor.  Intrusive  techniques,  which  we  have  looked  into,  are 
based  on  the  density,  the  chemical  properties,  and  the  dew  point  of 
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cesium.  Non-intrusive  techniques,  which  we  have  used,  depend  on  the 


optical  properties  of  the  vapor. 


The  density-based  technique  measures  the  weight  of  Cs  collected  in  a 


cold  trap  in  a  flowing  system  over  a  known  time  period.  The  chemistry 


method  reacts  collected  Cs  with  alcohol  to  form  CsOH  and  then  titrates 


the  CsOH  with  HCl  to  determine  the  eunount  of  Cs  collected.  The  dewpoint 


method  varies  the  temperature  of  a  specular  surface  which  is  in  contact 


with  the  Cs  vapor.  When  the  surface  temperature  equals  the  dew  point, 
condensation  can  be  observed  (e.g.,  with  em  increase  in  diffuse  scatter¬ 
ing  of  a  laser  beam) .  The  frequency  response  of  this  technique  will 
decrease  at  higher  N2  pressures  euid  at  lower  n^.^. 


The  three  optical  techniques  mentioned  below  involve  observing  a  Cs 
absorption  line  using  a  tungsten  halogen  lamp,  a  monochromator,  and  a 


photomultiplier  tube.  These  techniques  are  relatively  fast  and  can  be 
used  for  real-time  monitoring  of  n^g.  The  range  of  Cs  number  density 
over  which  the  technique  is  useful  depnds  on  the  absorption  path  length. 


the  monochromator  resolution,  and  the  line  chosen.  The  technique  of 


matching  actual  and  calculated  line  shapes  requires  both  good  spectral 
resolution  and  knowledge  about  broadening.  The  equivalent-width  tech¬ 
nique  needs  knowledge  about  the  broadening  but  less  resolution,  and  the 
integrated-absorption  technique  requires  good  spectrometer  resolution 


but  no  knowledge  of  broadening  to  determine  n^g. 


For  our  conditions,  the  absorption  line  shape  was  dominated  by 


nitrogen  pressure  broadening.  In  this  case  the  line  shape  is  Lorentzian 


and  the  absorption  cross  section  is  [1.56] 


0(0,)  =  - li  ( - ^ ^ 
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Here  f^2  denotes  the  oscillator  strength  for  absorption  from  the  ground 
to  the  first  excited  level  of  Cs,  o)  is  the  angular  frequency  (radians/ 


sec),  e  is  the  charge  of  one  electron,  m^  is  the  mass  of  an  electron. 


c  is  the  speed  of  light,  e-  is  the  permittivity  of  free  space,  and  Au; 
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The  value  of  Aai  for  N2-CS  broadening  was  inferred  from  Miles  [1.57] 
and  Chen  [■'.58].  A  good  source  for  oscillator  strengths  is  Fabry 
[1.59].  When  this  cross  section  is  used  in  Beer's  law, 

I(a),L)/I(tij,0)  =  exp[-n^(J(io)L]  , 

the  width  of  the  absorption  profile,  l(a>,L)  can  be  much  greater  than 
Ab).  Here  L  denotes  the  path  length  in  meters,  1(0), x)  is  transmitted 
light  intensity  at  the  specified  frequency  and  position,  euid  n^  denotes 
the  ground  state  cesium  number  density. 

The  expression  for  the  absorption  cross  section  can  be  used  to  match 
the  measured  absorption  profile  with  a  theoretical  profile.  A  sophisti¬ 
cated  routine  could  optimize  both  the  parameters  n^  and  Au)  for  the 
fit.  A  fit  is  especially  useful  in  checking  for  anomalies  due  to  aero¬ 
sols,  dimer  band  absorption,  or  systematic  errors. 

The  expression  for  the  cross  section  can  be  used  to  evaluate  an 
equivalent  width  defined  as  [1.56] 

W  =  /  {1  -  exp(-n  a{a))L]  }du)  . 

line 

Evaluating  the  equivalent  width  in  the  optically  thin  and  optically 
thick  limits  gives  expressions  for  the  ground  state  Cs  number  density. 
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Between  these  two  limits  a  curve  of  growth  technique  is  necessary  to 
establish  the  relationship  between  ni»  Au,  emd  W.  The  equivalent  width 
is  simply  a  measure  of  total  absorption,  and  it  can  be  obtained  directly 
with  a  monochromator  of  resolution  which  is  coarse  relative  to  the  Cs 
line.  This  technique  avoids  the  need  to  scan  the  absorption  line,  but 
it  is  susceptible  to  errors  due  to  effects  which  would  change  the  ab¬ 
sorption  profile  from  the  expected  one. 


By  integrating  Beer's  law  over  the  frequency  range  contributing  to 
the  line,  it  can  be  shown  [1.60]  that  the  ground  state  Cs  number  density 


is  given  by 
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This  integration  approach  requires  a  means  to  integrate  the  spectral 
transmission,  but  no  information  about  the  line  shape  is  needed. 

B.3  Light  Sources 

Photoionization  of  the  Cs  by  Xe  flashlamp  radiation  occurs  through 
both  one-  and  two-proton  processes.  The  single-photon  process  converts 
a  ground  state  6s  electron  to  a  free  electron  as  represented  by 


Cs(6^S^/25  ® 

The  photon  energy,  hv^,  must  exceed  the  ionization  potential  of  Cs, 
i.e.,  3.89  eV  or  318  nm.  Two-photon  ionization  first  excites  a 

ground  state  electron  with  a  resonant  photon;  a  second  photon  then 
ionizes  the  excited  atom.  The  two  steps  are 


Cs(6^S./j)  +  hv- 


Cs  (6  P3/2,1/2) 


Cs*  (6^P2/2, 1/2^  ^^^3  ® 


The  resonance  wavelengths  2u:e  X2  =  852  nm  (J  =  3/2)  or  X2  =  894  nm  (J  = 
1/2),  and  ionization  from  the  6p  state  can  occur  with  any  photon  of 
wavelength  X^  ^  508  nm  ( J  =  3/2)  or  X^  ^  494  nm  ( J  =  1/2) .  The  cross 
sections  for  ionization  from  the  6s  and  6p  levels  are  approximately 
10  cm  and  10  cm  ,  respectively.  The  single-photon  process  has 
been  included  in  our  Icinetics  calculations;  the  ionization  rate  is  pro¬ 
portional  to  the  light  intensity.  The  contribution  to  ionization  rate 
from  the  two-photon  process  was  estimated  to  be  significant;  this  con¬ 
tribution  is  proportional  to  the  square  of  light  intensity  [1.61]. 


The  experimental  arrangement  of  Fig.  1.7  has  the  sin5)lifying  feature 
that  the  light  source  is  positioned  away  from  the  heat  of  the  reaction 
zone.  Placing  the  source  outside  the  cell  necessitated  directing  the 
light  through  a  sapphire  window  at  the  end  of  the  Inconel  tube.  To  pro¬ 
duce  a  well-collimated,  intense  photon  flux,  we  used  a  5-atm  xenon  short 
arc  lamp  and  an  ellipsoidal  reflector. 

The  flashlamp  dissipated  a  maximum  of  300  Joules  in  4  ms,  corre¬ 
sponding  to  an  average  electrical  input  power  of  75  kw.  A  circuit  dia¬ 
gram  of  the  pulsed  power  supply  is  shown  in  Fig.  1.17.  The  circuit 
design  was  determined  by  two  principal  factors:  the  very  low  lan^  impe¬ 
dance  combined  with  the  need  to  keep  the  pulse  length  no  longer  than  a 
few  milliseconds,  and  the  requirement  of  a  high-voltage  (30  kV)  trigger 
pulse  to  initiate  the  discharge.  The  current  pulse  trace  of  the  final 
circuit,  shown  in  Fig.  1.18,  was  obtained  with  0.010F  of  storage  capac¬ 
itor  charged  to  375  V.  The  peak  current  was  2.2  kA  and  the  total  pulse 


Fig.  1.17.  Fla.shlamp  circuit. 
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Fig.  1.18.  Flashlamp  pulse  shape.  Broadband  light  output  follows 
the  current. 


length  was  4  ms.  The  energy  input  could  be  varied  by  changing  the  volt¬ 
age  of  the  DC  power  supply  used  to  charge  the  capacitor.  Although  the 
kinetics  calculations  indicate  the  need  for  a  shorter  pulse  length,  the 
present  power  supply  was  deemed  to  be  satisfactory  for  performing 


exploratory  experiments. 


B .  3 . 2  Linear  Lancs 


The  low  impedance  of  the  short  arc  lan^)  resulted  in  a  pulse  duration 
which  was  about  as  long  as  the  expected  recombination  times.  A  shorter 
pulse  can  be  produced  by  crowbarrinq  the  short  arc  lan^)  or  by  using  a 
linear,  low-pressure  lamp  with  greater  lamp  impedance.  The  linear  lamp 
typically  would  have  a  length  of  10  cm  and  a  diameter  of  0.4  cm.  Light 
vrauld  be  produced  throughout  this  volume.  Light  from  the  linear  lamp 
could  be  coupled  to  a  reaction  cell  by  placing  the  lamp  parallel  to  and 
outside  a  cylindrical  cell  with  transmissive  walls  (e.g.,  a  Lucalox 
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tube),  or  by  placing  the  lamp  coaxially,  inside  a  cell  of  arbitrary  wall 
material.  The  former  arrangement  could  employ  commercially  available 
flashlamps  which  must  be  operated  below  400®C  to  protect  their  glass-to- 
metal  seals.  The  lamps  would  be  cooled  by  em  annular  flow  of  air 
between  the  lamp  emd  a  quartz  tube  sheath.  Experiments  in  our  labora¬ 
tory  have  demonstrated  the  feasibility  of  such  cooling.  The  latter 
arrangement  would  require  a  lamp  which  could  withstand  the  high  tempera¬ 
ture  and  alkali  environment  of  the  reaction  cell. 

A  lamp  capable  of  serving  in  the  coaxial  configuration  has  been 
built  and  tested.  A  lucolox  tube,  4  nm  I.D.  and  26"  long,  was  fitted 
with  tungsten  electrodes  which  formed  a  4"  central  gap.  The  lamp  ends 
were  sealed  with  o-rings  and  the  lamp  was  filled  with  200  torr  of  argon. 
The  lamp  circuit  was  similar  to  the  short  arc  circuit  of  Fig.  1.17,  ex¬ 
cept  that  the  discharge  capacitor’s  characteristics  were  5  yF  and  10  kV, 
the  trigger  transformer  saturated  inductance  was  18  yH,  and  no  lumped 
series  resistance  was  used.  The  lamp  was  tested  in  an  oven  with  a 
center  temperature  of  1050®C  and  cold  ends.  The  capacitor  energy  was  22 
Joules  and  the  pulse  duration  was  30  ysec.  This  lamp  had  a  pulse  dura¬ 
tion  much  shorter  than  the  expected  recombination  times.  In  the  case  of 
single-photon  ionization,  the  final  electron  number  density  depends  on 
the  Icunp  energy,  which  was  a  factor  of  12  smaller  for  the  linear  lamp 
compared  to  the  short  arc.  In  the  case  of  two-photon  ionization,  the 
final  electron  number  density  depends  on  the  product  of  the  square  of 
lamp  intensity  and  pulse  duration.  This  product  was  about  equal  for  the 
two  lamps  we  have  built. 

B.4  piagnostics 

The  number  density  of  reactants  (electrons  or  cesium  ions),  or 
alternatively,  the  number  density  of  the  products  (cesium  atoms)  must  be 
measured  in  order  to  determine  the  recombination  reaction  rate.  In  this 
section  various  pertinent  diagnostic  techniques  are  discussed.  These 
techniques  aure  based  on  electromagnetic  wave  propagation  through  the 
plasma,  determination  of  the  plasma  electrical  impedance  or  optical 
spectroscopy.  The  cesium  ions  are  generally  not  useful  for  diagnostic 
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purposes  because  they  have  a  small  effect  on  the  electrical  characteris¬ 
tics  of  the  plasma  compared  to  the  electrons/  and  the  ions  are  difficult 
to  access  spectroscopically.  The  plasma  to  be  diagnosed  oem  be  charac¬ 
terized  by  the  following  representative  conditions:  a  constant  nitrogen 

number  density/  nj^  z  of  10  cm  ,  and  initial  cesium  number  density/ 

1  ^  8 
”Cs'  10''^  cm“^,  and  an  afterglow  electron  number  density/  n^i  of  10 

to  10^^  cm”^;  the  electron  collision  frequency,  v,  is  1.3x10^^  sec“^  and 

independent  of  temperature  [1.62]/  and  the  inverse  of  the  characteristic 

time  for  recombination  on  nitrogen  is  10^  sec”S  the  temperature  is 

1000®C  and  the  plasma  is  contained  in  a  cylinder  of  5-cm  diameter  and 

10-cm  length.’  Values  of  plasma  frequency/  oip  =  (ngC  '  ,  and 

d.c.  electron  electrical  conductivity,  Oq  =  UgO  /ung,  are  listed  in 

Table  1.5  for  the  broadest  anticipated  range  of  n^  during  recombination. 


Table  1.5 

Plasma  Afterglow  Electrical  Characteristics 

10®  10®  10^®  10^’  10^^ 
5.6x10®  1.8x10®  5.6x10®  1.8x10^®  5.6x10' 

2.2x10"^  2.2x10""*  2.2x10"®  2.2x10"^  2.2x10" 


(Up,  sec 
a„, mho/m 


B.4.1  Microwave  Diagnostics 

The  plasma  dimensions  must  be  much  larger  than  the  wavelength  of 
electromagnetic  radiation  which  is  used  to  diagnose  the  plasma  if  plane 
wave  propagation  is  assumed.  A  technique  which  uses  a  resonant  cavity 
generally  requires  the  plasma  dimensions  to  be  on  the  order  of  a  wave¬ 
length.  If  the  wavelength  range  of  useable  radiation  is  10  cm  (the 
length  of  the  plasma  cylinder)  or  shorter,  the  corresponding  angular 
frequency  range  is  (u  1.8x10*®  sec"*.  Plane  electromagnetic  waves 
propagating  in  the  z-direction,  in  a  uniform  medium  have  the  phase  fac¬ 
tor  exp [ j ((ut-6z)  “  az]  .  Here  a  amd  B  are  respectively  the  attenuation 
and  phase  coefficients  which  are  given  by  the  following  expressions 
[1.63]  : 
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For  the  plasma  conditions  of  Table  1.5,  it  is  generally  true  that  >> 
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01^.  In  this  limit 
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The  plane  wave  will  be  altered  by  the  plasma  if 
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This  condition  is  satisfied  only  when  o)  <  v  ,  and 
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The  last  relation  in  conjunction  with  Table  1.5  establishes  the  range  of 


useful  probe  frequencies,  2.4x10®  <  o)  <  2.4x10^®  sec”  .  Recall  that  the 


size  of  the  plasma  raises  the  lower  limit  on  probe  frequency  to  1.8x10 


sec”  .  Therefore  the  range  of  electron  number  density  which  can  be 
diagnosed  by  electromagnetic  wave  propagation  diagnostics,  based  on  mea¬ 


surable  attenuation  or  phase  shift  and  plasma  size,  is  7.5x10''  <  ng  < 


1.0x10''^.  Only  a  small  part  of  the  anticipated  range  of  n^  (see  Table 


1.5)  is  included,  and  the  included  range  is  too  small  to  measure  appre¬ 


ciable  recombination. 
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The  electrical  character  of  the  plasma  afterglow  is  dominated  by  the 
electrons;  ions  play  a  negligible  role.  Applied  electromagnetic  fields 
and  collisions  with  neutrals  control  the  behavior  of  the  electrons. 
Assuming  the  applied  fields  are  electric,  and  proportional  to  exp(jcjt), 
the  complex  plasma  electrical  conductivity  is  given  by  [1.63] 


e  1  -  jo)/v 

m  V  .  ^  ,  ,  .2 
e  1  +  (m/v) 


When  the  wavelength  of  the  applied  field  is  much  longer  than  the  length 
of  plasma  filling  the  space  between  opposed  electrodes  across  which  the 
field  is  applied,  the  plasma  may  be  treated  as  a  lumped  impedance.  Tak¬ 
ing  the  electrode  spacing  to  be  one  centimeter,  an  impedance  measurement 
can  be  made  when  oi  <<  1.8x10^^  sec~^.  If  an  impedance  measurement  is 
made  in  the  afterglow,  where  v  =  1.3x10^^  sec“^,  then  u)  <<  v,  and  the 
plasma  conductivity  is  essentially  the  d.c.  conductivity,  =  UgC  /mgV. 
Measurement  of  the  d.c.  conductivity  yields  information  about  the  ratio 
Ug/v,  but  data  are  available  [1.62]  to  enable  independent  determination 
of  the  electron  neutral  collision  frequency.  The  electrical  circuit 
which  measures  the  plasma  conductivity  may  be  directly,  capacitively,  or 
inductively  coupled  to  the  plasma. 

Electrodes  contact  the  plasma  in  a  directly  coupled  conductivity 
measurement.  Account  must  be  taken  of  the  sheath  impedance  between  the 
electrode  and  the  plasma.  Although  theory  exists  for  spherical  Langmuir 
probes  in  collision-dominated  plasmas  [1.64],  the  accuracy  of  the  theory 
for  such  probes  is  still  in  question  [1.65].  To  avoid  calculation  of 
the  sheath  impedance,  a  four-electrode  probe  was  fabricated  and  tested. 
As  shown  in  Fig.  1.19,  the  probe  consists  of  a  linear  2urray  of  three 
electrically  insulated,  cylindrical  (wire)  electrodes  with  a  short 
length  of  the  free  ends  uninsulated  and  in  contact  with  the  plasma.  The 
inside  surface  of  the  reaction  cell  wall  serves  as  the  fourth  electrode. 

The  design  and  operation  of  the  probe  have  been  described  in  an 
earlier  report  [1.66].  The  basic  idea  is  that  a  measured  current,  I,  is 
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Fig.  1.19.  Conductivity  probe  electrode  configuration. 


driven  through  the  plasma  between  the  outer  electrodes  consisting  of  the 


pin  on  the  axis  of  the  cell  and  the  reaction  cell  wall,  while  the  volt¬ 


age  6V  picked  up  between  the  two  off-center  inner  pins  is  recorded.  By 
measuring  the  voltage  between  the  inner,  floating  pins,  sheath  drops  can 
be  canceled,  and  the  probe  interpretation  is  simplified.  A  detailed 
schematic  of  the  probe  circuitry  is  shown  in  Fig.  1.20. 


Given  the  concentric  arrangement  of  the  current-carrying,  outer 


electrodes,  the  voltage  difference  V  =  V(a)  -  V(b)  between  the  outer 


electrodes  can  be  calculated  (assuming  plasma  resistance  dominates 


sheath  drops,  and  neglecting  end  effects)  from  the  relation 


b  b 

-  /  dV(r)  =  /  pJ  dr  . 
a  a 


Here  a  and  b  are  the  radii  of  the  on-axis  and  outer  electrodes,  respec¬ 


tively,  p  is  plasma  resistivity  (ohm  meters),  J  is  current  density 


(amps/m^),  and  r  is  the  radial  coordinate.  The  result  of  the  integra¬ 


tion  is 


V(a)  -  V(b) 
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Fig.  1.20.  Conductivity  probe  circuitry. 


*■ 


The  voltage  difference  6V  between  the  inner  electrodes  is  given  by  a 
similar  expression.  Inserting  values  for  the  actual  radii,  the  loca¬ 
tions  of  the  inner  electrodes,  and  the  electrode  length  L,  the  following 
values  for  the  effective  resistances  (in  ohms)  between  the  electrodes 
are  obtained: 


iY 

31 


22.0  p; 


35V 

31 


4.34  p  . 


With  reference  to  Fig.  1.20,  let  R  represent  the  sum  series  resis¬ 
tance  of  the  1  Mfi  potentiometer  and  the  10  kO  or  30  0  current  resistors, 
let  Rpj^asma  represent  dV/dl,  and  let  E  represent  the  voltage  applied  to 
the  outer  electrodes  (assumed  to  be  constant).  By  differentiating  the 
circuit  equation. 


the  effect  of  a  change  in  the  plasma  resistance  ( such  as  a  sudden  de¬ 
crease  during  photoionizaticn)  on  the  measurable  parameters  V,  6V,  and  I 
can  be  predicted: 
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The  probe  cem  be  rvin  under  two  limiting  regimes*  In  the  constant- 
voltage  regime,  R  <<  Rpi^sma  R  is  set  for  30  fi) ,  and  only  the 

current  is  affected  by  a  change  in  Rpiastaa’  constant-current 

regime  the  inequality  is  reversed  (generally  R  is  set  to  1  MJO  ,  and  only 
V  and  6V  are  affected  by  a  change  in  I^iasma* 

Testing  of  the  four-electrode  probe  has  not  been  completed,  but  some 
general  observations  can  be  stated.  The  capacitive  susceptance,  ojC, 
between  electrodes  should  be  kept  much  less  than  the  conductance,  1/R, 
so  that  the  probe  measures  conduction,  not  displacement  currcr  ,*  This 
condition  applied  to  two  electrodes  in  a  uniform  plasma  is 


0)  «  1/(RC)  =  /e  =  (D  /v  * 
dc  o  p 

For  the  lowest  ojp  of  Table  1.5,  this  condition  is  met  when  w  <<  2x10® 
sec”^.  For  the  lowest  of  Table  1.5,  the  signal  source  impedance  of 
the  probe  was  36V/3I  =  2x10^  n.  This  large  soxxrce  impedance  could  limit 
frequency  response  and  induce  noise  without  careful  circuit  design.  The 
high  source  impedance  required  good  electrical  insulation  of  the  elec¬ 
trodes,  which  was  difficult  to  achieve  in  the  high-temperature,  liquid 


metal  environment. 


Recombination  on  the  probe  surface  could  create 


local  plasma  conductivity  non-uniformity.  Photoemission  from  the  elec¬ 
trode  surfaces  during  photoionization  may  significantly  affect  the 
measurement. 

A  capacitively  coupled  measurement  operates  at  a  frequency  where  the 
capacitive  susceptance  of  the  sheaths  is  large  compared  with  the  total 
admittance  of  the  electrodes,  i.e.,  »  uC  +  VRpj^ggnja*  The 
electrodes  can  be  separated  from  the  plasma  by  a  dielectric  wall  if  this 


Inequality  holds  when  the  sheath  capacitance  Is  replaced  by  the  series 
capacitance  of  the  wall  and  sheath,  l.e.,  Cgj^g^^h 

^  ^sheath  Sijall^’  case  of  parallel  plate  electrodes  of  area  A 

and  separation  d,  where  Impedance  Is  dominated  by  the  Inter-electrode 
capacitance,  C,  auid  uniform  plasma  conductivity,  current  through 

the  electrodes,  1  exp(jajt),  due  to  an  applied  voltage,  V  exp(jait).  Is 
given  by  I  =  ( VA/d)  ( jaie^  +  cr^j^).  A  high-frequency  bridge  can  be  used  to 
separate  the  conduction  and  displacement  currents  to  determine  the 
plasma  conductivity  [1.67] .  In  the  actual  circumstance  of  a  non-uniform 
plasma  conductivity,  the  current  Is  given  by 


1  =  va[  / 
0 


dx 


{ jtoe  + 
o 


dc 


t1 


-1 


The  plasma  conductivity  spatial  dependence  determines  If  the  elec¬ 
trode  conduction  current  Is  proportional  to  the  average  plasma  conduc¬ 
tivity.  Conductivity  In  the  afterglow  Is  controlled  by  volume  and 
surface  recombination,  and  arablpolar  diffusion.  We  have  mathematically 
modeled  the  spatial  and  temporal  behavior  of  the  afterglow  conductivity 
field,  and  we  have  used  the  result  to  numerically  integrate  the  expres¬ 
sion  for  I.  The  integral  always  yields  proportionality  between  conduc- 
tion  current  and  conductivity  in  the  limit  (a>eQ/aj|p)  >>  1.  For  the 
largest  conductivity  listed  in  Table  1.5,  this  limit  would  require  co  >> 
2x10^®  sec“^,  but  our  numerical  calculations  indicate  that  this  is  an 
unnecessarily  extreme  constraint  on  u. 

The  Impedance  measuring  circuit  may  also  be  coupled  to  a  cylindrical 
plasma  inductively  through  ^uclal  magnetic  fields  [1.68],  or  azimuthal 
magnetic  fields  [1.69].  Inductive  coupling  may  have  some  advantages 
over  capacitive  coupling  in  working  at  lower  applied  field  frequencies, 
and  in  working  with  a  conductive  plasma  container. 

B.4.3  Optical  Spectroscopy  Diagnostics 


Optical  diagnostics  are  attractive  because  the  diagnostic  hardware 
can  usually  be  located  outside  the  harsh  environment  of  the  reaction 
cell,  and  optical  access  to  the  afterglow  plasma  often  can  be  compatible 


with  access  for  photoionization.  Line  and  continuum  emission,  and 
atomic  absorption  have  been  considered  as  the  basis  for  measurement  of 
electron  number  density  and  temperature  in  our  plasma. 

The  specific  power,  P^_j^  (watts/m^),  of  line  emission  from  an  upper 
state,  u,  to  a  lower  state,  £,  of  the  cesium  atom  is  proportional  to  the 
upper  state  nundser  density,  n^,  and  can  be  written  [1.70] 

^u-Jl  ~  ' 


where  is  the  Einstein  coefficient  for  spontaneous  emission,  and 

hv^_l^  is  the  energy  difference  between  the  states.  If  the  upper  state 
is  chosen  to  be  close  to  the  continuum,  the  population  of  that  state 
will  be  in  equilibrium  with  the  free  electron  population  throughout 
recombination.  Assuming  the  ion  and  electron  number  densities  are 


equal,  n^  and  n^  are  then  related  by  the  Saha  equilibrium  relation. 
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A  >•  /»  ^ 
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where  g  denotes  the  degeneracy  of  the  cesium  ion  or  excited  state, 
is  the  energy  difference  between  the  excited  state  u  and  the  ground 
state  ion,  and  T  is  the  equilibrium  temperature  of  the  system.  Colli¬ 
sions  with  nitrogen  rapidly  bring  the  photo -electrons  to  the  equilibrium 
teaperature.  An  absolute  intensity  measurement  of  line  emission  can  be 
related,  by  the  preceding  equations,  to  electron  number  density. 


The  specific  power,  Pq_£,  of  continuum  emission  accompanying  capture 
of  free  electrons  of  speed  distribution  Fg(Ug)  to  a  bound  state,  £,  in 
the  cesium  atom  in  the  frequency  range  to  V2  is  given  by  [1.71] 


P 

e-£ 


U2 

/  hv 


I  a(u  )u  F  (u  )du 
e-£  e  e  e  e  e 


where  Ug  is  the  free  electron  velocity,  hVg_j^  is  the  kinetic  plus  poten¬ 
tial  energy  difference  between  the  free  and  bound  electron,  and  0(Ug)  is 
the  cross  section  for  radiative  capture.  Measurement  of  the  oontinuxim 
emission  can  yield  information  about  the  electron  number  density  and  the 
electron  velocity  distribution.  If  Fg  is  a  Maxwellian  distribution,  the 
electron  tenperature  can  be  determined. 
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Diagnostics  based  on  line  and  continuum  emission  share  several  com¬ 
mon  features.  Since  both  forms  of  emission  are  present  in  the  spectrum 
simultaneously/  care  must  be  taken  to  chose  regions  of  observation  where 
overlap  is  minimized.  Additional  interference  o£m  come  from  the  overlap 
of  continuum  bands.  Although  the  photoionization  source  should  be 
extinguished  during  recombination  measurements,  the  optical  detector 
must  be  protected  from  d^unage  during  photoionization.  These  emission 
measurements  require  absolute  intensity  calibration;  relative  intensity 
measurements  can  determine  temperature,  but  not  electron  number  density. 
The  measured  intensity  is  proportional  to  the  square  of  electron  number 
density  when  ion  and  electron  number  density  are  assumed  to  be  equal. 
The  ratio  of  signal  to  shot-noise  is  proportional  to  the  square  root  of 
the  signal.  To  observe  a  factor  of  10  decay  in  n^,  a  factor  of  10 
decrease  in  signal-to-noise  must  be  tolerated.  At  low  electron  number 
density  the  emission  signal  at  the  start  of  recombination  is  near  the 
shot  noise  detection  limit  and  further  signal  reduction  may  not  be 
acceptable.  For  example,  we  have  calculated  the  signal-to-shot  noise 
ratio  for  continuum  emission  in  the  4700  A  to  5000  A  range  due  to  recom¬ 
bination  to  the  6P  level  of  cesium,  and  for  line  emission  at  6326  A  due 
to  the  10F  to  5D3/2  transition.  The  regions  of  observation  were  chosen 
to  minimize  line  and  cont^.’^aum  overlap  (1.72].  With  the  assumptions  of 
an  emitting  volume  of  20  cra^,  2%  optical  collection  efficiency,  25% 
detector  quantum  efficiency,  10“^  sec  temporal  resolution,  and  ng  = 
3x10^®  cm”^,  the  signal-to-noise  ratio  for  the  continuum  and  line  emis¬ 
sion,  respectively,  -.^ere  55  and  245.  Careful  optical  design,  perhaps 
utilizing  filters  for  spectral  selectivity  with  large  etendue,  is 
required  to  use  emission-based  techniques  to  diagnose  the  plasma  condi¬ 
tions  of  Table  1.5. 

Cesium  atomic  absorption  has  been  used  extensively  in  our  reservoir 
development  efforts.  When  substantial  photoionization  of  the  ground 
state  of  cesium  occurs,  a  corresponding  reduction  in  absorption  by  the 
ground  state  will  result.  It  has  been  shown  [1.73,  1.74]  for  conditions 
of  our  plasma  (where  the  number  density  of  excited  atoms  is  much  less 
th^ul  the  number  density  of  free  electrons,  and  the  mean  thermal  energy 
is  much  less  than  the  first  excitation  energy)  that  the  relaxation  times 


for  the  excited  levels  are  much  shorter  than  the  relaxation  time  for  the 
ground  level  or  for  the  free  electrons.  This  display  in  relaxation 
times  {i.e.,  partial  equilibrium)  allows  us  to  equate  the  negative  of 
the  time  derivative  of  the  ground  state  population  with  the  time  deriva¬ 
tive  of  the  electron  population.  Thus  the  recombination  rate  can  be 
determined  by  observation  of  ground  state  absorption,  provided  the  frac¬ 
tional  change  in  the  ground  level  population  is  not  too  small.  To  exam¬ 
ine  the  feasibility  of  this  approach,  a  photoionization  experiment  was 
performed  to  measure  the  fractional  degree  of  ionization,  under  the  same 
assumptions.  For  this  experiment,  the  monochromator  was  kept  at  a 
single  spectral  location,  as  shown  in  Fig.  1.21,  and  the  change  in 
absorption  during  photoionization  was  monitored.  A  flashlamp  served  as 
the  source  both  for  the  absorption  measurements  and  for  photoionization. 

The  data  were  interpreted  by  inverting  Beer's  law  to  get  the  follow¬ 
ing  expression  for  ground  state  cesium  number  density  as  it  varied  with 
time, 

n,(t)  =  {fn[Io(t)/I(t)]}/[La(u)]  . 

Here  L  is  the  optical  path  length,  assumed  to  be  constant  during  the 
flash,  o(a))  is  the  optical  cross  section  for  absorption,  assumed  not  to 
be  a  function  of  time  or  location.  The  quantity  I(t)  is  the  time- 
varying  transmitted  source  spectral  intensity  amd  lQ(t)  is  the  source 
spectral  intensity  without  absorption.  Note  that  IqC^)  also  varied  with 
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Fig.  1.21.  One-beam  ionization  experiment. 


time.  By  manipulating  this  expression,  the  following  expression  for 
relative  change  in  number  density  is  derived: 

n(t^)  -  n(t2)  An[I(t2)/lQ(t2)] 

^TtT)  ''  ■  An[I(tJ/I  (t,)]  • 

1  1  O  I 

Here,  t^  is  typically  the  time  immediately  preceding  photoionization, 
and  t2  refers  to  any  time  during  which  the  cesium  was  not  in  equilib- 


First  an  absorption  scan  was  taken  with  the  tungsten  lamp  and  the 
spectral  locations  for  transmission  fractions  equal  to  0.2,  0.8,  and  1.0 
were  identified.  Then  the  tungsten  halogen  lan^)  was  replaced  by  the 
xenon  short-arc  flashlanp  and  the  cesium  vapor  was  photoionized.  The 
light  which  passed  through  the  vapor,  during  the  flash,  was  monitored  at 
a  single  wavelength  by  the  monochromator.  The  data  indicated  that 
ground  state  depletion  occurred,  but  scatter  due  to  pulse -to-pulse  vari¬ 
ation  of  the  flashlanp  output  made  the  results  inconclusive. 

Two-beam  experiments  which  utilize  orthogonal  photoionization  and 
probe  beams  to  remove  the  measurement  dependence  on  pump  beam  variations 
were  also  performed.  No  indication  of  ground  state  depletion  was  ob¬ 


served  in  these  experiments. 


However,  the  interpretation  of  these 


ground  state  absorption  measurements  were  con5)licated  by  absorption  in 
cesium  vapor  which  is  not  within  the  photoionizing  beam,  and  by  attenua¬ 
tion  of  the  photoionizing  beam  before  it  crossed  the  probe  beam. 

B.5  Conclusions 

The  work  described  above  leads  to  the  following  conclusions  with 
respect  to  the  design  of  a  new  experiment.  Reaction  cell  materials 
could  be  polycrystalline  alumina,  stainless  steel,  or  Inconel.  The 
metals  should  be  plated  on  the  interior  with  nickel  to  prevent  the  out- 
gassing  of  alloys.  These  naterials  are  relatively  inexpensive,  avail¬ 
able,  oxidation -resistant,  and  easy  to  use,  in  addition  to  being  compat¬ 
ible  with  cesium  vapor  to  1200®C.  Window  naterial,  to  separate  the 
cesium  vapor  and  the  photionization  source,  could  be  fused  silica,  Luca- 
lox,  or  single -crystal  alumina.  The  latter  two  materials  could  with¬ 
stand  the  reaction  cell  environment  while  the  silica  must  be  protected 
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with  a  heat-pipe  configuration.  Sealing  should  be  accomplished  with 
elastomeric  o-rings  which  are  protected  from  high  temperature  and  alkali 
attack  with  a  heat-pipe  configuration.  Metalized  and  brazed  seals, 
although  quite  competent,  are  much  more  expensive  and  difficult  to  use 
than  o-ring  seals. 

The  heat-pipe -oven  design  obviates  many  of  the  gas-handling  diffi¬ 
culties  inherent  in  the  flowing  system  design.  The  main  advantages  are 
elimination  of  aerosol  in  the  high-temperture  region,  low  bulk  veloci¬ 
ties,  and  continuous  purification  of  the  alkali  metal  vapor  and  reactant 
gas.  A  jpotential  difficulty  of  the  heat-pipe  approach  lies  in  the  gen¬ 
eration  of  aerosols  in  the  region  of  vapor  condensation.  These  aerosols 
are  problematic  only  if  photoionizing  or  diagnostic  light  must  pass 
through  them,  and  if  extinction  is  unacceptably  high. 

The  recombination  rate  is  most  easily  measured  in  the  high-collision- 
frequency,  low -plasma -frequency  conditions  of  the  afterglow  with  a 
lumped  impedance  conductivity  measurement.  Directly,  capacitively ,  or 
inductively  coupled  measurements  are  all  possible.  The  latter  two 
coupling  approaches  alleviate  concerns  of  plasma  perturbation  through 
surface  recombination  on  the  probe  or  photoemission,  and  of  operation  in 
the  alkali  environment  by  coupling  through  the  reaction  cell  wall. 
Optical  diagnostics  should  be  explored  as  a  check  on  the  recombination 
rate  measurement  and  as  a  possible  measure  of  the  electron  temperature. 
The  success  of  optical  diagnostic  techniques  will  depend  on  the  plasma 
emission  signal  level  and  on  the  extent  of  photo ionization,  quantities 
which  are  best  determined  experimentally. 

A  design  for  the  next -generation  experiment  which  draws  on  the  above 
features  would  employ  a  polycrystalline  alumina  2 "-diameter  tube.  A 
4-mm  diameter  Lucalox  flashlamp  like  the  high-temperature  version  we 
have  previously  built  would  be  mounted  coaxially.  A  split  wick,  over¬ 
pressured  heat-pipe -oven  design  would  be  used  having  cold  ends,  a  wick 
region  of  controlled  intermediate  temperature,  and  a  center  high- 
temperature  reaction  region.  A  conductivity  probe  could  be  capacitively 
coupled  through  the  ceramic  cell  wall.  Optical  diagnostic  access  would 
be  provided  through  fused  silica  end  windows  sealed  with  o-rings. 
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ae  Effects: 


Plasma-Electrode  Interaction 


2.2.1  Introduction 


The  operation  of  plasma  devices,  such  as  magnetohydrodynamic 
generators,  thermionic  converters,  and  magnetoplasmadynamic  thrusters, 
depends  on  their  ability  to  transfer  current  to  and  from  the  plasma  via 
electrodes,  yet  little  is  known  about  the  current  transfer  process  or  the 
electrical  boundary  layer  which  forms  in  the  near-electrode  region. 

Current  transfer  at  the  electrode  can  be  either  diffuse  or  constricted 
in  nature.  The  diffuse  mode  occurs  for  relatively  low  current  densities. 
In  this  mode,  the  boundary  layer  resistance  is  large,  and  a  substantial 
voltage  drop  occurs  across  the  layer.  The  resulting  power  loss  to  the 
plasma  decreases  the  overall  device  efficiency.  For  larger  current 
densities,  current  constrictions  form,  and  the  current  is  transferred  to 
the  electrode  through  one  or  more  highly  conductive  channels  or  arcs.  The 
high  temperatures  within  the  arc  rapidly  erode  the  electrode  surface, 
shortening  the  device  life. 

The  work  described  herein,  is  aimed  at  improving  the  present 
understanding  of  the  diffuse  current  transfer  mode,  for  the  case  of 
a  seeded  thermal  plasma  at  pressures  on  the  order  of  one  atmosphere, 
flowing  between  plane-parallel  electrodes.  In  addition,  the  dominant 
mechanisms  and  conditions  which  lead  to  the  transition  to  the  constricted 
current  transfer  mode  have  been  investigated. 


The  plasma-electrode  boundary  layer  consists  of  four  major 
overlapping  regions.  For  a  plasma  flowing  over  a  plane  electrode  (see 
Fig.  2.1)  ,  at  an  arbitrary  streamwise  (z)  location  one  finds 


i  )  The  Gasdynamic  Boundary  Layer,  of  scale 

c  ^oo  ~  Telectrode  i  n  —  2 

0  «  10  m, 

dT  {X)  /dx\^=o 


iV 


which  includes  both  the  thermal  and  velocity  boundary  layers. 
Here,  Too  is  the  core  plasma  temperature  at  the  edge  of 
the  boundary  layer  and  TcUctrode  is  the  plasma  temperature  at 
the  electrode  surface.  In  this  region,  the  plasma  is  in 
thermodynamic  equilibrium  at  the  local  temperature  T  (x)  and 
the  electric  field  E(x)  =Jja(x)  is  established  where  '  is 
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Fig.  2.1.  Coordinate  system  for  plasma  flow  over  a  plane 
electrode . 


the  current  density  and  a (x)  is  the  equilibrium  electrical 
conductivity; 


ii)  The  Ionization  Nonequilibrium  Layer,  of  scale 


tR  =  \J‘2.DalP(Te)  n|o 


Ion  Recombination 
Length 


w  10”^  m. 


Here,  P  (T^)  «  10"^°  m®/sec  is  the  three-body  recombi¬ 
nation  coefficient  divided  by  the  plasma  concentration  in 
the  freestream,  no,,.  T*  is  the  electron  temperature,  and 
Da  is  the  ambipolar  diffusion  coefficient.  In  this  region, 
the  generation  rate  of  electron-ion  pairs  exceeds  the 
recombination  rate  to  supply  net  fluxes  of  charged  particles; 


Hi)  The  Sheath,  of  scale 

In  this  layer,  charge  neutrality  is  not  preserved  and 
large  electric  fields  develop  which  must  be  calculated  from 
Poisson's  equation.  The  sheath  matches  the  fluxes  from  the 
ionization  nonequilibrium  layer  to  the  emission,  absorption, 
and  recombination  processes  at  the  electrode  surface;  and 


iv)  The  Knudsen  Layer,  of  scale 

ti  =  iCi/fi)  ~  4  •  (Ce/t-e) 


Mean  Free 
Path 


w  10 


m, 


where  free  molecule  flow  takes  over  from  a  continuum 
description.  While  the  mean  free  path  and  Debye  length  are  of 
the  same  order,  the  sheath  is  normally  several  Debye  lengths 
thick,  so  that  the  sheath  is,  for  the  most  part  collision 
dominated. 

As  indicated  above,  the  scale  lengths  are  ordered  according  to  the 
relation  4  ~  ^  ~  Ad  These  divisions  are  quite  arbitrary,  and 

no  distinct  boundaries  exist. 

2.2.2  Research  Objectives 

1.  To  develop  analytical  and  computer  models  and  to  predict  the  diffuse 
mode  current-voltage  characteristics  for  anodes  and  cathodes  under 
a  wide  range  of  operating  conditions. 

2.  To  use  the  computer  model  to  investigate  the  conditions  whereby 
electron  temperature  nonequilibrium  effects  become  dominant  and 
lead  to  the  transition  from  a  diffuse  to  a  constricted  mode  of 
current  transfer. 

2.2.3  Status  of  Research 

A  1-D,  continuum  model  has  been  developed  to  examine  the  DC  discharge 
characteristics  for  a  collision-dominated  weakly  ionized,  thermal  plasma 
contacting  a  plane  electrode. 

A  novel  feature  of  this  formulation  is  that  a  single  set  of  governing 
equations  is  solved  throughout  the  entire  boundary  layer  region,  hence 
a  solution  is  obtained  which  is  valid  through  both  the  ionization 
nonequilibrixim  layer  and  the  sheath.  Solutions  are  found  over  a  wide 
range  of  plasma  and  electrode  conditions;  our  goal  has  been  to  examine  the 
nature  of  the  plasma-electrode  interaction  instead  of  simply  modeling 
one  specific  experiment  or  set-up. 

The  current  approach  is  in  contrast  to  asymptotic  methods  where 
differing  approximations  are  made  for  each  of  the  overlapping  sublayers, 
separate  sets  of  equations  are  solved  for  each  sublayer,  and  the 
solutions  are  matched  at  the  boundaries.  The  latter  method  may  only 
be  used  in  the  limit  of  some  small  parameter  (the  sheath  thickness  in 


this  case),  and  it  is  complicated  by  the  difficulty  in  choosing  the 
(non-physical)  boundary  positions. 


The  charged  particle  behavior  is  governed  by  the  electron  and  ion 
conservation  equations  for  mass  and  momentum,  together  with  Poisson's 
ecjuation  and  the  electron  energy  equation.  The  dependent  variables 
include  the  electron  and  ion  n\amber  densities  (rie  (x)  and  n,’(x)),  the 
electron  and  ion  particle  fluxes  {Fgix)  and  Fix)),  the  electric  field 
{E(x)),  and  the  electron  temperature  (Te(x)).  Additionally,  the  electric 
field  may  be  integrated  to  obtain  the  scalar  potential  (^(i)).  The 
independent  variable  (x)  is  the  normal  distance  from  the  electrode 
surface. 

The  basic  equations  were  developed  subject  to  the  following 
assumptions : 

1)  One-dimensional,  steady-state  system  with  a  plane  electrode 
contacting  a  semi-infinite  plasma  in  the  half-space 
0  <  X  <  oo; 

2)  Ions  and  background  neutral  particles  are  at  the  same 
temperature  TneuUai  (*)  =  Ti  (x)  =  T  (x)  ; 

3)  Net  production  (or  loss)  of  charged  particles  fi(x),  is 
included  via  rate  expressions  for  electron  impact  ionization 
and  three-body  recombination  on  electrons; 

4)  Charged  particle  inertia  is  neglected  due  to  the  highly 
collisional  nature  of  the  plasma; 

5)  There  is  no  external  magnetic  field,  and  such  fields  due  to 
plasma  currents  are  neglected; 

6)  A  net  current  to  or  from  the  electrode  is  allowed,  and 
is  assumed  to  be  due  to  an  applied  potential  between  the 
electrode  and  the  distant  plasma;  and 

7)  The  electrode  surface  is  assumed  to  absorb  incident  electrons 
and  ions.  No  electron-  or  ion-emission  mechanism  is 
specifically  imposed,  but  such  emission  as  is  required  by  the 
solution  is  assumed  to  be  provided  by  thermionic  emission  of 
electrons  or  surface  ionization  of  positive  ions. 

To  keep  the  model  as  general  as  possible,  all  variables  are  written 
in  non-dimensional  form  (denoted  with  an  overbar) .  Then  the  governing 


e  *,«'#  n't 


equations  are  written: 


Mass  Conservation 


^  =  n=  {T^/%xp[e,  (1-^)]  -  neO^}  (2.1a) 


dr, 

dx  " 


Momentum  Conservation 


drij  ^  f ^  ^  A 

dx  T  [dx  J  1+fi  T 

dn,  ^  ^  / dT,  ^  £, 

dx  T,[dx  )  l  +  fiT, 


Poisson's  Equation 


dE  1  _  _ 

-jz:  =  (ni  -  n,) 

dx 


Electron  Energy  Conservation 


(fcT,)  +  ^  r,E 

b  dx  y  dx  J  I  +  n  dx  1  +  A* 


r^en,{T,-T)  -£<-^n 
1  +  /I  1  +  M 


with  the  two  ancillary  relations 


Current  Conservation 


Potential 


^ (I)  =  -  f  E(i)  df 

Jo 


(2.16) 


{2.2a) 


(2.26) 
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The  variables  are  non-dimensionalized  as  follows: 


X  =  xjtR 

Hi  (I)  =  m  (x)  /tir 
Fiix)  =  Fid)  /Tr 
E(x)  =  B(x)  /Er 
(X)  =T^(x) 


h{.x)  =  h(x)  /hR 

Tieix)  =ne(x)  /riR 
r«  (X)  =Fe(x)  /Fr 
^(X)  =  #  (Z)  J^R 

T(x)  =T(x)lTe^ 


J  =  J/Jr 


with  the  reference  quantities: 


(2.7) 


Ir  =  Ion  Recombination  Length 


—  y/^  ^a/ P  (^eji ) 

(2.8o) 

Zero  Current,  Far-Field,  Equilibrium 

riR  = 

Charged  Particle  Concentration 
s  tii  g  (X  — ►  oo,  J  =  0) 

(2.86) 

njt  =  0(TgJ 

(2.8c) 

Fr  =  2  DaTlR/iR 

(2.8d) 

Ten  =  Zero  Current,  Far-Field,  Electron  Temperature 

(2.8c) 

Er  =  kTen/elR 

(2.8/) 

^R=kTeJt 

(2.85) 

Tr  —  <^R  Er 

(2.8/1) 

crR  =  nRe  (fii  +  fig) 

(2.81) 

where  Z?a  «  {Di  iie  + De  fii)  /  Ifii+fie)  is  the  ambipolar  diffusion  coefficient  and 
^  (Tg)  =  1.09  X  10~^°Te  ^^^m®/sec  is  the  tnree-body  (electron)  recombination 
coefficient . 

Nine  free  parameters  must  be  specified:  1)  the  ion/electron  mobility 
ratio  ii,  =  Hi/fie  I,  a  function  of  the  type  of  ion  and  background  gas; 
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2)  e  =  Xjd/Ir  ^  1/  the  ratio  of  the  Debye  length  A/>  =  y/cok  Tgj, / hr  to 
the  ion  recombination  length;  3)  0  =  trie  i/en^/ [rrin  0  (Te^)  ratio  of 

the  elastic  to  ionizing  collision  frequencies  weighted  by  the  average 
factional  energy  loss  in  electron-neutral  collisions,  where  m„  is  the 
mass  of  the  neutral  particles  and  Uen  is  the  electron-neutral  collision 
frequency;  4) the  total  current  density  J,  where  J  <  0  for  cathodes 
and  J  >  0  for  anodes;  5)  the  non-dimensional  ionization  energy  , 
£,•  =  Ci/kTeni  where  €<  is  the  ionization  energy  for  the  ionizing  species 
and  k  is  Boltzmann's  constant;  6)  T  (x)  =  T  (x)  the  non-dimensional 

temperature  of  the  ion  and  background  neutral  particles,  which  is 
normalized  on  the  reference  temperature  and  7)  boundary  values  for 

the  quantities  Hi  (0) ,  fie  (0)  and  either  Te  (0)  or  dTe.ld'x\^. 


The  governing  equations  comprise  a  seventh- order,  coupled,  non-linear 
set  of  ordinary  differential  equations,  so  seven  boundary  conditions  are 
required.  Four  values  are  known  in  the  far  field,  *-+oo.  They  are  the 
particle  fluxes  Aco  A®'  electric  field  Ego,  and  the  electron 
temperature  .  Three  remaining  boundary  conditions  are  set  at  the 
electrode  surface,  5=0,  and  include  the  charged  particle  concentrations 
(^(0)  ,  ne(0)  )  and  either  the  electron  temperature  TgiO)  or  the  electron 
temperature  gradient  {dTe/dx\o)  . 


To  evaluate  the  far-field  conditions,  we  recognize  that  far  from  the 
electrode  surface  the  plasma  is  uniform  and  the  electrons  and  ions  are 
in  equilibrium.  The  charged  particle  concentrations  are  assumed  equal, 
n.  (I  — ►  oo)  =  (5  — »  oo)  =  rioo,  and  may  be  calculated  from  the  Saha  equation 


where  Tg^  •=  Tg  (x  — ♦  oo,  J) /Tg^, .  Equation  <2. 9)  is  necessary  because  ng,i(i) 
and  Tg  (5)  are  normalized  on  their  values  for  5  — »  oo  and  J  =  0,  that  is 


eoo  '*00 


1  for  J  =  0 .  Both  Ago  ”«j  raay  be  >  1  for  J  /  0  as  a  result  of 


Joule  heating  (J  •  E)  .  The  flux  and  electric  field  boundary  conditions 


are  then 


—  1  +  u  — 


(2.10O) 
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Eoo-J- 
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(2.106) 


(2.11) 


At  the  electrode  surface,  5=0,  for  a  species  which  is  perfectly 
absorbed  (no  reflection  or  emission) ,  the  usual  continuum  boundary 
condition  for  a  diffusion-type  equation  is  to  let  n,-,e(0)  =0.  As  is 
well  known,  this  boundary  condition  is  strictly  invalid  because,  for 
a  non-zero  flux  at  the  boundary  (7’(0)  =/  O),  the  mean  velocity  becomes 
infinite.  The  diffusion  equation  fails  within  a  distance  on  the  scale 
of  the  mean  free  path,  and  a  kinetic  description  is  necessary  to  treat 
the  Knudsen  layer  where  the  velocity  distribution  changes  from  a  nearly 
isotropic  form  to  a  purely  one-sided  form  as  the  boundary  is  approached. 
Such  an  analysis  [2.1]  shows  that  the  exact  solution  is  well  modeled 
by  the  diffusion  equation  with  the  boundary  condition  n(0)  =0  applied 
at  a  position  approximately  one  mean  free  path  outside  the  physical 
boundary.  Within  the  continuum  approximation,  the  mean  free  path  is 
implicitly  assumed  to  be  much  smaller  than  any  other  length  scale  of 
interest,  and  the  boundary  condition  is  normally  applied  right  at  the 
physical  boundary.  In  the  present  problem,  use  of  the  boundary  condition 
^*,<(0)  =  0  for  an  absorbing  electrode  implies  that  the  mean  free  path  is 
much  smaller  than  the  sheath  thickness.  This  requirement  is  satisfied 
for  many  plasmas  of  practical  interest. 

Letting  rigiO)  =0  does,  however,  lead  to  difficulty  with  regard  to 
the  electron  temperature  and  the  electron  energy  equation.  As  rig  —*  0, 
the  electron  thermal  conductivity  goes  to  zero  and  the  electron  energy 
equation  terms  for  thermal  conduction  and  collisional  energy  exchange 
with  the  background  neutral  particles  both  vanish,  n*  =  0  is  a  regular 
singular  point  for  this  equation,  since  the  thermal  conduction  terra 
contains  the  highest  derivative  in  the  electron  energy  equation.  The 
fundamental  problem  is  that  for  a  perfectly  absorbing  electrode  where 
ne(0)  =  0,  the  electron  velocity  is  one-sided  and  the  concept  of 
temperature  is  no  longer  valid.  This  problem  and  the  electron  temperature 
boundary  conditions  will  be  addressed  more  completely  when  we  discLiss 
the  solutions  which  include  the  electron  energy  equation  ,eq.  (2.4)  .  Use 
of  the  condition  rigiO)  =0  is  quite  acceptable  within  the  isothermal  limit 
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where  Tf(x)  is  assumed  to  be  known,  and  the  electron  energy  equation  may 
be  neglected. 

When  there  is  electron  emission  or  ion  emission  from  the  electrode, 
the  boundary  condition  generally  corresponds  to  fii,e(0)  >  0.  The  value  of 
fii,e(0)  can  be  related  to  the  emission  current  density  by  the  following 
approximate  argument,  which  is  discussed  on  a  dimensional  basis,  for  the 
case  of  electrons. 

The  electron  velocity  distribution  at  any  point  may  be  represented  as 
the  sum  of  an  outwardly-directed  half-Maxwellian  of  concentration  rigix) 
and  an  inwardly-directed  half-Maxwellian  of  concentration  n~  {x) ,  so  that 


Tie  (i)  =  nl  (I)  +  rig  (x) 


(2.12) 


The  outwardly  and  inwardly  directed  random  fluxes  are  then  n*  (x)  [ce  (x)  /2] 
and  n~  {x)  [c*  (x)  /2]  where  c,.  (z)  =  \/8  fcT*  (x)  /jr  me  is  the  mean  thermal  speed. 


Hence  the  net  flux  is 


Te  (I)  =  (fie  (x)  -  n~  (X)) 


Fromeqs.  (2.12-2.13)  we  find 


n.U)  -  +  — 


A  (X) 

e(l)/2) 


/  ,  A(x)  • 

(Ce(X)/2).- 


(2.13) 


(2.140) 


(2.14t) 


Using  the  definition  Feix)  =ne(x)Ve(x),  where  Ve  (x)  is  the  electron 
mean  velocity,  the  condition  that  n*  and  ri~  >  0  requires  that 
Ive  (x) I  $  Ce  (x) /2,  a  requirement  which  is  usually  met  for  conditions  of 
interest.  At  the  electrode  surface,  x  =  0, 


^  1  A(0)  1 

fie  (  0  )  =  ~  fie  (  0  )  +  — r 7 • 

'  2  L  '  (Ce(0)/2)  J 


The  positively  directed  flux  at  the  electrode  surface  is 
C(0)  =<(0)  ^4^  =  1  fle(0)^^  +  re(0) 
and  the  emitted  current  density  is 


e  r  ( 0 ) 

ne(0)  -^^  +  re(0) 


(2.15) 


(2.16) 


(2.17) 


Note  that  at  positions  other  than  i  =  0,  it  is  not  possible  to  identify  F* 
as  due  to  emitted  electrons  because  of  the  interchange  of  F*  and  F~  due 
to  collisions.  An  equivalent  expression  for  the  electron  emission  flux 
was  derived  by  Kalnavarns  and  Bienkowski  [2.2]  by  a  different  approximate 
argument . 

Physically,  J*  (0)  is  the  electron  emission  saturation  current  density 
given  by  the  Richardson-Dushman  law 


Jg  —  -^^eleeirode  /^Fglggtrode)  (2.18) 


where  is  the  surface  work  function,  TgUctrode  is  the  electrode 

temperature,  and  A  «  1.2  x  10®  A/m^K^.  Equation  (2.18)  is  applicable  when 
the  electric  field  at  the  electrode  is  zero  or  decelerating  (E(0)  >  0). 
When  there  is  an  accelerating  field  at  the  electrode  {E(0)  <  0),  field 
enhanced  emission  due  to  the  Schottky  effect  [2.3]  is  applicable; 


Jg - ^J'eUetrode  ~  |  ~  sj  ^^  ( 0  )  /4  JT  to  /A  rejeclrode  ^ 


JgO  65{p  ^0.44  s/  E  (0)  /Tgtggfrod^ 


(2.19) 

(2.190) 


where  j:,  is  the  emission  given  by  eq.  (2.18) •  In  practice,  extremely 
large  electric  fields  are  required  to  significantly  enhance  the  electron 
emission. 


C.  Solutions 

The  above  equations  and  boundary  conditions  completely  define  the 
behavior  of  the  electrons  and  ions  for  the  diffuse  discharge  with 
or  without  charged  particle  emission  by  the  electrode.  The  solutions 
to  these  equations  describe  the  spatial  profiles  for  the  charged 
particle  concentrations  and  fluxes,  the  electric  field,  and  the  electron 
temperature  as  a  function  of  the  nine  free  parameters.  This  is  a  notable 
cjuantity  of  information.  It  is  important  to  keep  in  mind  that  from 
a  practical  standpoint,  the  cur rent -voltage  characteristics  are  what 
is  usually  measured,  and  that  we  are  really  trying  to  determine  how 
the  current-voltage  characteristics  depend  on  the  plasma  and  electrode 
conditions  (the  free  parameters)  . 
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Our  approach  was  first  to  simplify  the  governing  equations  to  a 
level  for  which  an  analytical  solution  could  be  found.  The  simplifying 
assumptions  were  then  relaxed  in  three  steps,  and  the  effect  of  each 
assumption  on  the  solution  was  observed. 

The  problem  was  solved  in  the  following  four  limits: 

1)  Isothermal,  Quasi-Neutral  Plasma  (no  Sheath)  —  Poisson's  equation 
and  the  electron  energy  equation  are  neglected,  and  the  mass  ad 
momentum  conservation  equations  are  simplified. 

2)  Isothermal  Plasma  —  The  electron  energy  equation  is  neglected  but 
Poisson's  equation  is  retained  and  the  sheath  is  included  in  the 
solution . 

3)  Non-Isothermal  Plasma,  Simplified  Energy  Equation  —  The  electron 
energy  equation  is  included  in  a  simplified  form  which  does 
not  include  energy  transfer  due  to  conduction,  convection 
and  ionization/recombination.  The  background  gas  is  assumed 
isothermal,  but  the  electron  temperature  may  vary. 

4)  Non-Isothermal  Plasma,  Complete  Set  of  Ecpiations  —  The  complete 
system  of  equations  defined  by  eqs.  (2. 1-2. 4)  is  solved  for  an 
isothermal  background  gas. 

C.l  Isothermal.  Quasi-Neutral  Solution  (Analytic) 


If  we  consider  an  isothermal  plasma  where  the  electrons,  ions,  and 
neutral  particles  all  remain  at  the  same  tenperature,  the  electron 
energy  equation  is  redundant  and  may  be  neglected.  This  approximaticn 
is  valid  so  long  as  the  net  current  density  is  small  and  there  is  little 
Joule  heating  of  the  electrons,  or  the  collisional  coupling  between  the 
electron  and  the  background  neutral  particles  is  strong  (0  — »  oo)  . 

We  also  require  the  plasma  conditions  be  such  that  the  sheath  is 
infinitely  thin,  that  is,  e— >  0.  In  this  case,  ri,  (x)  =  rie  (x)  =  n(x) 
throughout,  Poisson's  equation  may  be  neglected,  and  the  solution  does 
not  include  the  sheath. 


The  mass  and  momentvim  conservation  equations  simplify  to 

dTi  _  _ 

— r  =  h  =  n  (1  -  n^) 
ax 

dTe 

dx  " 


and 


(2.20a) 

(2.206) 
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Eliminating  the  fluxes  between  (2.20)  and  (2.21)  yields 


(n£?) 

X 

d  (nE) 

X 


<fn 

dx^ 

d^n 

dx^ 


4m  r=. 


l+^t 


[n  (1  -  ri^)] 


4 

rni 


which  may  be  added  to  produce  a  single  nonlinear  equation  for  n: 


d^  n  _2 

— +  2n(l-n=')  =0, 
dx 


(2.210) 

(2.21b) 

(2.22a) 

(2.22b) 

(2.23) 


subject  to  the  boundary  conditions  n(0)  =  0  (no  particle  emission  at  the 
electrode)  and  n -+ 1  for  x  oo. 

The  solution  to  eq.  (2.23)  is 


n (I)  =  tanh (i) , 


(2.24) 


shown  in  Fig.  2.2.  The  net  generation  rate  is  given  by 

h(x)  =  sech^  (X)  tanh(x)  (2.25) 


which,  as  shown  in  Fig.  2.2  is  zero  for  i  =  0  and  x  — ^  oo.  The  center  of  the 
ionization  nonequilibrium  region  is  located  at  x  =  tanh“^ (3“^/^ )  ~  0.658, 
where  the  net  generation  rate  has  a  maximum  value  of  (2/3)  ^/2  ~  0.385.  We 
shall  find  later  that  the  magnitude  of  the  total  net  generation  rate  in 
the  boundary  layer  determines  the  saturation  current  for  both  anodes  and 
cathodes  in  the  absence  of  emission.  We  note,  for  reference,  that  both 
n  (x)  and  h  (i)  are  independent  of  (x  and  J,  and  that 


Total  Net  -n  -  1 

=  /  n  (X)  dx  =  -. 
Generation  Rate  Jq  2 


(2.26) 


The  other  dependent  variables  are  found  by  back-substituting  the 
solution  for  n  (x)  into  eqs.  (2.20-2.21)  .  The  electric  field  is  given  by 


E(x) 


-(m 


sech^ (x) 


coth (x) 


(2.27) 


which  may  be  separated  into  two  components:  an  ambipolar  field  which 
decreases  exponentially  with  x,  and  a  resistive  component  which  is 
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Fig.  2.2.  Profiles  of  normalized  charged  particle  concentrations 
and  net  generation  rate. 


necessary  to  drive  the  current  through  the  plasma.  Both  components  have 
a  simple  pole  at  *  =  0  due  to  the  coth  (x)  factor. 


The  ambipolar  field  is  always  negative,  i.e.,  directed  towards 
the  electrode.  The  resistive  field  is  in  the  same  direction,  for  a 
cathode  (J  <  0)  ,  but  in  the  opposite  direction  for  an  anode  (J  >  Q) . 
If  J  is  positive  and  large  enough,  the  field  at  the  anode  can  become 
positive.  Conventionally,  in  gas  discharge  terminology,  this  is  termed  a 
transition  from  a  negative  anode  fall  to  a  positive  one.  The  condition 
for  this  to  happen  is 


J>Ji  = 


1-^ 

1+/1 


»  1. 


(2.28) 


The  potential  is  given  by 

fix)  =  (- - ^  I  In  tanh  (i)  —  J  In  sinh(i) 

\  1  -"/i/ 


(2.29) 
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with  the  limiting  forms 


^(x)— ►  — ——J  ln(i)  for  X 1  (2.30o) 

Ll  +  A* 

and 

?  (X)  -+ -7  (x  -  In  2)  forx»l  (2.306) 

where  the  zero  of  potential  has  been  taken  at  x  — ►  00  for  J  =  0  because 
of  the  logarithmic  singularity  at  x=0.  Potential  profiles  for  both 
cathodes  and  anodes  are  shown  in  Fig.  2.3,  which  includes  cases  for  which 
the  anode  fall  is  negative  (0  >  J  >  J/)  and  positive  (J  >  J/  w  1)  .  The 
curves  are  asymptotic  to  the  broken  lines  ,eq.  (2.306),  whose  slopes 
represent  the  electric  field  E^o  =  J  in  the  distant  plasma. 

The  solutions  for  the  electron  and  ion  fluxes  are  also  composed  of  two 
terms 

Pi  (X)  =  --  sech^  (i) - — 

Te(x)  = sech^  (I)  +  .  (2.316) 

2  L  . 

The  first  term  represents  the  equal  ambipolar  fluxes  to  the  electrode; 
this  flux  is  independent  of  the  current  density  and  originates  from  the 
net  generation  of  electron-ion  pairs  in  the  ionization  nonequilibrium 
layer.  The  second  term  represents  the  fluxes  associated  with  the  current . 
The  saturation  behavior  of  the  electrode  current -voltage  characteristics 
may  be  explained  by  considering  eqs.  (2.31)  in  more  detail. 

For  a  cathode  {J  <  0)  ,  the  ion  flux  is  always  directed  towards  the 
cathode  and  increases  with  current  magnitude.  However,  the  electron 
flux  to  the  cathode  decreases  with  increasing  current  magnitude  until 
J  <  — 2/i/  {1  +  fi)  ,  at  which  point  there  must  be  a  net  electron  flux  leaving 
the  cathode.  This  means  that  the  cathode  must  emit  electrons  for  currents 
greater  than  the  relatively  small  value  Jcr  given  by 

_  -2fx 

=  (2.32) 

1  +M 

Jc  is  interpreted  as  the  cathode  saturation  current  in  the  absence  of 
electron  emission  by  the  cathode. 

For  J  =  Jc,  the  current  at  the  cathode  surface  is  carried  entirely  by 
the  arriving  ion  flux,  which  is  only  a  factor  (l+/x)  times  the  equal 


Profiles  of  normalized  potential  for  several  values 
of  normalized  current  density:  J  <  0,  cathode; 
J  >  0,  anode.  Broken  lines  indicate  the  asymptotes 
corresponding  to  the  uniform  electric  field  in  the 
distant  plasma. 


electron  and  ion  fluxes  arriving  at  the  electrode  for  zero  current. 
In  previous  analyses  of  the  cathode  sheath  in  the  absence  of  electron 
emission  [2.4,5],  it  has  often  been  assumed  that  the  maximum  ion  flux 
which  can  be  drawn  from  the  plasma  to  the  electrode  is  the  random  one-way 
flux,  nooC$/4.  The  results  derived  here  indicate  otherwise. 

For  a  given  current  density  J ,  the  electron  flux  in  the  distant  plasma 
must  always  be  larger  that  the  ion  flux  by  the  ratio  (l/fi)  . 


(2.33) 


Thus  for  a  cathode  (J  <  0)  there  is  an  ion  flux  drawn  out  of  the  distant 
plasma  toward  the  cathode,  and  a  much  larger  electron  flux  drawn  away 
from  the  cathode  into  the  plasma.  By  mass  conservation. 


=  r^iO) 


+  lidx. 
Jo 


(2.34) 


At  the  critical  current  density  Jct  in  the  absence  of  electron  emission, 
fe  (0)  =0,  so  that 


ndx. 

Jo 


(2.35) 


Physically,  this  means  that  for  J  =  Jc  the  total  net  generation  rate  of 
electrons  in  the  ionization  nonequilibrium  layer  is  used  to  supply  the 
electron  flux  in  the  distant  plasma.  As  discussed  previously,  the  total 
net  generation  rate  is  given  by 


Total  Net  1 

=  /  h(x)  dx  =  — 
Generation  Rate  Jo  2 


(2.26) 


and  is  independent  of  J .  Thus  the  maximiim  electron  flux  in  the  distant 
plasma,  in  the  absence  of  electron  emission,  is  = l/2,  corresponding 
to  a  maximum  current  of  (-J)  =  2p/  (1  +  /i)  . 

It  is  clear  that  the  current  limitation  is  caused  not  by  a  limitation 
on  the  flux  of  ions  which  may  be  drawn  to  the  cathode  from  the  plasma,  but 
by  a  limitation  on  the  total  net  generation  rate  of  electrons-ion  pairs 
in  the  ionization  nonequilibrium  layer  which  supply  the  electron  flux  to 
the  distant  plasma. 
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In  the  case  of  an  anode,  a  similar  current  limitation  phenomenon 
arises.  For  an  anode  (J  >  0)  the  electron  flux  is  always  directed  towards 
the  electrode  and  increases  with  increasing  current  density.  The  ion 
flux  decreases  with  increasing  current  until  for  J  >2j  (1  +  ,  a  current 

approximately  twice  J/  for  the  transition  to  a  positive  anode  fall,  there 
is  a  net  flux  of  ions  leaving  the  anode.  This  implies  that  the  anode  must 
emit  ions  for  currents  greater  than 


Ja  = 


(2.36) 


For  J  =  Jj^  the  current  at  the  anode  surface  is  carried  entirely  by  the 
arriving  electron  flux  which  is  a  factor  (1  + /i) //i  ~  l//i  ^  1  times  the 
arriving  electron  flux  in  the  aibsence  of  current. 

The  physical  reason  for  this  limitation  parallels  that  for  the 
current  limitation  in  the  case  of  a  cathode.  In  the  distant  plasma,  a 
given  current  implies  an  ion  flux  =  (l+/i)  J/4  directed  away  from 
the  anode.  In  the  absence  of  ion  emission  by  the  anode,  the  flux  in 
the  distant  plasma  must  originate  from  ions  generated  in  the  ionization 
nonequilibri\am  layer.  Hence,  the  maximum  ion  flux  in  the  distant  plasma 
is  limited  by  the  total  net  generation  rate  of  electron-ion  pairs  in  this 
layer,  and  the  maximum  current  which  corresponds  to  this  maximum  ion  flux 
is  given  by  J  a  . 

Thus,  the  anode  current  limitation  is  not  set  by  a  limitation  on  the 
flux  of  electrons  that  can  be  drawn  to  the  anode  from  the  plasma,  but  by  a 
limitation  on  the  ion  current  in  the  distant  plasma  set  by  the  total  net 
generation  rate  of  ions  in  the  ionization  nonequilibrium  layer.  Cheng 
et  al.  [2.6]  recognize  this  current  limitation,  and  derive  an  equivalent 
expression  for  the  value  of  saturation  current. 

To  summarize,  in  the  absence  of  charged  particle  emission  by  the 
electrode,  current  saturation  occurs  due  to  a  limitation  on  the  flux  of 
charged  particles  drawn  away  from  the  electrode  in  the  distant  plasma, 
set  by  the  finite  supply  of  particles  generated  in  the  ionization 
nonequilibrium  region. 

In  general,  a  fraction  of  the  repelled  species  (electrons  for 
cathodes,  ions  for  anodes)  which  are  generated  in  the  ionization 
nonequilibrium  layer  are  drawn  to  the  electrode,  with  the  remaining 
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particles  drawn  into  the  plasma.  The  saturation  currents  Jc  ~  +  ft) 

and  Jji  =  2/  (1  +/i)  are  the  limiting  cases  whereby  the  total  net  generation 
rate  in  the  ionization  nonequilibrium  layer  is  used  to  supply  the 
repelled  particle  flux  in  the  distant  plasma  (directed  away  from  the 
electrode) .  A  third  notable  case  occurs  for  the  absence  of  current, 
J  =  0,  where  the  total  net  generation  rate  supplies  the  equal  ion  and 
electron  fluxes  to  the  electrode,  since  there  is  no  net  flux  in  the 
distant  plasma.  We  have  also  found  that  the  sign  of  the  electric  field  at 
the  surface  of  the  electrode  changes  for  J  =  J/  =  {1  —  n)  /  (1  +  fi)  ~  1,  which 
is  approximately  half  the  value  of  the  anode  saturation  current. 

Note  that  these  critical  currents  Jc,  Ja  and  J/  depend  only  on  the 
mobility  ratio  /i.  In  the  limit  of  /i  — ►  1,  the  anode  and  cathode  saturation 
behavior  becomes  symmetric  in  J  (i.e.,  —*  |Jd  — +  1),  and  the  surface 

electric  field  changes  sign  for  J  =  0 . 

While  much  has  been  learned  from  the  quasi-neutral  solution,  it  is 
unable  to  describe  the  sheath,  and  hence  cannot  be  used  to  calculate 
current-voltage  characteristics,  for  which  solutions  which  include 
Poisson's  equation,  are  required. 

C.  Isothermal  Solution  (Numerical) 

In  this  part,  we  consider  e/0.  Poisson's  equation  is  included  in 
the  set  of  governing  equations,  and  the  solution  includes  the  sheath. 
We  still  assume  an  isothermal  plasma  for  which  the  electrons,  ions  and 
neutral  particles  all  remain  at  the  same  temperature  (Te(x)  =  T  (i)  =  1)  , 
and  we -neglect  the  electron  energy  equation. 

Letting  (T*  (r)  =  T (x)  =  1)  simplifies  the  continuity  and  momentum 
equations  significantly  from  eqs.  (2.1-2. 2) .  The  governing  equations 
reduce  to 


Mass  Conservation 


i£i 

dx 

dT, 

dx 


=  h  =  Tie  ( 1  -  Me  n,  ) 

=  n 
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(2,37a) 


i2.31b)  \ 


Momentum  Conservation 


dUi 

dx 

dUe 

dx 


--n,E-^T, 

l+fi 


(2.380) 

(2.386) 


and 


Poisson's  Equation 


dx 


1  _  _ 

=  —  (ni  -  rie) 


(2.3) 


along  with  the  five  free  parameters  n,-  (0) ,  rie  (0) ,  J,  e  =  Xo/f-R  and  fi  -  \iil\ie  • 

These  equations  were  solved  niomerically  by  the  method  of 
quasilinearization,  combined  with  orthonormalization  and  superposition. 
Use  of  a  numerical  method  required  that  the  far-field  boundary  (i  — ♦  oo)  be 
located  at  a  finite  value  of  x.  We  chose  a  position  located  sufficiently 
far  into  the  plasma  (at  x-Xf,  say  0.658)  to  ensure  quasi-neutrality, 
and  well  past  the  ionization  nonequilibrium  layer,  i.e.,  where  n  1 . 
Equations  (2.27,2.31)  are  then  used  with  i  =  ij,  (<^  typically  5)  for  the 
boundary  values  of  the  dependent  variables  Eco^I'ea,  Aa,  • 

A  solution  database  has  been  generated  for  a  wide  range  of  the 
parameters  J ,  e,  and  n,  and  for  several  values  of  ion  and  electron 
emission  by  the  electrode.  Our  aim  has  been  to  observe  how  charged 
particle  emission  by  the  electrode  and  the  presence  of  uhe  sheath  modify 
the  conclusions  drawn  from  the  quasi-neutral  solution. 

The  results  have  been  organized  into  three  major  areas:  (i)  an 
electrically  isolated  (floating)  electrode  (J  =  0)  with  and  without 
electron  emission;  (ii)  a  cathode  (J  <  0)  with  and  without  electron 
emission;  and  (iii)  an  anode  (J  >  0)  with  and  without  ion  emission. 

We  shall  examine  only  the  results  for  the  critical  currents  discussed 
in  the  previous  section.  For  illustrative  purposes,  we  choose  to 
discuss  solutions  for  c=10"^.  This  permits  both  the  sheath  and  the 
ionization  nonequilibrium  layer  to  be  displayed  together  on  a  linear 
scale.  Atmospheric  pressure  plasmas  of  practical  interest  tend  to  have 


much  smaller  values  of  c,  so  we  have  computed  solutions  in  the  range 
10"^  <  c  <  10“^.  Comments  made  about  the  solutions  for  e-  10“^  apply 
equally  well  for  smaller  values  of  e,  unless  stated  otherwise.  The 
solutions  are  also  calculated  for  /i=  1/300,  but  all  of  the  dependent 
variables  are  quite  insensitive  to  the  value  of  fi. 


This  is  the  familiar  case  of  a  non-emitting,  perfectly  absorbing, 
insulating  surface  or  floating  electrode  contacting  the  plasma.  Profiles 
of  He,  tit  and  H  are  compared  with  their  quasi-neutral  values  nqti  and  Aqif 
in  Fig.  2.4  (a) ,  with  the  corresponding  profiles  of  /’«,  A,  E  and  <P,  shown 
in  Figs .  2.4  (b)  and  (c)  . 

All  of  the  electrons  and  ions  generated  in  the  ionization 
nonequilibri\im  layer  flow  to  the  surface  and  constitute  equal  (negative) 
fluxes  Fg  and  F,  at  every  point.  These  fluxes  increase  with  decreasing  i. 
Both  He  and  n,  are  significantly  lower  than  Hqu  except  that  Hi  exceeds 
Hqn  in  a  very  thin  diffusion  layer  very  close  to  the  surface.  Throughout 
the  electrical  boundary  layer  n,-  >  He,  so  there  is  a  positive  ion  sheath 
which  merges  smoothly  with  the  ionization  nonequilibrium  layer.  Because 
the  profiles  of  n*  and  Hi  differ  from  that  of  Hqs,  the  peak  of  the  net 
generation  rate  profile  h  is  significantly  shifted  away  from  the  surface, 
but  retains  the  same  general  shape.  In  fact,  the  total  net  generation 
rate  is  only  very  slightly  changed  from  the  quasi-neutral  value  of  1/2. 

If  we  (arbitrarily)  define  the  outer  edge  of  the  ion  sheath  ij,  as  the 
point  where 


A  = 


In,-  -  ne\ 
(Hi  +  fie)  /2 


0.1, 


(2.39) 


then  we  find  Zj,  =  0.85.  Hence,  since  c-  0.1,  the  sheath  thickness  in  terms 
of  Debye  lengths  (of  the  distant  plasma)  is  8.5.  It  is  found  that  as  c 
decreases,  the  sheath  thickness  increases  as  measured  in  Debye  lengths, 
reaching  a  value  of  36.3  for  e=10~^.  This  result  justifies  the  use  of 
continuiim  theory,  even  when  the  mean  free  path  for  electrons  and  ions  is 
comparable  to  or  greater  than  the  Debye  length,  because  the  sheath  is 
collision  dominated. 
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The  potential  decreases  toward  the  electrode  and  the  total  potential 
drop  is  —7.33  {kToo/e)  for  c-  10“^.  The  magnitude  of  this  potential  drop 
increases  slowly  with  decreasing  e,  reaching  a  value  of  —10.6  (fcToo/e)  for 
6=  10“^. 

C.2.2  Electron-Emitting.  Floating  IJ = 0)  Electrode 

When  the  above  conditions  are  modified  to  allow  electron  emission 
(ne(0)  >0),  an  electron  sheath  forms  immediately  against  the  electrode 
surface,  and  merges  smoothly  into  the  ion  sheath  further  out.  The  ion 
sheath  is  weakened  as  compared  with  the  non-emitting  case,  that  is  |ni  —  ne| 
decreases  everywhere,  and  the  potential  fall  decreases  in  magnitude. 

With  increasing  electron  emission,  the  width  and  strength  of  the 
electron  sheath  increases  while  the  strength  and  width  of  the  ion  sheath 
decreases.  As  a  result,  the  total  potential  drop  decreases  dramatically; 
the  electrode  floats  up  towards  the  plasma  potential,  becoming  equal  to 
it  for  n(0)  =1.  This  behavior  is  consistent  with  the  use  by  Langmuir  [2.7] 
and  others  [2.2, 8, 9]  of  a  thermionically  emitting  probe  to  measure  the 
plasma  potential.  For  a  given  value  of  e,  there  exists  a  critical  value 
of  electron  emission,  measured  by  «««,„  (0)  f  equal  to  0.223  for  e=10“^, 
at  which  the  surface  field  changes  sign  from  negative  to  positive  and  a 
potential  minimum  occurs  at  the  electrode  surface.  As  ne(0)  is  increased 
above  the  critical  value,  the  potential  minimum  moves  out  from  the 
electrode  surface  and  the  depth  of  the  potential  well  increases.  There 
is  also  a  minimum  in  Ue  associated  with  the  potential  minimum,  due  to 
the  fact  that  electrons  approximately  satisfy  the  Boltzmann  relation 
n  oc  exp^.  It  should  also  be  noted  that  although  the  electrode  potential 
floats  up  to  the  plasma  potential  for  ne(0)  =  1,  the  plasma-sheath  is 
not,  by  any  means,  a  region  of  uniform  potential  and  charged  particle 
concentration,  as  is  clear  from  Fig.  2.5. 

To  our  knowledge,  this  work  represents  the  first  satisfactory 
quantitative  analysis  of  an  emitting,  floating  electrode,  apart  from  the 
very  approximate  treatment  originally  done  by  Langmuir.  Figure  2.6  shows 
the  total  potential  drop  as  a  function  of  electron  emission  for  various 
values  of  e.  It  indicates  that  except  for  very  small  values  of  emission, 
the  potential  drop  is  independent  of  c,  and  that  while  the  electrode 
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potential  approaches  the  plasma  potential  for  relatively  small  values  of 
He  (0) ,  it  only  tends  to  exact  equality  as  ne(0)  tends  to  unity. 


Fig.  2.6.  Total  potential  drop  A^total  as  a  function  of  ne{0),  for 
various  t.  Parameter  values  ft  =  1/300,  J  =  0,  rij  (0)  *  0  . 


One  other,  unexpected,  result  should  be  mentioned.  With  increasing 
electron  emission,  there  is  a  significant  shift  of  the  region  of  net 
generation  tovards  the  electrode  surface,  because  of  the  enhanced  value 
of  He  in  the  electron  sheath.  However,  despite  the  dramatic  change  in 
the  profile  of  h,  the  total  net  generation  rate  remains  essentially 
invariant  at  the  value  1/2  given  by  the  quasi-neutral  approximation. 


C  ■  2  ■  3  Non-Emitting  Cathode 

We  now  consider  what  happens  when  a  non-emitting  electrode  draws 
a  net  current,  or  equivalently,  is  biased  more  negative  with  respect 
to  the  plasma  than  its  floating  potential  for  J-0.  Due  to  resistive 
losses,  the  plasma  is  no  longer  at  a  uniform  potential,  so  we  define 
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the  total  potential  drop  in  the  layer,  A^totait  as  the  integral  of  the 
space-charge-induced  electric  field.  In  this  way,  the  current-dependent 
electric  field  needed  to  overcome  the  finite  plasma  conductivity  is 


neglected.  Hence, 


/■**  —  — 

-  /  [EiX)  -  Eoo]  dx 
Jo 

—  +  Xb  Eqo 


(2.40) 


(2.40a) 
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where  the  zero  of  potential  is  ta)cen  at  the  electrode  surface  (#(0)  -  0)  . 
This  is  shown  graphically  in  Fig.  2.7. 


Fig.  2.7.  Total  potential  drop  in  the  boundary  layer. 


As  the  current  is  increased,  the  ion  sheath  becomes  wider  and  stronger 
because  the  electron  concentration  has  a  very  small  value  for  greater 
distances  from  the  surface.  At  the  same  time,  the  (negative)  potential 
fall  and  sheath  field  increase  and  tend  to  very  large  values  as  the 
saturation  current  Jc  is  approached.  As  was  noted  earlier,  in  practical 
terms,  the  saturation  current  density  for  non-emitting  cathodes  is  quite 
small . 

Figure  2 . 8  shows  the  solution  near  saturation,  for  J  «  Jc  “  “2/i/  (1  +/i)  . 
A  thic)?,  strong  ion  sheath  drives  the  ionization  nonequilibrium  region 
away  from  the  electrode  surface  and  into  the  plasma.  There  is  negligible 


net  generation  within  the  strong  part  of  the  sheath,  because  of  the  lack 
of  ionizing  electrons.  The  bulk  of  the  net  generation  occurs  in  the  weak 
part  of  the  sheath  where  it  merges  into  the  ionization  nonequilibrium 
layer.  This  tail  would  comprise  the  overlap  region  for  an  analysis 
matching  inner  and  outer  expansions,  and  illustrates  the  difficulty  in 
trying  to  match  approximate  solutions. 

All  of  the  dependent  variables  vary  smoothly  through  the  sheath,  the 
ionization  nonequilibrium  region,  and  into  the  uniform  distant  plasma. 
The  quasi-neutral  solution  (nqif  =  tanh  (i)  )  is  drawn  for  reference  in  the 
top  figure,  to  show  the  significant  effect  that  the  sheath  has  on  number 
density  profiles  throughout  the  entire  boundary  layer. 

The  electron  flux  at  the  cathode  surface  {x=  0)  is  essentially  zero; 
the  presence  of  the  sheath  only  perturbs  the  number  density  profiles, 
it  does  not  affect  the  total  net  generation  rate  in  the  layer.  This 
conclusion  is  valid,  independent  of  the  value  of  e. 

C.2.4  Electron-Emitting  Cathode 

When  electron  emission  from  the  cathode  is  allowed,  ne(0)  >  0,  it 
considerably  changes  the  appearance  and  behavior  of  the  sheath  as 
compared  to  a  non-emitting  cathode.  An  electron  sheath  is  formed 
immediately  in  front  of  the  cathode,  and  merges  smoothly  into  the  ion 
sheath.  The  emitted  electrons  weaken  (partially  neutralize)  the  ion 
sheath,  and  as  a  result,  the  magnitude  of  the  potential  fall  decreases 
for  a  given  current. 


As  ne(0)  is  increased  for  a  given  current  J,  the  electron  sheath 
becomes  stronger  and  thicker,  and  one  finds  a  critical  value  (i  =  0,  J) 
for  which  the  electric  field  at  the  cathode  surface  changes  sign  from 
negative  to  positive,  and  a  potential  minimum  occurs  at  the  electrode 
surface.  As  fie  (0)  is  increased  further,  the  potential  minimum  moves  out 
from  the  cathode  and  into  the  plasma. 

The  condition  whereby  E[0)  >  0  is  known  as  (electron)  space-charge- 
limited  emission,  to  distinguish  it  from  emission-  (or  temperature-) 
limited  emission  when  E{0)  <  0  and  there  is  no  potential  minimum.  In  the 
space-charge-limited  case,  the  full  emission  capability  of  the  cathode  is 
not  realized.  The  slower  electrons  in  the  half-Maxwellian  distribution 


of  emitted  electrons  are  reflected  by  the  retarding  electric  field  and 
are  returned  to  the  cathode.  Only  the  electrons  with  sufficient  energy 
to  overcome  the  potential  "hill"  presented  by  the  potential  minimum  are 
able  to  contribute  to  the  current  leaving  the  cathode. 

For  a  given  value  of  n<!(0)  <  (0,  J)  ,  the  cathode  is  emission 
limited.  If  ne(0)  is  held  fixed  and  the  current  magnitude  is  increased 
(made  more  negative) ,  the  system  remains  emission  limited  and  the 
(negative)  potential  fall  increases  rapidly. 

When  ne(0)  the  cathode  is  space-charge  limited.  As  the 
current  is  increased,  the  electron  sheath  thins,  the  potential  fall 
slowly  increases,  and  the  potential  minimum  becomes  shallower  and  moves 
baclt  toward  the  cathode  surface.  Eventually  a  value  of  J  is  reached  where 
the  cathode  becomes  emission  limited,  and  the  potential  fall  increases 
much  more  rapidly. 

This  behavior  is  summarized  in  Fig.  2.9,  where  the  cathode  potential 
fall  is  plotted  versus  (0)  with  J  as  .  parameter .  The  contour  separating 
the  space-charge-limited  and  emission-limited  operation,  where  E(0)  =  0 
is  shown  as  a  dashed  line. 

C.2.5  Non-Emitting  Anode 

Starting  with  a  floating,  non-emitting  electrode,  we  let  the  plasma 
draw  a  net  current  from  the  anode,  or  conversely,  the  electrode  is  biased 
more  positive  relative  to  the  floating  potential.  As  the  current  is 
increased,  the  flux  of  collected  electrons  rapidly  increases,  while  the 
flux  of  collected  ions  slowly  decreases;  the  ion  sheath  becomes  thinner 
and  weaker  with  increasirig  current;  and  the  negative  anode  fall  decreases 
i  n  magnitude . 

The  electron  flux  increase  is  large  compared  with  the  ion  flux 
decrease,  and  most  of  the  current,  is  carried  by  electrons  which  are  drawn 
from  the  plasma.  This  contrasts  with  the  behavior  of  the  cathode,  where 
the  ions  were  the  dominant  cha.ge  carrier  because  the  electron  flux 
decreased  with  increasing  current,  leaving  only  the  ions  to  carry  the 
current . 

As  J  approaches  Jj  =  (1  +  /i)  /(I  —  /i)  ,  the  ion  sheath  almost  completely 
dissapears,  an  electron  sheath  forms  at  the  anode  surface,  and  the 
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Fig.  ?.9.  Cathode  potential  fall  A^total  versus  ne(0)  with  current 
7  as  a  parameter,  showing  the  distinction  between  the 
space-charge-limited  and  emission-limited  regimes. 


electric  field  at  the  anode  surface  changes  sign  from  negative  to 
positive.  This  is  shown  in  Fig.  2.10. 


While  a  floating  electrode  (J  -  0)  differentiates  anodes  from 
cathodes,  one  can  argue  that  the  current  J  “  Ji  for  field  inversion 
is  the  natural  division  determining  the  plasma-electrode  behavior. 
This  asymmetry  is  strictly  a  result  of  the  unequal  electron  and  ion 
mobilities.  As  was  stated  earlier,  in  the  limit  ft  —*  1,  Jj  0,  and  hence 
the  field  inversion  and  floating  electrode  conditions  merge. 


For  J  >  Jj  the  anode  electron  sheath  becomes  thicker  and  stronger 
and  the  positive  potential  fall  increases.  Figure  2.11  shows  that  near 
saturation  (for  JmJy^  =  2/(l■\■ft)),  the  anode  boundary  layer  is  composed 
of  an  electron  sheath  adjacent  to  the  electrode  which  merges  into  the 
ionization  nonequilibrium  region,  followed  by  the  uniform  plasma.  As 
before,  the  presence  of  the  sheath  does  not  alter  the  value  of  the 
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saturation  current  found  from  the  quasi-neutral  solution. 


Comparing  Figs.  2.8  and  2.11,  at  saturation  the  anode  sheath  is 
clearly  much  thinner  than  the  cathode  sheath,  and  the  net  generation 
region  is  very  nearly  identical  in  size  and  position  to  that  of 
the  quasi-neutral  solution.  As  the  current  is  increased  for  a 
non-emitting  anode,  the  charged  particle  profiles  remain  similar  to 
their  quasi-neutral  values.  The  opposite  is  true  for  cathodes;  as  the 
current  increases  to  a  cathode,  the  charged  particle  concentration 
profiles  deviate  strongly  from  those  of  the  quasi-neutral  solution. 

C .  2 . 6  lon-Emittino  Anode 

The  possibility  of  ion  emission  from  an  anode  is  relevant  to  the  cases 
of  the  emitter  of  a  thermionic  converter,  and  the  anode  of  a  noble  gas  MHD 
generator,  seeded  with  an  alkali  metal  and  having  hot  electrodes  made  of 
a  refractory  metal .  Basically,  ion  emission  is  significant  only  when  the 
work  function  of  the  electrode  is  higher  than  the  ionization  of  the  seed, 
and  significant  surface  ionization  can  occur.  Moreover,  ion  emission 
from  an  anode  only  makes  a  significant  difference  for  J  >  Jj,  when  there 
is  an  electron  sheath  and  the  ion  emission  (partially)  neutralizes  it, 
reducing  the  positive  potential  fall  and  allowing  more  electron  current 
to  be  drawn  from  the  plasma. 

The  effect  is  analogous  to  the  effect  of  electron  emission  from 
a  cathode  which  partially  neutralizes  the  ion  sheath,  reducing  the 
potential  fall  and  allowing  more  current  to  flow. 

To  summarize,  Poisson's  equation  has  removed  the  singularity  for 
the  electric  field  and  potential  at  the  electrode  surface  which  was 
found  with  the  quasi-neutral  solution.  As  a  result,  we  have  been  able 
to  calculate  the  current-voltage  characteristics  for  both  cathodes  and 
anodes  with  and  without  charged  particle  emission  by  the  electrode.  We 
have  also  found  that  tor  non-emitting  electrodes,  the  sheath  does  not 
alter  the  saturation  behavior  predicted  by  the  quasi-neutral  solution, 
so  long  as  the  electrons  remain  at  the  same  temperature  as  the  background 
neutral  particles . 
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Fig.  2.11.  Profiles  of  the  dependent  variables  for  a  non- 
emitting  anode  close  to  saturation.  Parameter  values 
f- 10-^7- 2,  n*(0)  -n,(0)  -0. 
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C.3  Non-Isothermal  Solution.  Simplified  Energy  Equation  (Numerical) 

Results  from  the  quasi-neutral  and  isothermal  models  show  that  the 
ionization  and  recombination  reactions  play  an  important  role  in  the 
current  saturation  behavior. 

A  necessary  condition  for  true  current  saturation  is  that  the  total 
net  generation  rate  in  the  boundary  layer  must  be  independent  of  current . 
This  condition  is  satisfied  for  the  quasi-neutral  solution  and  for 
non-emitting  electrodes  with  the  isothermal  model,  primarily  because  we 
have  assumed  that  the  electron  temperature  remains  constant.  However, 
for  many  plasma  conditions.  Joule  dissipation  {J  •  E)  may  heat  (or 
cool)  the  electrons,  whence  the  electron  temperature  profile  becomes 
a  function  of  the  current  density.  The  ionization  and  recombination 
reaction  rates  vary  with  the  electron  temperature  such  that  heating  the 
electrons  also  increases  the  total  net  generation  rate.  The  end  result 
is  that  by  increasing  the  electron  temperature  and  total  net  generation 
rate.  Joule  heating  can  change  the  value  of  the  saturation  currents,  or 
may  remove  the  saturation  phenomenon  entirely. 

To  investigate  this  effect,  the  electron  energy  equation  is  added  to 
the  model  so  that  the  electron  temperature  may  differ  from  the  ion  and 
neutral  particle  temperatures,  which  we  still  assume  to  be  equal.  For 
the  moment,  we  consider  a  simplified  form  of  the  electron  energy  equation 
given  by 


-3en^  (Te-T) ,  (2.41) 

which  balances  the  rate  of  energy  exchange  between  the  electrons  and 
the  electric  field  (Joule  heating/cooling),  and  the  rate  of  collisional 
energy  exchange  between  the  electrons  and  the  background  neutral 
particles.  The  free  parameter  6  is  a  measure  of  the  collisional  coupling 
between  the  electrons  and  the  background  neutral  particles.  For  O  -*  oo, 
the  electrons  become  perfectly  coupled  to  the  neutral  particles,  the 
electron  temperature  and  neutral  temperature  become  equal,  and  the 
solutions  go  over  to  the  isothermal  solution  of  the  previous  section. 

An  equivalent  expression  for  eq.  (2.41)  has  been  used  in  previous 
analyses  [2.10,11]  and  we  include  it  in  eq.  (2.41)  to  form  a  basis  for 
comparison  with  the  solutions  for  the  full  electron  energy  equation. 


Neglected  here,  are  the  electron  energy  equation  terms  which  account 
for  thermal  conduction,  energy  loss  by  the  electrons  due  to  the  net 
generation  rate  of  electron-ion  pairs,  and  energy  convected  by  the 
electron  flux.  These  terms  will  be  considered  in  the  next  section. 

The  electron  temperature  is  now  a  dependent  variable  of  the 
problem,  hence  the  temperature-dependent  forms  of  the  mass  and  momentum 
conservation  equations,  eqs.  (2. land!. 2) ,  are  solved  with  Poisson's 
equation,  eq.  (2.3),  and  the  electron  energy  equation,  eq.  (2.41).  In 
review : 
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Poisson's  Equation 

dE  1  _  _ 

=  *7  (ni  -  Me)  (2.3) 

dx  e' 

Using  an  algebraic  form  for  the  electron  energy  equation  does  not 
change  the  order  of  the  system  of  equations  from  that  of  the  isothermal 
problem.  Hence,  five  boundary  conditions  are  required;  and  Eoo 

in  the  far-field  and  ne^<(0)  at  the  electrode  surface.  However,  solving 
eq.  (2.41)  for  the  electron  temperature. 


Te  =  T  - 


FeE 

3  0ne 


(2.42) 


shows  that  T*  has  a  singularity  for  n*  — »  0 .  This  is,  in  part,  due  to 
the  fact  that  Joule  dissipation  and  electron  collisions  with  neutral 
particles  are  no  longer  the  dominant  electron  energy  exchange  mechanisms 


when  rte  — ►  0.  Also,  ng  =  0  is  only  a  reasonaOsle  assumption  for  distances  on 
the  order  of  one  mean  free  path  from  the  surface  of  a  perfectly  absorbing 
electrode.  In  this  region,  the  velocity  distribution  for  the  electrons 
becomes  completely  one-sided  and  the  translational  electron  temperature 
becomes  undefined. 

To  avoid  this  difficulty,  we  examine  solutions  for  0  yne(O)  1  (i.e., 
for  very  low  levels  of  electron  emission)  and  submit  that  the  solutions 
convey  the  same  practical  results  as  one  would  find  for  a  perfectly 
absorbing  electrode  (ne(0)  =  0)  .  The  results  have  been  divided  into  two 
groups;  those  for  weak  electron  emission  by  the  electrode  and  those  for 
strong  emission  by  the  electrode. 

An  electrode  is  termed  weakly-emitting  when  the  the  electric  field 
at  the  surface  is  negative  for  J  =  0,  and  the  emitted  electrons 
are  accelerated  into  the  plasma  (the  electrode  is  emission  limited) . 
Strongly-emitting  electrodes  are  space-charge  limited  for  J  = 0,  i.e,, 
the  electric  field  at  the  surface  is  positive  and  most  of  the  emitted 
electrons  are  returned  to  the  electrode. 

As  with  the  isothermal  results,  for  the  purpose  of  discussion  it  is 
convenient  to  use  a  moderate  value  of  e,  e.g.  e  =  Xd/Ir  =  10“^,  so  that  the 
structure  of  the  sheath  and  the  ionization  non-  equilibrixjm  region  can 
both  be  resolved  on  a  linear  scale.  All  solutions  are  for  the  typical 
value  of  fi  =  1/300,  but  as  before,  the  solutions  are  quite  insensitive  to 
fi  when  fi  1 .  A  value  of  Cj  =  16.7  has  been  chosen,  which  is  consistent  with 
a  potassium-seeded  plasma  at  3000  K. 

In  addition,  we  consider  only  anodes  {J  >  0)  in  the  results  presented 
here , 


Current-voltage  characteristics  are  plotted  in  Fig.  2.12  for  an  anode 
with  fie  (0)  =  0.02  and  ri<  (0)  =  0  . 

The  solid  lines  are  solutions  for  5  x  10^  <  0  <  5  x  lO"*  and  the 
dotted  line  is  the  solution  to  the  isothermal  model  (neglecting  Joule 
dissipation)  .  From  this  figure,  we  see  that  as  0  — »  oo,  the  isothermal 
result  is  reproduced  as  expected.  However,  for  smaller  values  of  0  (less 
collisional  coupling)  ,  Joule  heating  produces  a  significant  electron 


Fig.  2.12.  Current-voltage  characteristics  for  a  weakly -emitting 
anode  (ne(0)  =0.02,  nj(0)  =0)  with  e-10“^,  £<  =  16.7 

and  1/300.  The  solid  lines  are  solutions  for 
5  X  10^  <  0  <  5  X  10^  and  the  dotted  line  is  a  solution 
for  the  isothermal  model  (0  — ►  oo)  .  Vertical  arrows 
indicate  the  point  where  the  incremental  resistance  of 
the  plasma-electrode  boundary  layer  equals  zero,  i.e., 
dA^total/^J  =  0  . 


temperature  elevation  at  large  currents.  The  resulting  increase  in  the 
net  generation  rate  of  electron-ion  pairs  decreases  the  incremental 
resistance  of  the  boundary  layer  as  the  current  increases,  where  we 


define 


dA^iotal  Incremental  Resistance  of  the 
cf  J  Boundary  Layer 


(2.43) 


Eventually,  the  incremental  resistance  becomes  negative  and  the 
potential  drop  decreases  with  increasing  current.  We  interpret  this  to 
be  the  initial  stage  in  the  transition  from  a  diffuse  to  a  constricted 


mode  of  current  transfer.  For  each  curve,  a  vertical  arrow  indicates  the 
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point  where  the  incremental  resistance  of  the  plasma-electrode  boundary 
layer  equals  zero,  i.e.,  dA^totaifdJ^O. 

In  Fig.  2.13,  the  total  net  generation  rate  is  plotted  vs  J  for  the 
Scune  solutions  shown  in  Fig.  2.12.  The  solid  lines  are  solutions  for 
5  X  10^  <  O  <  5  X  10^  and  the  dotted  line  is  the  solution  for  the  isothermal 
model  (0  — »  oo)  . 
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Fig.  2.13.  Total  net  generation  rate  vs  7  for  a  weakly-emitting 
anode  (ne(0)  -0.02,  nj(0)  -0)  with  c«10~^,  =  16.7 

and  ~  1/300.  The  solid  lines  are  solutions  for 
5  X  10^  <  0  <  5  X  10^  and  the  dotted  line  is  a  solution 
for  the  isothermal  model  (0  — ♦  oo)  . 

For  a  weakly-emitting,  floating  electrode  (71^(0)  -ne(0)  -0,J-0),  an 
ion  sheath  forms  and  the  electron  temperature  in  the  sheath  is  decreased 
below  the  background  neutral  temperature  as  shown  in  Figs.  (2.14)  and 
(2.15).  When  the  electrode  is  emission  limited  (0  <  J  <  0.8,  in  this 
case) ,  the  electrons  move  against  the  field  which  decelerates  and  cools 
them.  This  also  decreases  the  total  net  generation  rate  to  a  value  below 
that  of  the  isothermal  solution.  The  ion  sheath  decreases  in  strength 
as  the  current  approaches  J  -  0.8  from  below.  For  J  -  0.8  the  ion  sheath 
essentially  disappears  (see  Fig.  2.16)  .  When  J  >  0.8  an  electron  sheath 
forms  adjacent  to  the  anode  surface  as  shown  in  Fig.  2.17.  The  electron 
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sheath  creates  a  positive  electric  field  near  the  anode  surface  and  the 
electron  temperature  is  raised  above  the  neutral  background  temperature. 
The  temperature  elevation  also  increases  the  total  net  generation  rate 
in  the  layer,  hence  for  J  >  0.8  here,  the  total  net  generation  rate  is 
larger  than  the  value  of  l/2  found  with  the  quasi-neutral  solution. 


Fig.  2.14.  Ion  and  electron  number  density  and  net  generation  rate 
spatial  profiles  for  fie(O)  -0.02,  fii(O)  -0,  c  =  10~\ 

Ci  -  16.7,  1/300,  ©  -  2  X  10^  and  7-0. 

C.3.2 _ Stronalv-Emittina  Anode 

Strong  electron  emission  by  the  anode  (^*(0)  >  0,^(i-0,J-0)  >  0) 
substantially  modifies  the  low  current  sheath  structure  discussed  above. 
The  emitted  electrons  neutralize  the  ion  sheath  which  would  norm  lly  be 
present  in  the  eibsence  of  emission.  An  electron  sheath  forms  adjacent 
to  the  anode  surface,  and  a  positive  electric  field  is  set  up  which 
accelerates  the  electrons  drawn  from  the  plasma,  heating  them.  The 
heated  electrons  raise  the  total  net  generation  rate  above  the  value 
calculated  with  the  isothermal  model. 


The  total  net  generation  rate  increases  with  increasing  current  and 
the  boundary  layer  incremental  resistance  decreases  with  increasing 
current.  Current-voltage  characteristics  are  shown  in  Fig.  2.18  for 


Fig.  2.15.  Electron  temperature  spatial  profile  for  ne(0)  ~  0.02, 
ni(0)-0,  «  -  10-^  Ci  -  16.7,  /I  =  1/300,  ©  =  2  X  10‘  and 


Fig.  2.16.  Ion  and  electron  number  density  and  net  generation  rate 
spatial  profiles  for  ne(0)  -0.02,  ^(0)  =0,  c=10 
Cj  -  16.7,  /i  =  1/300,  ©  -  2  X  10“  and  7=  0.8. 
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Fig.  2.17,  Ion  and  electron  number  density  and  net  generation  rate 
spatial  profiles  for  ne(0)  =0.02,  nj(0)  =0,  f=10~^, 
li  -  16.7,  n  =  1/300,  e  =  2  X  10^  and  7=  2.0. 

an  anode  with  ne(0)  =0.4.  Comparing  Figs.  (2.13)  and  (2.18),  it  is 
observed  that  for  a  given  value  of  the  collision  parameter  0,  increasing 
the  electron  emission  from  an  anode  decreases  the  current  at  which  the 
transition  from  a  diffuse  to  constricted  discharge  will  begin. 

These  results  demonstrate  that  Joule  dissipation  significantly 
enriches  the  possible  solutions  for  the  plasma-electrode  boundary  layer. 
As  with  the  simpler  models,  it  is  clear  that  the  total  net  generation 
is  the  key  variable  to  consider,  as  it  has  a  direct  bearing  on  the 
current-voltage  characteristics.  In  the  next  section,  we  investigate 
the  importance  of  the  remaining  (heretofore  neglected)  electron  energy 
equation  terms. 

C.4  Non-Isothermal  Solution.  Complete  Set  of  Equations  (Numerical) 

In  this  part,  we  replace  the  simplified  electron  energy  equation 
of  the  previous  section  with  a  more  complete  form  of  the  electron 
energy  equation  which  includes  additional  terms  for  electron  thermal 
conduction,  electron  energy  loss  due  to  the  net  generation  rate  of 
electron-ion  pairs,  and  energy  convected  by  the  electron  flux.  The 
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Fig.  2.18.  Current-voltage  characteristics  for  a  strongly- 
emitting  anode  (ne(0)  =  0.4,  ^<(0)  “  0)  with  c  =  10“^, 
Ci  =  16.7  and  /i  =  1/300.  The  solid  lines  are  solutions  for 
5  X  10^  <  0  <  5  X  10^  and  the  dotted  line  is  a  solution 
for  the  isothermal  model  (0  -+  oo)  . 


energy  equation  is  written: 


6  d  dTe\  5n  d  --  2n  - — 

5  dx  dx)^  \  +  ^dx  ^  1  +  M 


:~0n,  IT^-T)  -ur^n 
1  +  p  1  +  M 


This  second-order,  nonlinear  ordinary  differential  equation  replaces 
the  algebraic  energy  equation  of  the  previous  section,  so  two  additional 
boundary  conditions  are  required.  We  use  eq.  (2.42)  to  find  the  electron 
temperature  in  the  far-field  as  one  boundary  condition.  The  second 
boundary  condition  is  taken  at  the  electrode  surface,  where  we  assume 
that  for  an  electrode  with  strong  electron  emission,  TeiO)  —*  Tdectroder 
and  for  weak  electron  emission  dTe/dx\o  oc  E{0)  .  The  latter  condition  is 
based  on  a  nonequilibrium,  kinetic  analysis  by  Chung  [2.12]  of  the  near 
electrode  region.  This  form  of  the  energy  equation  also  has  a  singularity 
for  Te  when  n*  — ♦  0,  so  we  approximate  the  solution  for  a  non-emitting 


electrode  with  a  small  but  non-zero  value  for  n*  (0)  . 


Figure  2.19  shows  the  calculated  current-voltage  characteristics  for 
an  anode  with  strong  electron  emission,  and  compares  solutions  for  the 
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simple  and  complete  forms  of  the  electron  energy  equation.  Very  little 
difference  is  observed.  A  slightly  larger  potential  drop  is  calculated 
for  the  full  solution  for  currents  near  the  transition  to  a  constricted 
mode  discharge.  Equivalently,  the  solution  for  the  simple  form  of  the 
energy  equation  tends  to  slightly  underpredict  the  current  at  which 
constriction  begins. 


J 

Fig.  2.19.  Comparison  of  the  current -voltage  characteristics  for 
the  simple  electron  energy  equation  (solid  line)  and 
the  full  electron  energy  equation  (broken  line)  for  a 
strongly-emitting  anode. 

To  understand  why  the  current-voltage  characteristics  differ,  one 
must  examine  the  magnitude  of  each  of  the  terms  in  the  electron  energy 
equation.  As  shown  in  Fig.  2.20,  when  the  current  is  large  the  collisional 
loss  and  Joule  dissipation  are  in  balance  far  from  the  electrode.  Near 
the  electrode  surface,  the  electrons  gain  energy  by  Joule  heating  and 
convection,  and  lose  energy  by  conduction.  Note  that  in  this  figure, 
the  positive  sense  represents  an  energy  loss  by  the  electrons.  Joule 
dissipation  adds  energy  to  the  electrons,  hence  it  is  negative. 


The  additional  energy  loss  by  thermal  conduction,  not  included  in  the 
simple  model,  tends  to  decrease  the  electron  temperature  in  the  sheath. 
For  a  given  current  and  electron  emission  capability,  the  net  generation 


rate  is  slightly  less  for  the  full  solution  than  for  the  simple  energy 
equation  solution.  This  is  why,  as  the  transition  to  a  constricted 
discharge  begins,  the  total  potential  drop  for  the  full  solution  falls 
off  more  slowly  than  the  total  potential  drop  found  with  the  simple 
energy  equation. 
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Fig.  2.20.  Magnitude  of  the  electron  energy  equation  terms  for  a 
strongly-emitting  anode  at  high  current . 

The  mechanisms  for  energy  exchange  are  quite  different  for  small 
currents,  as  shown  in  Fig.  2.21,  for  a  floating  electrode  with  strong 
electron  emission.  Here,  the  ionization  and  collision  terms  are  dominant 
over  most  of  the  boundary  layer,  except  in  the  sheath.  In  the  sheath,  the 
energy  transfer  by  collisions  is  replaced  by  thermal  conduction,  as  was 
the  case  for  large  currents.  Also  note  that,  in  this  case,  the  electrons 
move  against  the  field  and  are  cooled  0),  so  that  the  total  net 

generation  rate  is  actually  decreased  below  the  quasi-neutral  value. 

To  summarize,  there  is  no  significant  difference  between  the 
isothermal  and  non-isothermal  results  for  an  anode  at  low  current,  but 
for  large  currents  the  increased  net  generation  rate  due  to  Joule  heating 
removes  the  current  saturation  behavior  found  with  the  isothermal 
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Currents  and  Fluid  Mechanics 


MHD  Secondary  Floy 

2.3.1  Introduction 

In  several  space  power  and  propulsion  systems  of  potential  interest 
to  the  Air  Force,  such  as  MPD  thrusters,  MHD  generators,  and  rail  guns, 
the  interaction  of  current-induced  forces  and  energy  transfer  with  the 
plasma  flow  can  be  critical  to  both  the  understanding  and  performance  of 
the  system.  The  channel  flow  of  a  plasma  through  a  transverse  magnetic 
field  is  influenced  by  the  Hall  effect,  which  causes  nonuniform  Lorentz 
forces  to  arise  in  the  direction  mutually  perpendicular  to  the  channel 
axis  and  the  magnetic  field.  The  forces  cause  secondary  flows  to 
develop  as  the  flow  progresses  down  the  region  of  electromagnetic-fluid 
interaction.  The  theoretical  treatment  of  this  problem  has  included  the 
work  of  Fay  [3.1],  Broer  et  al.  [3.2],  Sato  [3.3],  Tani  [3.4],  Apollon- 
skii  and  Kos'kin  [3.5],  Sastry  and  Bhadram  [3.6],  Liu  et  al.  [3.7],  and 
Ishikawa  and  Umoto  [3.8];  in  particular,  for  conditions  which  would  be 
typical  of  a  large-scale  MHD  power  generator  these  secondary  flows  were 
predicted  to  be  of  substantial  magnitude  in  the  calculations  reported  by 
Maxwell  et  al.  [3.9],  Bityurin  et  al.  [3.10],  and  Doss  and  Ahluwalia 
[3.11].  The  calculated  magnitude  of  these  flows  for  a  large-scale  de¬ 
vice  ranges  from  10  to  30%  of  the  bulk  velocity,  depending  on  conditions, 
and  the  predicted  consequences  include  the  development  of  significant 
asymmetries  in  the  profiles  of  mean  axial  velocity  and  of  temperature, 
with  resultant  effects  on  wall  heat  transfer  and  electrode  voltage 
drops.  In  addition,  Demetriades  et  al.  [3.12]  predicted  that  these 
secondary  flows  could  lead  to  magneto-aerothermal  Instabilities,  which 
could  include  boundary  layer  separation  and/or  electric  field  breakdown 
near  the  electrodes. 

Indirect  evidence  for  the  existence  of  MHD-induced  secondary  flow, 
reviewed  by  Girshick  and  Kruger  [3.13],  Included  measurements  of  temper¬ 
ature  profiles  obtained  by  James  and  Kruger  [3.14]  and  of  electrode 
voltage  drops  obtained  by  Barton  [3.15],  which  indicated  that  MHD 
effects  caused  the  anode  (top)  wall  boundary  layer  of  the  channel  to  be 
colder  than  expected  while  the  cathode  (bottom)  wall  was  hotter  than 
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expected.  It  was  hypothesized  that  these  results  could  be  ascribed  to  a 
secondary  flow  mechanism.  Demetriades  et  al.  [3.12]  conjectured  that 
evident  arc  deimage  to  the  center  of  the  amode  wall  in  the  downstream 
half  of  the  channel  used  in  the  High  Performance  Demonstration  Experi¬ 
ment  at  the  Arnold  Engineering  Development  Center  [3.16]  was  a  manifes¬ 
tation  of  the  magneto-aerothermal  instability  mentioned  above.  McClaine 
et  al.  [3.17]  reported  axial  velocity  profiles  measured  using  a  travers¬ 
ing  stagnation  pressure  probe  mounted  in  the  diffuser  of  the  Avco  Mk  VI 
generator.  The  profile  measured  between  the  top  and  bottom  walls  was 
found  to  be  strongly  skewed  by  MHD  forces,  presumably  as  a  consequence 
of  secondary  flow. 

The  work  reported  here  consisted  of  experiments  which  included  the 
first  direct,  quantitative  measurements  of  the  secondary  flow  field  in 
an  MHD  channel.  Additionally,  the  effect  of  secondary  flow  on  the 
plasma  momentum,  thermal  and  electrical  behavior  was  studied  by  making 
measurements  of  the  profiles  of  mean  axial  velocity,  turbulence  inten¬ 
sity  (axial  and  transverse)  and  plasma  voltage,  and  by  measuring  elec¬ 
trode  surface  temperatures. 

These  experiments  were  performed  under  conditions  which  were  not 
ideal,  in  that  the  flow  in  a  combustion-driven  MHD  channel  is  inherently 
less  "clean"  than  in  a  carefully  arranged  laboratory  flow  involving  low- 
temperature  air  or  water.  Factors  such  as  combustion  nonuniformities, 
cooled  walls,  and  irregularities  in  the  electrode-insulator  wall  struc¬ 
ture  inevitably  affect  the  flow.  Nevertheless,  the  present  experiments 
were  designed  to  minimize  extraneous  effects  to  the  extent  possible. 
Also,  the  channel  was  run  under  non-slagging  conditions,  and  the  mea¬ 
surements  were  made  under  different  orientations  of  the  magnetic  field, 
and  with  the  magnet  turned  off,  in  order  to  isolate  the  effects  on  the 
flow  field  of  Lorentz  forces  as  opposed  to  non-ideal  conditions. 

2.3.2  Research  Objectives 

The  goal  of  the  work  in  this  area  was  to  obtain  experimental  infor¬ 
mation  which  would  help  describe  the  three-dimensional  flow  field  in  an 
MHD  device.  Specifically,  the  major  objectives  which  guided  this  work 
were  the  following: 
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1. 


To  study  the  effect  of  electromagnetic  body  forces  on  the  secondary 
flow  field  including  investigation  of  the  magnitudes  of  these 
secondary  flows  and  how  the  vortex  structure  is  affected  by  applied 
magnetic  fields  and  discharges. 

2.  To  study  the  effect  of  MHD  interaction  on  mean  axial  velocity  pro¬ 
files. 

3.  To  study  the  effect  of  MHD  interaction  on  turbulence. 

2.3.3  Status  of  Research 

A.  Theoretical  Background 

A  schematic  representation  of  an  MHD  channel  is  pictured  jji  Fig.  3.1. 
The  channel  is  a  rectilinear  duct  whose  axis  points  in  the  x-direction. 
An  external  magnet  supplies  a  uniform  magnetic  induction  B,  pointing  in 
the  z-direction.  The  magnetic  Reynolds  n\amber  is  assumed  small,  thus 
the  induced  magnetic  field  is  negligible  in  comparison  to  the  applied 
field. 


Fig.  3.1.  Schematic  representation  of  an  MHD  channel. 


The  notion  of  charged  particles  through  the  magnetic  field  estab¬ 
lishes  a  uxB  Faraday  field  (in  the  laboratory  frame  of  the  reference) 
which  points  in  the  y-direction.  In  an  MHD  generator  this  field  is 
utilized  by  drawing  current  throu^  electrodes  placed  in  the  top  and 
bottom  walls,  these  electrodes  being  connected  through  an  external 
load.  In  addition  to  the  Faraday  current  Jy,  the  Hall  effect  for  a 
partially  ionized  gas  causes  the  current  density  vector  J  to  have  a 
local  component  in  the  x-direction.  This  can  be  expressed  by  consider¬ 
ing  the  conductivity  as  a  tensor  and  writing  Corn's  Law  in  the  form  (see, 
e.g.,  Mitchner  and  Kruger  [3.18]) 

J  =  OS'  ,  (3.1) 

where  E'  is  the  electric  field  in  a  frame  of  reference  moving  with  the 
mean  mass  velocity  of  the  fluid  and,  assuming  ion  slip  to  be  negligible, 

1  -R  0 

a  =  — ~ — r  BIO.  (3.2) 

1+6“^  0  0  1 

For  simplicity,  it  is  assumed  that  the  Hall  current  is  every¬ 
where  non-negative,  as  pictured  in  Fig.  3.2.  In  a  combustion  MHD  chan¬ 
nel  the  walls  must  be  cooled  because  of  neterials  constraints,  and 
consequently  the  plasma  is  much  more  conductive  in  the  core  than  near 
the  walls.  Furthermore,  the  uxB  interaction  goes  to  zero  at  the  walls 
because  of  the  no-slip  condition.  Thus,  to  first  order,  the  profile 
Jjj(z)  would  appear  approximately  as  indicated  in  Fig.  3.2. 

The  fluid  momentum  equation  under  the  stated  conditions  can  be 
written 

pg=-Vp+VT+JxB,  (3.3) 

where  p  is  the  mass  density,  p  the  thermodynamic  pressure  and  t  the  vis¬ 
cous  stress  tensor.  with  B  =  the  presence  of  a  non-zero  Hall  cur¬ 
rent  causes  the  JxB  Lorentz  force  to  have  a  coitponent  in  the  y-direction. 
Taking  the  curl  of  (3.3)  to  obtain  the  fluid  vorticity  equation,  the 
following  result  is  obtained  for  the  vorticity  Q; 
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transport  hot,  fast  core  plasma  toward  the  bottom  wall  while  sweeping 
relatively  cold,  slow,  side-wall  fluid  into  the  boundary  layer  of  the 
top  wall. 

As  the  axial  vorticity  and  the  distribution  of  are  coupled 
through  (3.5),  the  actual  structure  of  the  secondary  flow  field  could  be 
considerably  more  complicated  than  the  simple  pattern  suggested  in  Fig. 
3.2;  indeed,  might  not  increase  monotonically  from  the  walls  to  the 
core,  in  which  case  the  axial  vorticity  could  change  sign  in  some 
regions  of  the  cross-plane.  It  should  also  be  pointed  out  that 
commercial  Faraday  MHD  generators  would  suppress  the  flow  of  net  Hall 
current  by  using  segmented  electrodes  separated  by  insulators.  However, 
the  Hall  effect  would  still  cause  the  local  value  of  to  be  in  general 
non-zero;  only  the  integrated  value  of  over  the  cross-plane  would 
vanish,  and  then  only  in  the  absence  of  axial  leakage  currents.  A  typi¬ 
cal  profile  of  under  these  circumstances  and  of  the  resulting  six¬ 
cell  secondary  flow  field  are  shown  in  Fig.  3.3.  For  the  conditions  of 
envisioned  commercial-scale  MHD  generators  the  strength  of  these  second¬ 
ary  flows  could  be  significant,  in  spite  of  the  suppression  of  net  Hall 
cxjrrent;  this  case  was  discussed  by  Swean  et  al.  [3.19]  as  well  as  by 
several  of  the  authors  cited  earlier. 


Description  of  Experiments 


A. 1.1  Flow  Train  and  Run  Conditions 


MHD  secondary  flow  was  investigated  in  series  of  experiments  con¬ 
ducted  at  the  High  Temperature  Gasdynamics  Laboratory  at  Stanford  Uni¬ 
versity.  The  M-2  flow  train,  shown  in  Fig.  3.4  was  used.  This  con¬ 
sisted  of  a  combustor  with  a  nominal  2-MW  thermal  input  rating,  a  large 
mixing  plenum,  a  nozzle  through  which  the  combustion  products  were 
accelerated  into  a  duct  having  a  square  cross-section  measuring  5.1  cm 
on  a  side,  a  diffuser,  and  an  exhaust  system  Including  a  scrubber  and  a 
stack.  The  square  duct  consisted  of  a  run-in  section  and  the  active 
channel  itself,  which  was  enclosed  in  a  2.6-T  water-cooled,  copper-coil 
magnet.  Physical  dimensions  of  the  flow  train  are  given  in  Table  3.1, 
including  distances  from  the  nozzle  exit  to  the  measurement  locations . 


Table  3 . 1 


1 


Channel  Geometry 


Component 


Dimensions  ( in  cm) 


Combustor 

D 

= 

8.9,  L  = 

36.8 

Plenum 

H 

22.9,  W  = 

17.8 

:,  L  =  48.3 

Nozzle 

H 

= 

W  =  5.f>8, 

L  = 

2.5 

Run-in  section 

H 

= 

W  =  5.08, 

L  = 

74.8 

Active  channel 

H 

= 

W  =  5.08, 

L  = 

74.9 

Extender  section 

H 

= 

W  =  5.08, 

L  = 

20.3 

Diffuser 

H 

= 

13.3,  W  = 

6.4 

expanding  to 

Transfer  tulae 

D 

= 

25.4,  L  = 

290 

Vertical  exhaust/scrubber  system 

Distances  from  nozzle  exit  (cm): 

to  active  channel  entrance 
to  center  of  first  electrode 
to  conductivity-measurement  plane 
to  temperature -measurement  plane 
to  velocity-measurement  plane 
to  center  of  last  electrode 


50.8 

64.1 

98.4 

108.2 

108.9  (21.4  hydraulic  diam. ) 

109.9 


Electrodes:  width  =  5.08  cm,  length  =  1.81  cm. 


Run  conditions  are  listed  in  Table  3.2.  The  plasma  consisted  of 
the  combustion  products  of  ethyl  alcohol  burned  in  pure  oxygen.  The 
fuel  was  seeded  with  potassium  hydroxide.  The  total  flow  rate,  54.4  g/s, 
at  a  calculated  nozzle  exit  temperature  of  2750  K,  corresponded  to  a 
Reynolds  number  (based  on  hydraulic  diameter)  of  1.3x10^;  the  calculated 
nozzle  exit  velocity  was  171  itv/s.  This  value  has  been  used  throughout 
to  characterize  the  bulk  velocity  Uj,f  although  it  should  be  noted  that 
the  actual  bulk  velocity  decreases  somewhat  in  the  downstream  direction 
as  the  plasma  cools.  At  the  velocity  measurement  plane,  which  was 
located  21.4  hydraulic  diameters  downstream  of  the  nozzle  exit,  the 
boundary  layers  of  this  turbulent  flow  were  fully  developed  (or  at  least 
nearly  so)  in  the  absence  of  MHD  interaction. 
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Table  3.2 


Run  Conditions 


Fuel 

Ethyl  alcohol  (C-HcOH) 

N2/O2  ratio 

0 

Stoichiometry 

1.05  (fuel  rich) 

Seed 

Potassium  hydroxide 

(KOH) 

Potassium  mass  fraction  in  reactants 

2.1% 

Total  flow  rate 

54.4  g/s 

Reynolds  number  based  on  hydraulic  diameter 

1.3x10^ 

Calculated  core  conditions  at  nozzle  exit: 

Velocity 

~  170 

Tenperature 

2750  K 

Mass  density 

0.122  kg/m^ 

Mach  number 

0.17 

Electrical  conductivity 

14  S/m 

Electron  mobility 

Channel  conditions: 

0.37  m^/vs 

Wall  temperatures: 

Insulator  (estimated) 

~  1900  K 

Electrodes 

~  1000  K 

Magnetic  induction 

2.4  T 

Electron  Hall  parameter 

0.8 

Applied  Hall  current 

Magnetic  interaction  parameter  at 
velocity  measurement  plane. 

0,  3.6,  7.6,  9.3,  15 

.6 

based  on  applied  Hall  current 

0,  0.4,  0.85,  1.1,  1 

.74 

A.  1.2  Electrical  Configuration 

An  important  goal  of  the  experiments  was  to  isolate  the  effect  on 
the  fluid  flow  of  the  Hall  current,  since  the  Faraday  current  itself  has 
important  effects  because  it  induces  an  axial  Lorentz  force  given  by 
JyB.  It  was  also  considered  desirable  to  establish  conditions  which 
allowed,  insofar  as  possible,  for  a  straightforward  interpretation  of 
the  measured  secondary  flow  field.  Accordingly,  the  channel  was  con¬ 
nected  in  the  "Hall  configuration",  with  a  pair  of  electrodes  at  the 
upstream  end  connected  through  an  external  load  to  a  pair  of  electrodes 
at  the  downstream  end,  as  shown  in  Fig.  3.5.  A  controlled  current  was 
driven  through  the  plasma  by  means  of  a  bank  of  batteries,  and  by  varying 
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Fig.  3.5.  Electrical  configuration  used  for  the  min  experiments. 


the  load  resistance.  Thus  the  experimental  situation  corresponded  to 
the  case  illustrated  in  Fig.  3.2,  with  a  positive  net  Hall  current. 

A  non-dimensional  parameter  which  indicates  the  strength  of  the 
electro-mgnetic-f  luid  interaction  is  the  Stuart  number,  or  mgnetic 
interaction  parameter,  S.  For  this  configuration  an  appropriate  defini¬ 
tion  of  S  is 


(3.6) 


where  is  the  applied  Hall  current,  L  the  distance  from  the  first 

loaded  electrode  to  the  measurement  plane,  m  the  mass  flow  rate  and  U 
the  mean  axial  velocity.  The  interaction  parameter  in  this  form  repre¬ 
sents  the  ratio  of  the  Lorentz  twisting  force  to  the  axial  inertia  of  the 
flow.  The  interaction  parameter  values  used  are  listed  in  Table  3.2. 

A. 1.3  Laser-Doppler  Anemometry  System 

Plasma  velocities  were  measured  using  laser  Doppler  anemometry. 
The  3JDA  system,  shown  schematically  in  Fig.  3.6,  used  a  dual-beam, 
single-color,  backscatter  configuration.  The  optical  train  consisted 
primarily  of  TSI  modular  components.  Measurements  of  the  x-  and 
y-directed  velocity  components  were  nade  separately  by  rotating  com¬ 
ponents  of  the  optical  train  by  90®.  For  the  y-directed  measurements  a 
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Fig.  3.6.  Schematic  of  laser-Doppler  anemometry  optical  train. 

Bragg  cell  was  used  to  frequency  shift  one  of  the  tr'"smitted  beams  by 
40  MHz  so  as  to  allow  discrimination  of  positive  and  negative  veloci¬ 
ties.  Tight  spatial  and  spectral  filtering  was  used  to  maximize  the 
signal -to -noise  ratio  of  the  collected  scattered  light;  this  included  a 
narrow-line  (10-A  half-width)  interference  filter,  the  purpose  of  which 
was  to  filter  the  highly  radiative  background  emanating  from  the  hot, 
potassium-laden  plasma. 

Suitable  light-scattering  centers  were  obtained  by  seeding  the  fuel 
with  zirconia  particles  having  a  mean  diameter  of  0.7  ym.  These  were 
mixed  into  an  alcohol -glycerine  slurry  and  then  injected  into  the  fuel 
line  by  means  of  a  mechanical  injector  utilizing  a  slowly-descending 
piston.  Optical  access  to  the  plasne  was  provided  by  an  open  slit  cut 
into  the  channel  side  wall  and  located  at  the  leading  edge  of  the  down¬ 
stream  electrodes. 


The  photomultiplier  tube  signal  was  processed  by  a  TSI  1990A 
counter-type  processor  interfaced  to  a  laboratory  computer.  Typical 
data  rates  were  ~500  Hz;  total  sample  size  was  chosen  in  each  case  so 
that  the  sampling  interval  was  about  one  minute. 

The  major  sources  of  experimental  uncertainty  in  these  measurements 
were  the  imperfect  stationarity  in  conditions  ( including  combustion  tem¬ 
perature  and  applied  current)  and  the  difficulty  in  perfectly  aligning 
the  intersecting  beams  for  the  transverse  velocity  measurements.  The 
latter  problem  is  always  important  in  secondary  flow  measurements, 
because  a  slight  misalignment  introduces  a  portion  of  the  much  higher 
axial  velocity  into  the  transverse  velocity  measurement  (see  [3.20]  for 
a  brief  review  of  this  question  in  relation  to  secondary  flow  measure¬ 
ments  made  by  hot-wire  anemometry  as  well  as  by  LDA) .  In  the  present 
experiments  the  problem  was  exacerbated  by  the  fact  that  turning  the 
magnet  on  or  off  caused  a  slight  shift  of  the  entire  channel;  a  calibra¬ 
tion  procedure  was  devised  to  correct  for  the  effect  of  this  occurrence 
on  the  measured  transverse  velocities.  A  detailed  uncertainty  analysis 
by  Girshich  [3.21]  estimated  the  cumulative  uncertainty  (to  10:1  confi¬ 
dence)  in  the  mean  velocity  measurements  as  ±  0.016  for  v/Uj^  in  the  MHD 
cases  and  ±  0.0055  in  the  non-MHD  case,  and  ±  0.045  for  u/Uj^  ,  where  v 
and  u  are,  respectively,  the  y-  and  x-directed  components  of  the  mean 
velocity.  The  uncertainty  in  spatial  location  was  estimated  as  ±  0.11 
for  y/H,  ±  0.031  for  z/W  in  the  measurements  of  secondary  flow  and 
±  0.042  for  z/W  in  the  measurements  of  axial  velocity,  where  H  and  W 
are,  respectively,  the  channel  height  and  width.  The  major  contribu¬ 
tions  to  this  uncertainty  were  the  finite  size  of  the  measurement  volume 
and  the  uncertainty  in  locating  the  wall. 


A. 1.4  Conductivity  Measurements 


The  effect  of  secondary  flow  on  the  plasma  conductivity  profile  was 
investigated  by  driving  a  Faraday  discharge  across  a  set  of  electrodes 
whose  external  circuitry  was  separate  from  that  used  to  drive  the  Hall 
current.  This  set  consisted  of  three  electrode  pairs  located  just  up- 
stre£un  of  the  optical  port.  These  electrodes  were  used  only  for  the 
conductivity  measurements,  and  were  open-circuited  at  all  other  times. 


Near-electrode  voltage  drops  in  both  the  x-  and  y-directions  were  ob¬ 
tained  fran  a  grid  of  iridium  voltage  pins  inserted  through  a  channel 
side  wall.  Voltage  distributions  were  obtained  over  a  range  of  values 
for  Jy/  and  under  various  applied  conditions,  including: 

(1)  applied  Faraday  current  only; 

( 2 )  both  Faraday  current  and  Hall  current ; 

(3)  Faraday  current  with  magnetic  field;  and 

(4)  Faraday  current.  Hall  current,  and  magnetic  field. 

For  cases  ( 3 )  and  ( 4 ) ,  measurements  were  made  with  the  magnetic 
field  pointing  in  both  the  positive-z  and  negative-z  directions.  The 
open-circuit  voltage  with  magnetic  field  was  measured  as  well.  Note 
that  only  case  (4)  would  be  expected  to  drive  secondary  flow. 

The  voltage  information  was  analyzed  using  an  approximate  model 
based  on  the  expression,  obtainable  from  the  generalized  Ohm's  law, 


E  -  uB  +  3E 


where  p  is  the  resistivity,  Ey  and  the  electric  field  components,  and 
u  the  axial  velocity.  The  overbar  indicates  an  average  over  a  region  of 
the  cross-plane  which  encompasses  the  width  W  in  the  z-direction  and 
ranges  in  the  y-direction  from  the  electrode  surface  to  the  height  of 
the  nearest  voltage  pin,  a  distance  of  6.2  mm,  or  H/8,  H  being  the  chan¬ 
nel  height.  The  major  assumptions  of  the  model  are  that  the  transverse 
conductivity  distribution  in  the  absence  of  secondary  flow  is  symmetric; 
that  the  near-electrode  voltage  gradients  are  much  greater  in  the  y- 
direction  th^u^  in  the  z-direction;  that  the  term  Jy  in  the  denominator 
of  (3.7)  is  constant  and  can  be  estimated  using  a  "current-constriction 
factor"  obtained  from  a  comparison  of  the  data  from  cases  (1)  and  (3)  in 
the  above  list;  and  that  6  is  constant.  The  data  from  the  first  case 
were  also  used  to  estimate  a  constant  value  for  the  cathode  surface- 
sheath  voltage  drop,  which  was  subtracted  from  the  measured  value  of  Ey 
for  the  near-cathode  data. 


All  electrodes  were  fitted  with  thermocouples  located  close  to 
their  surfaces.  Of  particular  interest  were  the  temperatures  of  the 
electrode  pair  just  upstream  of  the  velocity  measurement  plane,  for 
cases  in  which  these  electrodes  were  not  passing  current.  A  nuiidt>er  of 
measurements  were  recorded  in  which  a  case  without  secondary  flow  (i.e., 
without  the  simultaneous  presence  of  magnetic  field  and  Hall  current) 
was  immediately  followed  or  preceded  by  a  case  with  secondary  flow 
(i.e.,  both  were  present).  The  change  in  surface  temperature  between 
these  two  cases  could  then  be  attributed  to  altered  heat  transfer  as  a 
result  of  secondary  flow.  As  with  the  other  measurements,  data  were 
collected  for  both  the  positive-z  and  negative-z  orientations  of  the 
magnetic  induction. 


sriaental  Results 


Transverse  Velocity  Measurements 


B.1.1  Won-MHD  Case 

In  order  to  have  a  base  case  with  which  to  compare  the  measurements 
with  MHD  effects  present,  measurements  of  the  y-directed  velocity  were 
made  for  the  case  of  zero  magnetic  induction  and  zero  current.  The 
results  are  shown  in  Fig.  3.7.  The  peak  measured  transverse  velocity 
was  0.05  Ujj.  This  is  significantly  higher  than  reported  measurements 
for  secondary  flows  driven  purely  by  the  interaction  of  turbulence 
stresses  with  the  corners  of  a  rectangular  duct.  Several  factors  may  be 
involved  in  this  discrepancy,  including  combustion  nonuniformities  and 
other  nonidealities  related  to  the  combustor,  density  gradients,  and  the 
heterogeneous  nature  of  the  channel  walls.  (The  top  and  bottom  walls 
had  alternating  sections  of  brick  and  of  electrodes,  whereas  the  side 
walls  had  only  brick.) 


The  data  appear  to  indicate  a  top-bottom  asymmetry,  in  that  flow 
toward  the  corners  appears  at  the  bottom  of  the  cross-plane  but  not  at 
the  top.  This  may  possibly  be  ascribable  to  top-bottom  density  nonuni¬ 
formity  from  the  combustor.  It  is  also  possible  that  flow  toward  the  top 
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Fig.  3.7.  Measurements  of  the  mean  y-directed  velocity 
without  magnetic  field. 


corners  would  have  been  observed  had  measurements  been  made  closer  to 


these  corners. 


B.1.2  MHD  Case  with  S  =  1.1 


The  value  of  the  interaction  parameter  used  in  the  first  series  of 
experiments  was  1.1  at  the  plane  of  the  velocity  measurements.  This 
value  is  comparable  to  envisioned  values  in  large-scale  channels.  Mea¬ 
surements  of  the  y-directed  velocity  were  made  with  the  magnetic  field 
pointing  in  the  positive-z  and  negative-z  directions.  Reversing  the 
magnet  polarity  reverses  the  sign  of  the  Lorentz  force  and  so  should  be 
approximately  equivalent  to  rotating  the  channel  by  180®  about  its  axis. 


The  measurements  for  the  case  of  B  pointing  in  the  positive-z 
direction  are  shown  in  Fig.  3.8.  In  this  figure  and  in  the  figures  fol¬ 
lowing,  the  secondary  flow  appears  stronger  near  the  left  side  wall  than 
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Fig.  3.8.  Measurements  of  the  mean  y-directed  velocity 
with  the  magnetic  field  pointing  in  the 
positive-2  direction  (S  =  1.1). 


near  the  right;  it  is  believed  that  this  is  an  artifact  caused  by  a  con¬ 
sistently  erroneous  measurement  of  the  wall  location.  Shifting  all  the 
measurement  locations  to  the  left  by  an  eunount  &(z/W)  *=  0.035  would 
produce  a  symmetric  result,  which  is  physically  what  we  would  e;q>ect. 
This  systematic  error  is  additional  to  the  uncertainty  in  spatial  loca¬ 
tion  discussed  above. 

In  any  case,  the  results  suggest  the  presence  of  two  large,  counter¬ 
rotating  secondary  flow  cells,  with  the  flow  in  the  core  directed  down¬ 
ward  in  agreement  with  the  first-order  argument  Illustrated  in  Fig.  3.2. 
The  peeOc  measured  velocities,  found  in  the  side-wall  region,  measured 
0.12  Ujj;  as  the  magnitude  of  the  transverse  velocity  was  increasing 
steeply  toward  the  side  wall,  it  is  reasonable  to  assume  that  signifi¬ 
cantly  higher  velocities  would  have  been  found  had  measurements  been 
made  closer  to  the  walls. 
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The  vorticity  is  concentrated  in  the  top  half  of  the  channel,  with 
the  vortex  cells  closing  at  about  y/H  =  0.35.  An  explanation  for  this 
behavior  is  suggested  by  a  three-dimensional,  numerical  simulation  per¬ 
formed  by  Maxwell  et  al.  [3.22]  for  conditions  which  were  simileur  to 
those  of  the  present  experiments.  In  these  calculations,  the  secondary 
flow  field  at  the  upstream  end  of  the  interaction  region  has  a  symmetric 
structure  similar  to  that  in  Fig.  3.2.  As  the  flow  progresses,  hot 
plasma  from  the  core  has  time  to  be  swept  sideways  and  then  upwards  as 
it  follows  a  secondary  flow  path.  This  has  the  effect  of  increasing 
3Jjj/3z  in  the  upper  side-wall  regions,  which  in  turn  intensifies  the 
vorticity  in  these  regions.  The  coupling  between  and  axial  vor¬ 

ticity  causes  the  secondary  flow  cells  to  migrate  toward  the  top  walls 
in  the  downstream  half  of  the  channel.  In  this  respect,  the  simulation 
and  the  measurements  are  in  good  agreement.  The  quantitative  agreement 
is  also  reasonably  good,  although  precise  comparisons  would  require  a 
more  detailed  specification  than  is  available  of  the  side-wall  tempera¬ 
ture  distributions  and  of  the  upstream  turbulence  history.  Both  of 
these  factors  were  shown  in  the  simulation  to  have  a  substantial  influ¬ 
ence  on  the  secondary  flow  field:  the  side-wall  temperature  because  it 
affects  9Jjj/32;  the  turbulence  intensity  because  of  its  cross-plane 
smoothing  effect. 

The  measurements  in  the  region  at  the  bottom  of  the  cross-plane 
indicate  a  flow  which  was  everywhere  directed  away  from  the  bottom  wall. 
A  possible  explanation  for  these  measurements  is  discussed  in  connection 
with  the  temperature  results. 

In  Fig.  3.9,  the  results  for  the  case  with  the  magnet  polarity 
reversed  are  shown.  As  anticipated,  the  velocity  field  was  upside-down 
compared  to  the  results  in  Fig.  3.8.  A  few  measurements  were  made  in 
this  case  closer  to  the  side  walls,  resulting  in  a  peak  transverse- 
velocity  measurement  of  0.15 
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Fig.  3.9.  Measurements  of  the  mean  v-directed  velocity 
with  the  magnetic  field  pointing  in  the 
negative-z  direction  (S  =  1.1). 

B.2  Effects  of  Secondary  Flow 

The  measured  secondary  flow  field  was  certainly  strong  enough  to 
have  an  effect  on  several  aspects  of  the  plasma  momentum  and  thermal 
behavior.  This  section  reports  the  results  of  the  measurements  of  mean 
axial  velocity,  turbulence  intensity,  electrical  conductivity,  and 
electrode  surface  temperature. 

B.2.1  Mean  Axial  Velocity 

Measurements  of  the  mean  axial  velocity  with  the  magnetic  field 
pointing  in  the  positive -z  direction  are  shown  in  Fig.  3.10;  these 
results  correspond  to  the  secondary  flow  field  shown  in  Fig.  3.8.  The 
distoicion  of  the  velocity  profile  was  dramatic:  the  peak  velocity  was 
forced  downward  from  the  center  of  the  cross-plane  to  a  position  of  about 
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Fig.  3.10.  Contour  map  of  the  measured  mean  x-directed  velocity 
(normalized  by  Uj^)  with  the  magnetic  field  pointing 
in  the  positive-z  direction  (S  =  1.1). 
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y/H  =  0.3.  The  constant  velocity  lines  near  the  top  wall  follow  the  i 

vortex  pattern,  as  slow-moving  fluid  was  swept  upward  along  the  side  | 

wall  and  then  pushed  downward  along  the  z-centerline.  The  axial  veloc¬ 
ity  profiles  along  the  z-centerline  for  three  conditions  of  the  magnetic 
field  are  shown  in  Fig.  3.11.  As  expected,  the  profile  was  symmetric 
when  Lorentz  forces  were  absent;  with  Lorentz  forces  present,  the 
profile  became  strongly  skewed  according  to  the  direction  of  the  force. 


These  results  are  similar  to  the  measurements  reported  by  McClaine 
et  al.  [3.17].  In  their  experiments,  performed  in  the  larger  Avco  Mk  VI 
channel,  the  conditions  included  a  magnetic  induction  of  4T,  an  interac¬ 
tion  length  of  2.5  m,  and  peak  velocities  of  about  1000  m/s.  Unfortu¬ 
nately,  the  electrical  loading  which  was  used  during  these  particular 
measurements  was  not  reported;  their  tests  encompassed  a  wide  range  of 
loading  conditions  in  both  Faraday  and  diagonal  configurations. 
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Fig.  3.11. 
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Measurements  of  the  mean  x-directed  velocity  along  the 
z-centerline ,  for  three  conditions  of  the  transverse 


Lorentz  force  F  { S  =  1.1). 


While  the  axial  velocity  profiles  measured  in  the  present  experi¬ 
ments  were  strongly  skewed,  there  was  no  evidence  of  the  occurrence  of  a 
magneto-aerothermal  instabilitiy.  This  is  not  remarkable  in  view  of  the 
fact  that  the  experiments  were  designed  to  suppress  the  Faraday  current 
while  Imposing  a  Hall  current.  In  the  magneto-aerothermal  case,  the 
Hall  cxirrent  is  driven  by  the  conductivity  nonuniformity  encountered  by 
the  Faraday  current;  this  nonuniformity  is  then  amplified  by  secondary 
flows  driven  by  the  Hall  current  [3.12].  In  the  present  experiments, 
the  local  Hall  current  is  still  coupled  to  secondary  flow,  but  the 
regulation  of  the  total  Hall  current  may  have  precluded  the  emergence  of 
an  instability. 
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B.2.2  Tarbulence  Intensity 

Previous  measurements  of  turbulence  intensity  in  the  Stanford  M-2 
channel  have  found  relatively  high  values  of  the  center-point  axial  tur- 
bulence  intensity  u^/u^,  where  u^  is  the  turbulence  velocity  (the  stan¬ 
dard  deviation  of  the  velocity-probability-distribution  function)  and  u^ 
is  the  center-point  mean  velocity.  For  example,  Reis  and  Kruger  [3.23] 
measured  values  ranging  from  6.6%  to  9.6%  for  a  variety  of  non-MHD  flow 

cases;  in  the  present  experiments,  several  measurements  of  u  /u  yielded 

c  c 

values  ranging  from  5.8%  to  6.3%.  For  comparison.  Mailing  and  Whitelaw 
[3.20]  in  their  water-flow  experiment  measured  values  of  4-5%  at  an 
axial  location  where  the  boundary  layers  had  merged.  It  has  been 
conjectured  that  the  relatively  high  axial  turbulence  in  the  Stanford 
channel  is  caused  by  combustion  nonuniformities.  The  measured  value  in 
the  non-MHD  case  of  the  y-directed  core  turbulence  v,  again  normalized 
by  Uj,,  was  3.75%;  this  is  virtually  identical  to  the  value  reported  by 
Melling  and  Whitelaw. 

Measurements  of  the  y-directed  turbulence  intensity  in  the  non-MHD 
case  are  shown  in  Fig.  3.12,  in  which  the  measured  values  of  v  have  been 
normalized  by  Uj^.  As  expected,  the  distribution  was  symmetric  and  |v| 
was  greater  near  walls  parallel  to  y  than  near  walls  normal  to  y. 

The  measurements  made  in  the  MHD  case,  with  the  magnetic  field 
pointing  in  the  positive-z  direction,  are  shown  in  Fig.  3.13.  The 
region  of  minimum  turbulence  (whose  magnitude  was  somewhat  less  than  in 
the  non-MHD  case  because  of  turbulence  damping  by  the  magnetic  field; 
cf .  [3.23] )  was  shifted  downward  by  Lorentz  forces. 

The  measurements  for  u/U,  ,  shown  in  Fig.  3.14,  exhibited  the  saune 

b 

effect;  these  same  measurements  are  shown  in  Fig.  3.15  normalized  in¬ 
stead  by  the  local  mean  velocity  u.  In  these  measurements,  there 


appears  to  be  a  region  beneath  the  center  of  the  top  wall  where  the 
turbulence  was  greater  than  its  value  closer  to  the  wall;  this  might  be 
explainable  in  terms  of  secondary-flow  convection. 


B.2.3  Electrical  Conductivity 


The  results  of  the  conductivity  measurements  for  the  MHD  cases  are 
summarized  in  Fig.  3.16,  in  which  the  near-electrode  resistivity  defined 
by  eq.  (3.’’)  has  been  normalized  by  the  value  of  the  core  conductivity 
measured  without  magnetic  field.  The  term  "anode  region"  and  "cathode 
region"  refer,  respectively,  to  the  regions  adjacent  to  the  top  and 
bottom  electrodes  which  were  used  for  the  conductivity  measurements. 
The  abscissa  in  this  figure  represents  the  value  of  the  core  current 
density  associated  with  the  Faraday  discharge;  for  all  the  data  shown, 
the  Hall  circuit  was  active  as  well,  corresponding  to  Case  (4)  described 
earlier. 

The  decrease  in  resistivity  as  increased  was  an  expected  conse¬ 
quence  of  Joule  heating.  Otherwise,  these  results  indicate  that  second¬ 
ary  flow  caused  the  resistivity  distribution  to  become  significantly 
asymmetric.  The  boundary  layer  receiving  cold  side-wall  fluid  was  much 
more  resistive  than  the  boundary  layer  toward  which  hot  core  fluid  was 
directed;  the  average  value  of  the  ratio  of  the  higher  resistivity  to 
the  value  at  the  opposite  wall  was  1.8.  The  importance  of  these  results 
is  that  near-electrode  voltage  drops,  which  are  a  performance-limiting 
effect  in  combustion  MHD  devices,  can  be  significantly  altered  by  sec¬ 
ondary  flow. 

B.2.4  Electrode  Surface  Temperatures 

Thirteen  separate  measurements  were  made  of  the  effect  of  secondary 
flow  on  electrode  surface  temperatures,  as  described  earlier.  The 
results  are  shown  in  Fig.  3.17.  Each  point  in  this  figure  represents 
the  temperatxire  with  secondary  flow  (with  simultaneous  Hall  current  and 
magnetic  field)  minus  the  temperature  without  secondary  flow  (without 
current);  the  latter  temperature  was  typically  about  1000  K.  A  separate 
analysis  of  data  not  shown  in  the  figure  indicates  that  the  presence  of 
magnetic  field  alone  made  no  difference  in  the  temperature,  while  the 
presence  of  Hall  current  alone  caused  a  warming  of  both  the  top  and 
bottom  electrodes  by  an  average  of  about  10  K,  because  of  Joule  heating. 


The  effect  of  secondary  flow  in  the  positive-B  case  is  seen  to  be  a 
substantial  heating  of  the  bottom  electrode  and  a  smaller  but  still  sig¬ 
nificant  cooling  of  the  top  electrode — reversing  the  magnet  polarity 
reversed  the  results.  The  electrode  toward  which  hot  core  fluid  was 
pushed  was  heated  on  the  average  by  B5  K;  the  opposite  electrode  was 
cooled  by  an  average  of  32  K.  Using  the  fact  that  the  electrodes  were 
cooled  by  water  at  about  300  K,  it  appears  that  secondary  flow  caused  a 
12%  increase  in  the  heat  flux  from  the  plasma  to  the  heated  electrode 
and  a  5%  reduction  in  the  heat  flux  to  the  cooled  electrode. 

Aside  from  the  obvious  conclusion  that  secondary  flow  has  the 
expected  qualitative  asymmetric  effect  on  heat  transfer  to  the  top  and 
bottom  walls,  it  is  interesting  to  consider  the  surprising  result  that 
the  effect  at  the  heated  electrode  was  more  than  twice  as  great  as  the 
effect  at  the  cooled  electrode.  This  result  is  surprising,  first, 
because  the  secondary  flow  cells  appeared  to  be  concentrated  near  the 
cooled  electrode,  and,  second,  because  the  resistivity  measurements 
showed  the  opposite  trend,  with  a  stronger  effect  being  observed  in  the 
region  adjacent  to  the  cooled  electrode. 

A  possible  explanation  is  as  follows.  As  already  noted,  the  sec¬ 
ondary  flow  cells  not  only  sweep  cold  side-wall  fluid  into  the  top-wall 
boundary  layer,  but  also  eventually  convect  hot  core  fluid  around  a  loop 
into  the  same  boundary  layer  in  the  side-wall  region.  As  the  stainless- 
steel  electrodes  are  thermally  quite  conductive,  the  net  result  is  that 
the  cooling  effect  is  mitigated;  no  such  effect  occurs  for  the  boundary 
layer  of  the  bottom  wall,  toward  which  core  fluid  has  been  pushed  but 
which  has  less  vortical  mixing.  Distortion  of  the  turbulence  distribu¬ 
tion  could  also  be  a  factor  in  the  asymmetrical  heat  transfer.  Near¬ 
electrode  voltage  drops,  on  the  other  hand,  are  dependent  on  complex 
physical  processes  which  are  strongly  but  nonlinearly  temperature- 
dependent.  Since  the  calculated  resistivities  were  averaged  over  an 
area  described  earlier,  it  is  not  remarkable  that  the  net  quantitative 
effects  for  the  resistivity  distribution  could  be  quite  different  from 
that  for  heat  transfer,  although  the  qualitative  trends  were  the  same. 


The  temperature  results  suggest  a  possible  explanation  for  the 
secondary  flow  measurements  In  the  region  near  the  bottom  wall  (In  the 
posltlve-B  case).  The  substantial  warming  of  the  bottom  electrode  sug¬ 
gests  that  the  mass  density  near  the  bottom  wall  may  have  been  decreas¬ 
ing  In  the  axial  direction,  causing  a  net  reduction  In  the  axial  mass 
flux  which  would  have  to  be  balanced  by  a  flow  away  from  the  wall . 
Flow  away  from  the  bottom  wall  did  not  appear  In  the  simulation  of 
Maxwell  et  al.  [3.22].  However,  their  calculations  assumed  a  constant 
wall  ten^rature.  It  Is  possible  that  such  a  flow  effect  might  appear 
In  the  simulations  If  the  wall  temperatxire  were  coupled  to  the  flow 
field. 


B.3  Dependence  of  Secondary  Flow  and  Related  Effects  on  Interaction 

Parameter 

Since  secondary -flow  phenomena  are  highly  complex  and  nonlinear,  It 
Is  Important  to  determine  their  dependence  on  driving  parameters.  The 
simplest  theoretical  description  of  secondary  flow  suggests  that  the 
magnitude  of  the  secondary -flow  velocities  should  be  directly  propor¬ 
tional  to  the  interaction  parameter,  previously  defined  in  eq.  (3.6). 
Experiments  were  conducted  to  determine  the  actual  dependence,  and  the 
results  are  discussed  here. 

A  straightforward,  although  indirect,  measure  of  the  magnitude  of 
the  secondary  flow  is  the  displacement  of  the  axial  velocity  profile — as 
in  Fig.  3.11.  Figure  3.18  shows  the  effect  of  the  interaction  parameter 
on  this  displacement  of  the  axial  velocity  for  three  values  of  inter¬ 
action  parameter,  S  =  0.1,  0.4,  and  0.85.  It  is  notable  that  the  dis¬ 
placement  is  significant  even  for  the  modest  interaction  S  =  0.4. 
Figure  3.19  shows  the  corresponding  effect  on  turbulence  intensity  for 
the  same  Interaction  parameter  values.  It  Is  evident  from  these  data 
that  the  secondary  flow  does  scale  with  the  interaction  parameter,  at 
least  for  interaction-parameter  values  less  than  approximately  S  =  1 . 
For  interaction  parameter  values  above  this  level,  however,  the  magni¬ 
tude  of  the  secondary  flow  does  not  appear  to  change  markedly.  This  is 
seen  in  Fig.  3.20,  in  which  the  axial  profiles  have  been  replotted  with 
the  additional  case  S  =  1.74  included,  and  in  Figs.  3.21  and  3.22,  which 
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Fig.  3.22  Measurements  of  the  mean  y-directed  velocity  (normalized  by 
U^)  for  interaction  parameter  S  =  1.74. 

show  the  y-directed  velocity  mapped  over  the  cross  plane  for  the  two 
cases  S  =  0.85  and  S  =  1.74,  respectively.  It  is  seen  that  the  second¬ 
ary  flow  appears  to  "saturate"  for  high  interaction  parameters.  This 
apparent  saturation  was  not  anticipated  prior  to  the  experiments  actu¬ 
ally  being  conducted.  This  effect  may  result  at  least  in  part  from  the 
fact  that  the  secondary  flow  causes  a  major  redistribution  in  the  flow 
field  which,  by  altering  the  conductivity  distribution,  couples  back  to 
the  driving  Lorentz  forces  in  a  nonlinear  way. 

B.4  Summary  and  Conclusions 

Experiments  were  conducted  to  measure  MHD-induced  secondary  flow. 
The  results  included  the  first  direct,  quantitative  measurements  of 
secondary  flow  in  a  combustion  MHD  channel.  The  experiments  were 
performed  in  a  laboratory  channel,  with  dimensions  an  order  of  magnitude 
smaller  than  in  envisioned  full-scale  MHD  devices;  however,  the  use  of 
an  applied  Hall  current  and  of  a  relatively  low  axial  velocity  allowed 
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values  for  the  magnetic  interaction  parameter  based  on  the  Hall  current 
comparable  to  anticipated  values  in  the  core  of  large-scale  MHD  chan¬ 
nels  . 

Laser-Doppler  anemometry  measurements  of  transverse  velocities 
under  various  controlled  conditions  showed  that  Lorentz  forces  caused 
large-scale  secondary  flows  to  develop  over  a  downstream  channel  cross¬ 
plane.  Peak  measured  transverse  velocities  were  as  high  as  15%  of  the 
bulk  velocity.  The  structure  of  the  measured  secondary  flow  field  was 
qualitatively  in  accord  with  simple  models,  with  a  concentration  of  the 
vortices  toward  the  top  wall  explainable  in  terms  of  the  coupling 
between  axial  vorticity  and  nonuniformities  in  the  Hall  current  density. 

The  simplest  theoretical  description  of  secondary  flow  indicates 
that  the  magnitude  of  the  secondary  flow  velocities  should  be  directly 
proportional  to  the  interaction  parameter.  However,  the  experiments 
showed  that  secondary  flow  causes  a  major  redistribution  of  the  flow 
which  couples  back  to  the  driving  Lorentz  forces  in  a  nonlinear  way. 
The  magnitude  of  the  flow  appears  to  "saturate"  for  interaction  parame¬ 
ter  values  greater  than  approximately  S  =  1 . 

Several  effects  of  MHD  secondary  flow  were  also  investigated, 
including  the  skewing  of  the  cross-plane  distributions  of  mean  axial 
velocity  and  of  the  axial  and  transverse  components  of  turbulence 
intensity,  and  the  development  of  asymmetries  in  the  wall  heat  flux  and 
in  the  near-electrode  conductivity.  These  results  consistently  demon¬ 
strated  that  secondary  flow  exerted  a  profound  influence  on  the  plasma 
momentxim,  thermal,  and  electrical  behavior.  The  major  conclusion  is 
that  secondary  flow  is  an  important  effect  which  needs  to  be  accurately 
assessed  in  designing  full-scale  MHD  devices. 
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